Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,551 A→C 100% pseudogene (41/402 nt) ytfA → pseudogene, related to transcriptional regulators

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,417,5510AC99.6% 955.3 / ‑5.5 242pseudogene (41/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (109/132);  total (110/132)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 250 total aligned reads displayed.

GGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGA  >  CP009273/4417406‑4417694
                                                                                                                                                 |                                                                                                                                               
ggAGAAAATGACTTCCACGCCTTACACCGTACATTCATGTTGTTCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGAAGGTTATCCAGCCACAAAAACCCATTGAGAAAATTATTGATTTTAAATTGGctt                                                                                                                                               >  2:2395738/1‑148 (MQ=255)
 gagaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGctta                                                                                                                                              >  2:1237704/1‑148 (MQ=255)
   gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGcttatt                                                                                                                                            >  1:2218186/1‑148 (MQ=255)
    aaaaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCttatta                                                                                                                                           >  1:314537/1‑148 (MQ=255)
    aaaaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCttatta                                                                                                                                           <  1:2469047/148‑1 (MQ=255)
     aaaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                          <  2:1869692/148‑1 (MQ=255)
       aTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTatat                                                                                                                                        <  1:1846132/148‑1 (MQ=255)
          aCTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTCCATGGGCTTATTATATTGc                                                                                                                                     <  1:836146/148‑1 (MQ=255)
            ttCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTa                                                                                                                                   >  1:1297705/1‑148 (MQ=255)
              ccACGCCTTACACCGTACATTCTTGTTGTGCGTTATGTAATTATTCTTCGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAAt                                                                                                                                 <  1:1402438/148‑1 (MQ=255)
              ccACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAAt                                                                                                                                 <  2:1009642/148‑1 (MQ=255)
                aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATcc                                                                                                                               <  1:687938/148‑1 (MQ=255)
                 cGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCt                                                                                                                              <  2:371478/148‑1 (MQ=255)
                    cTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTgg                                                                                                                           <  2:633115/148‑1 (MQ=255)
                    cTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTgg                                                                                                                           >  1:3881/1‑148 (MQ=255)
                      tACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGtt                                                                                                                         >  1:697979/1‑148 (MQ=255)
                      tACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCGGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGtt                                                                                                                         <  2:696888/148‑1 (MQ=255)
                       acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTAGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttt                                                                                                                        >  2:2663175/1‑148 (MQ=255)
                        cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCAGTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGtttt                                                                                                                       <  1:3108328/148‑1 (MQ=255)
                           cGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                    >  2:347708/1‑148 (MQ=255)
                             tACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaaa                                                                                                                  <  1:129233/148‑1 (MQ=255)
                             tACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaaa                                                                                                                  >  1:226064/1‑148 (MQ=255)
                             tACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaaa                                                                                                                  <  1:1115302/148‑1 (MQ=255)
                              aCATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaaaa                                                                                                                 >  2:2215064/1‑148 (MQ=255)
                                aTTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAAtt                                                                                                               <  1:364114/148‑1 (MQ=255)
                                 ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtg                                                                                                              >  1:353367/1‑148 (MQ=255)
                                 ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtg                                                                                                              >  1:272959/1‑148 (MQ=255)
                                 ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtg                                                                                                              >  1:1012338/1‑148 (MQ=255)
                                  tCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTGTTTAAAAATtgt                                                                                                             >  2:427398/1‑148 (MQ=255)
                                       ttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTTCATa                                                                                                        <  1:2650471/148‑1 (MQ=255)
                                       ttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATa                                                                                                        <  1:2673765/148‑1 (MQ=255)
                                       ttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATa                                                                                                        >  2:286517/1‑148 (MQ=255)
                                        tgtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAg                                                                                                       <  2:3370798/148‑1 (MQ=255)
                                         gtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGt                                                                                                      <  2:1057302/148‑1 (MQ=255)
                                              gTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGaaaa                                                                                                 <  1:2674101/148‑1 (MQ=255)
                                                aTGTAATTATTCTACGCTGTCGGTATCGGTGTCTTCGTTAATCAGCAGGTTAGCCCGCCCAAAAAAGCCAGTGAGAAAATTATTGATTTTACATGGGCTTATTATATTTCTATTCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATc                                                                                               >  1:1195704/1‑148 (MQ=255)
                                                  gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                             >  1:28927/1‑148 (MQ=255)
                                                    aaTTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                           <  1:253483/148‑1 (MQ=255)
                                                    aaTTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                           <  2:603195/148‑1 (MQ=255)
                                                    aaTTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                           <  1:511494/148‑1 (MQ=255)
                                                       tattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTaa                                                                                        <  1:420213/148‑1 (MQ=255)
                                                            tACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgtg                                                                                   <  2:65250/148‑1 (MQ=255)
                                                              cGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGtgtg                                                                                 >  1:1659428/1‑148 (MQ=255)
                                                                  gTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                             <  2:1297753/148‑1 (MQ=255)
                                                                  gTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                             <  2:1297695/148‑1 (MQ=255)
                                                                   tCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTa                                                                            <  1:1211279/148‑1 (MQ=255)
                                                                   tCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTa                                                                            <  1:2284022/148‑1 (MQ=255)
                                                                       tATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                        >  2:3136238/1‑148 (MQ=255)
                                                                        aTCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaa                                                                       >  2:1792303/1‑148 (MQ=255)
                                                                          cGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaaa                                                                     <  2:105025/148‑1 (MQ=255)
                                                                            gtgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATc                                                                   <  1:139316/148‑1 (MQ=255)
                                                                               tATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcaca                                                                <  1:1643623/148‑1 (MQ=255)
                                                                                aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAg                                                               <  1:334589/148‑1 (MQ=255)
                                                                                aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAg                                                               >  2:895752/1‑148 (MQ=255)
                                                                                aTTCGTTAATCAGCAGGTGAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAg                                                               <  1:741337/148‑1 (MQ=255)
                                                                                  tCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCg                                                             <  2:3012382/148‑1 (MQ=255)
                                                                                   cGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGt                                                            <  1:2347226/148‑1 (MQ=255)
                                                                                        aTCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCtt                                                       <  2:1015926/148‑1 (MQ=255)
                                                                                            gcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaa                                                   >  2:3461/1‑148 (MQ=255)
                                                                                             cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaaa                                                  >  1:2912163/1‑148 (MQ=255)
                                                                                             cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaaa                                                  >  2:3490982/1‑148 (MQ=255)
                                                                                             cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaaa                                                  <  2:1885970/148‑1 (MQ=255)
                                                                                              agGTTAGCCAGCCAAAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAt                                                 <  2:1933790/148‑1 (MQ=255)
                                                                                               tgTTAGCCAGCCACAAAAAGCCATTGAGAAATTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTTGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAtt                                                <  2:1387282/147‑1 (MQ=255)
                                                                                               ggTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAtt                                                >  1:2065755/1‑148 (MQ=255)
                                                                                                     ccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                          >  2:2265839/1‑148 (MQ=255)
                                                                                                     ccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                          >  1:3231491/1‑148 (MQ=255)
                                                                                                     ccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                          >  1:1546255/1‑148 (MQ=255)
                                                                                                           acaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                    <  2:2462176/148‑1 (MQ=255)
                                                                                                            caAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACt                                   >  2:258863/1‑148 (MQ=255)
                                                                                                             aaaaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACtt                                  >  2:203543/1‑148 (MQ=255)
                                                                                                              aaaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTg                                 >  1:2657325/1‑148 (MQ=255)
                                                                                                              aaaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTg                                 <  2:1451098/148‑1 (MQ=255)
                                                                                                               aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGc                                >  1:2717588/1‑148 (MQ=255)
                                                                                                               aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGc                                >  1:135173/1‑148 (MQ=255)
                                                                                                                aaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCa                               >  1:1005967/1‑148 (MQ=255)
                                                                                                                  gCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGt                             >  2:1692940/1‑148 (MQ=255)
                                                                                                                   ccATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                            >  2:2141585/1‑148 (MQ=255)
                                                                                                                    cATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTa                           <  1:3060877/148‑1 (MQ=255)
                                                                                                                      ttGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAAt                         >  1:3130810/1‑148 (MQ=255)
                                                                                                                         agaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGa                      <  2:577959/148‑1 (MQ=255)
                                                                                                                         agaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGa                      >  1:534611/1‑148 (MQ=255)
                                                                                                                           aaaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAAt                    <  1:3047864/148‑1 (MQ=255)
                                                                                                                           aaaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAAt                    >  2:2122092/1‑148 (MQ=255)
                                                                                                                            aaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                   <  2:1612694/148‑1 (MQ=255)
                                                                                                                            aaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                   <  1:2272748/148‑1 (MQ=255)
                                                                                                                            aaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                   <  1:2129624/148‑1 (MQ=255)
                                                                                                                              attattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCt                 <  1:286517/148‑1 (MQ=255)
                                                                                                                               ttattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCtt                >  2:2905460/1‑148 (MQ=255)
                                                                                                                                tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTa               <  2:449312/148‑1 (MQ=255)
                                                                                                                                   tGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg            >  2:1088319/1‑148 (MQ=255)
                                                                                                                                   tGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg            >  2:1620965/1‑148 (MQ=255)
                                                                                                                                   tGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg            >  1:918238/1‑148 (MQ=255)
                                                                                                                                    gATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTgg           <  2:3302043/148‑1 (MQ=255)
                                                                                                                                     aTTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGa          <  2:226064/148‑1 (MQ=255)
                                                                                                                                      ttttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGAt         <  2:470257/148‑1 (MQ=255)
                                                                                                                                         tACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGaa      <  1:1029505/148‑1 (MQ=255)
                                                                                                                                            aTGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAggg   >  2:1585266/1‑148 (MQ=255)
                                                                                                                                            aTGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAggg   <  1:952253/148‑1 (MQ=255)
                                                                                                                                             tGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCATTATTGTATGAATGga  >  2:3262207/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGA  >  CP009273/4417406‑4417694

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: