Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,340,901 A→T 20.9% S33S (TCT→TCA gmr ← cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,340,9010AT20.9% 454.3 / 119.6 187S33S (TCT→TCAgmrcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base A (63/85);  new base T (19/20);  total (82/105)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 193 total aligned reads displayed.

TTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACA  >  CP009273/1340752‑1341050
                                                                                                                                                     |                                                                                                                                                     
tttCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTa                                                                                                                                                       <  2:873368/150‑1 (MQ=255)
tttCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTa                                                                                                                                                       <  2:5403/150‑1 (MQ=255)
     tCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGaaaa                                                                                                                                                  <  2:1600176/150‑1 (MQ=255)
      ccTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAg                                                                                                                                                 <  2:195072/150‑1 (MQ=255)
        tACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAgcg                                                                                                                                               >  1:1029655/1‑150 (MQ=255)
        tACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAgcg                                                                                                                                               >  1:3479408/1‑150 (MQ=255)
         aCCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAgcgc                                                                                                                                              <  1:885753/150‑1 (MQ=255)
         aCCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAgcgc                                                                                                                                              <  1:2225104/150‑1 (MQ=255)
           cAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCaa                                                                                                                                            >  1:2424249/1‑150 (MQ=255)
            aGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCaaa                                                                                                                                           >  1:3446935/1‑150 (MQ=255)
              aTGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAAc                                                                                                                                         <  1:2028698/150‑1 (MQ=255)
                    aGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgct                                                                                                                                   >  1:1455721/1‑150 (MQ=255)
                     gAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctg                                                                                                                                  >  1:2277552/1‑150 (MQ=255)
                     gAAAGATCGCTTTCATCGAAGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctg                                                                                                                                  >  1:3410903/1‑150 (MQ=255)
                       aaGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctgct                                                                                                                                <  2:1555107/150‑1 (MQ=255)
                         gATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCttt                                                                                                                              <  2:210386/150‑1 (MQ=255)
                           tCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTcc                                                                                                                            <  1:4072871/150‑1 (MQ=255)
                             gCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGt                                                                                                                          >  1:2587759/1‑150 (MQ=255)
                                ttCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCaa                                                                                                                       <  1:3298734/150‑1 (MQ=255)
                                  cATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAcc                                                                                                                     >  1:2195793/1‑150 (MQ=255)
                                   aTCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTGCGGTGGTAGAAAAGCGCAAAACATCGCTGCTATCCGTCAACCg                                                                                                                    >  2:1478885/1‑150 (MQ=255)
                                    tCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGc                                                                                                                   >  2:3844396/1‑150 (MQ=255)
                                     cGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGcc                                                                                                                  >  2:121769/1‑150 (MQ=255)
                                     cGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGcc                                                                                                                  >  2:2713411/1‑150 (MQ=255)
                                      gACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCa                                                                                                                 <  2:1381164/150‑1 (MQ=255)
                                       aCGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAg                                                                                                                >  1:2147115/1‑150 (MQ=255)
                                          gTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTAt                                                                                                             <  1:2085843/150‑1 (MQ=255)
                                           tCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATg                                                                                                            <  2:46723/150‑1 (MQ=255)
                                            cAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATgg                                                                                                           >  2:4261318/1‑150 (MQ=255)
                                            cAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATgg                                                                                                           <  1:3389680/150‑1 (MQ=255)
                                               aCTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAtt                                                                                                        >  2:2844616/1‑150 (MQ=255)
                                               aCTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAtt                                                                                                        <  2:135092/150‑1 (MQ=255)
                                                  catcatcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATtgtg                                                                                                     <  1:3151403/150‑1 (MQ=255)
                                                  catcatcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATtgtg                                                                                                     >  1:2490150/1‑150 (MQ=255)
                                                     catcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGa                                                                                                  <  2:3038763/150‑1 (MQ=255)
                                                      atcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt                                                                                                 <  2:3011662/150‑1 (MQ=255)
                                                      atcaGGCTGGAGGTGATGACCGTCATGTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt                                                                                                 <  2:832644/150‑1 (MQ=255)
                                                       tcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATc                                                                                                <  2:1009196/150‑1 (MQ=255)
                                                         aGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCg                                                                                              >  1:4224921/1‑150 (MQ=255)
                                                         aGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCg                                                                                              <  2:734907/150‑1 (MQ=255)
                                                            cTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCgaga                                                                                           >  2:4062805/1‑150 (MQ=255)
                                                             tGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGaa                                                                                          >  2:3574415/1‑150 (MQ=255)
                                                                  gTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTg                                                                                     <  1:841956/150‑1 (MQ=255)
                                                                         ccGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACg                                                                              >  1:3173482/1‑150 (MQ=255)
                                                                          cGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGt                                                                             <  1:1654051/150‑1 (MQ=255)
                                                                              aTTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgtg                                                                         <  2:3494871/150‑1 (MQ=255)
                                                                              aTTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgtg                                                                         <  1:2835254/150‑1 (MQ=255)
                                                                                 tCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGAc                                                                      <  2:216789/150‑1 (MQ=255)
                                                                                  cccTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACt                                                                     <  1:1437539/150‑1 (MQ=255)
                                                                                  cccTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACt                                                                     <  2:1816430/150‑1 (MQ=255)
                                                                                    cTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTcc                                                                   >  1:2074304/1‑150 (MQ=255)
                                                                                        tGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTaa                                                               <  1:2595990/150‑1 (MQ=255)
                                                                                             gcagcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGtt                                                          <  2:205925/150‑1 (MQ=255)
                                                                                                gcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTc                                                       <  2:2494766/150‑1 (MQ=255)
                                                                                                gcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTc                                                       <  2:3396635/150‑1 (MQ=255)
                                                                                                gcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTc                                                       <  2:3396892/150‑1 (MQ=255)
                                                                                                  cTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCat                                                     <  1:1205558/150‑1 (MQ=255)
                                                                                                   tGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatg                                                    >  2:621444/1‑150 (MQ=255)
                                                                                                    gTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatga                                                   <  2:1051782/150‑1 (MQ=255)
                                                                                                      tCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatgatg                                                 <  2:1800605/150‑1 (MQ=255)
                                                                                                      tCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatgatg                                                 <  2:3987964/150‑1 (MQ=255)
                                                                                                        ggCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCa                                               >  2:2367280/1‑150 (MQ=255)
                                                                                                          cAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCata                                             <  1:1275633/150‑1 (MQ=255)
                                                                                                          cAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCata                                             <  1:2755692/150‑1 (MQ=255)
                                                                                                           aGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCatat                                            >  1:2026753/1‑150 (MQ=255)
                                                                                                            gATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtc                                           <  1:2065197/150‑1 (MQ=255)
                                                                                                             aTAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcc                                          <  2:3252540/150‑1 (MQ=255)
                                                                                                              tAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcct                                         >  2:3862917/1‑150 (MQ=255)
                                                                                                              tAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcct                                         <  1:1546230/150‑1 (MQ=255)
                                                                                                              tAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcct                                         <  1:2663862/150‑1 (MQ=255)
                                                                                                               aaCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctc                                        >  2:2084611/1‑150 (MQ=255)
                                                                                                               aaCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctc                                        >  2:2423038/1‑150 (MQ=255)
                                                                                                                  tGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCgg                                     >  1:2809336/1‑150 (MQ=255)
                                                                                                                    cAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGtt                                   <  2:1958681/150‑1 (MQ=255)
                                                                                                                       tGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAAc                                <  2:2450986/150‑1 (MQ=255)
                                                                                                                        gTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACa                               <  1:2691169/150‑1 (MQ=255)
                                                                                                                        gTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACa                               <  1:3226561/150‑1 (MQ=255)
                                                                                                                        gTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACa                               <  1:520666/150‑1 (MQ=255)
                                                                                                                          aCGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGc                             <  2:3317058/150‑1 (MQ=255)
                                                                                                                           cGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                            <  1:3034313/150‑1 (MQ=255)
                                                                                                                            gATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACCTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCgg                           <  1:478331/150‑1 (MQ=255)
                                                                                                                            gATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCgg                           >  1:3690981/1‑150 (MQ=255)
                                                                                                                            gATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCgg                           <  1:664735/150‑1 (MQ=255)
                                                                                                                            gATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCgg                           >  1:1750703/1‑150 (MQ=255)
                                                                                                                             aTCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                          >  1:3605402/1‑150 (MQ=255)
                                                                                                                              tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAg                         >  2:2086898/1‑150 (MQ=255)
                                                                                                                               cAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGt                        >  1:3397145/1‑150 (MQ=255)
                                                                                                                               cAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGt                        >  1:4041441/1‑150 (MQ=255)
                                                                                                                                  gTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGTATTGTGCGATCCGAGAAAGATGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcg                     >  1:558309/1‑150 (MQ=255)
                                                                                                                                  gTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcg                     >  2:80953/1‑150 (MQ=255)
                                                                                                                                   ttCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcgc                    >  2:3547997/1‑150 (MQ=255)
                                                                                                                                       gtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGtt                <  1:2570268/150‑1 (MQ=255)
                                                                                                                                        tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta               <  1:3031344/150‑1 (MQ=255)
                                                                                                                                        tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta               <  1:18206/150‑1 (MQ=255)
                                                                                                                                           tttCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaat            <  1:4013264/150‑1 (MQ=255)
                                                                                                                                            ttCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaata           <  1:3038958/150‑1 (MQ=255)
                                                                                                                                             tCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataa          <  2:2614465/150‑1 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  2:632424/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  1:592177/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  1:793863/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACA  >  CP009273/1340752‑1341050

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: