Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 910,671 T→G 12.3% intergenic (‑358/+137) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273910,6710TG12.3% 351.9 / 11.9 159intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (104/32);  new base G (2/21);  total (106/53)
Fisher's exact test for biased strand distribution p-value = 5.84e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  CP009273/910522‑910820
                                                                                                                                                     |                                                                                                                                                         
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                           <  2:1686888/150‑1 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                           >  1:4509207/1‑150 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                           >  1:4509163/1‑150 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                           >  1:4049125/1‑150 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                          >  2:489910/1‑150 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                         >  1:1042246/1‑150 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                         >  1:2530302/1‑150 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                       >  2:3709587/1‑150 (MQ=32)
     aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg                                                                                                                                                     <  2:1810237/150‑1 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                   <  1:2988088/150‑1 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                                  <  2:3718062/150‑1 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTGTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                                 <  1:2606419/150‑1 (MQ=17)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                            <  2:930734/150‑1 (MQ=32)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                            <  1:302359/150‑1 (MQ=32)
             tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGGCGCCCCCAgc                                                                                                                                          >  1:1194355/1‑150 (MQ=12)
                aTGTGGGCTGAGGGCGCGTTTTACAGGCATGTTACCTCGCGCCCGACCCGGGAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcg                                                                                                                                       <  2:3774509/150‑1 (MQ=17)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGGCGCCCCCTGCGCGGTccat                                                                                                                                 >  1:551507/1‑148 (MQ=12)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                               >  2:4287782/1‑150 (MQ=32)
                           gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGcc                                                                                                                            >  1:165797/1‑150 (MQ=32)
                           gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGcc                                                                                                                            >  1:2391552/1‑150 (MQ=17)
                              cgcgTTTAACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                         >  2:3857486/1‑150 (MQ=17)
                               gggTTTTACTGGGTTTTTCCCCCGCGCCCCACCGGGTCGCGGTATGGGTTATGTGGGGGCCCCCGCGCTGCGCCTTTTTGGTCCGCCCGCCCCGGGGGGGGCGGGGGGCGGCGGCTTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTg                                                                                                                        <  2:1868872/148‑1 (MQ=2)
                                 gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                      >  2:3882103/1‑150 (MQ=33)
                                    ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAGCGCGCTCCGTGCGCCCGTGGCAat                                                                                                                   >  1:3248690/1‑150 (MQ=255)
                                       cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                                >  1:2942178/1‑150 (MQ=38)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCATGCGCCAGTGGCAATATGt                                                                                                               >  1:1238166/1‑150 (MQ=255)
                                           cttggtACCTCGCGCCCGACCCTTTAGCTTTTTGTTATATGCCGTGGTCACCGCGCAGGGGCTTTTTTGTTCGAAAGCGCCCGGGGGTCATGGGGGGGGTGGGTTTGGGCCGCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGc                                                                                                            <  2:1484244/145‑1 (MQ=2)
                                            aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTTGCAATATGTTTCt                                                                                                           >  1:3084409/1‑150 (MQ=255)
                                             ttttACCTCGCCCCCGCCGGGTTCGCCTTAGGGTATTCGTGGGGGGCCCCGCGCGGGGGTTTTTTTGTTCGCCTGCGCGGGGGGGGGCGGGGGGGGGCGGCGTTGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                          <  2:2933542/148‑1 (MQ=255)
                                             tGTTCCCTCGCGCCCGCCCGGTTCCCCTTATGGTATACGTCGTGGCCCCCGCGCCGGGGCTTTTTTGTCCGCCTGCCCCGGGGGGTGCGGGGGGCGGCGGCTTTGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                          <  2:725611/150‑1 (MQ=255)
                                             tGTTCCCTCGCGCCCCACCGGGTCGCGTTCTGGTATTCGTGGGGGCCACCGCGCCGCGGCTTTTTTGTTCGCCCGCGCCGGGGGGGGCGGGGGGGGGCGGCGTTGGGCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                          <  2:1880113/150‑1 (MQ=255)
                                                tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAATTCGAGGGCACCGCGAAGCGGCGATTTTGTTAGCCAGAGCCCCGGGGTGCAGGGGGCGGCGGCGATGTGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAAAATTTTTCTTAGc                                                                                                       >  2:2336548/1‑150 (MQ=11)
                                                tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGCTTTAGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCCCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                                       >  2:2336549/1‑150 (MQ=255)
                                                      gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCACCTTGCGCCAGTGGCAAAATGTTTCTTATCTCATGa                                                                                                 >  1:4363702/1‑150 (MQ=14)
                                                         cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAg                                                                                              >  1:3686728/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  1:1355743/1‑150 (MQ=255)
                                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                            >  2:266782/1‑150 (MQ=255)
                                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                            >  2:785366/1‑150 (MQ=255)
                                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                            >  2:801451/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  1:1107740/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGAGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2237142/1‑150 (MQ=255)
                                                              cccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         >  2:2627941/1‑150 (MQ=255)
                                                                  gTAGCCTTATGGGATAAGTCGTGGGCCCCGCGCGGCGGCTATTTTGTTCGCCTGAGCCGGGGGTTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                     <  2:3045904/150‑1 (MQ=17)
                                                                    aGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaacaa                                                                                   >  2:2940165/1‑150 (MQ=255)
                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCGGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCATAGCTCATGAAAGGCGCGCAACAAga                                                                                 >  1:176018/1‑150 (MQ=17)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                >  2:1448902/1‑150 (MQ=255)
                                                                                tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCa                                                                       >  2:2863438/1‑150 (MQ=255)
                                                                                tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCa                                                                       >  2:981714/1‑150 (MQ=255)
                                                                                 aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      >  2:1591333/1‑150 (MQ=255)
                                                                                  cgTCGTGGCCACCGCGCTGCGGTTTTTTTGTCCGCCTGCGCGGGGGGGTGCTGGGGGCGGGGGGTTTGGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                                                                     <  2:1046653/149‑1 (MQ=18)
                                                                                  aGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCCCCCCCAGCGCGCGCCTTGCGCCAGTGGCAATAAGTTGGTTAGCTCAGGAAAGGGGCGCAACAAAAAGTTGAATCAAc                                                                     >  1:388455/1‑150 (MQ=17)
                                                                                     cGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCAATAAGTTGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGAATCAACata                                                                  >  2:1563597/1‑150 (MQ=255)
                                                                                      gtggccACCGCGCTGCGGTTTTTTTGTCCCCTTGCGCGGGGGGGTTTTGGGGGCGGGGGCTTTTGCCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                 <  2:1798972/145‑1 (MQ=255)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGAGTGGGCGCCCCCAGCGCCCTCCTTGCGCCAGTGGCAATATGTTTGATATGTCAAGTAAGGAGCGCAACAACATGAAGAAAAAACataa                                                                 >  1:921174/1‑149 (MQ=14)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCGCATGTAAGGCGCGCAACAAGATGATGAATAAACatat                                                                 >  1:920794/1‑150 (MQ=18)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACaaat                                                                 >  1:3330088/1‑150 (MQ=255)
                                                                                       aGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGTAAGGAGCGCAACAAGTTGTTGAATCAACAtata                                                                >  2:3752588/1‑150 (MQ=255)
                                                                                           cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                            >  2:2815857/1‑150 (MQ=255)
                                                                                           cACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAGCGAGCAACTTGCGCCAGTTGCAATATGTTGCTTAGCTCAAGAAACGAGCGCAACAAGATGAATAAACAACATATaacaa                                                            >  2:2815823/1‑150 (MQ=18)
                                                                                                   cgcggcgTTTGTGTGCGCCTGGGCGGGGGGGGGCGGGGGGGGGGGGCGTTTGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                                    <  2:2584392/149‑1 (MQ=25)
                                                                                                       gtttttttGTGCGCCTGCGCGGGGGGGTGCGGGGGGGGGCGGCTTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                <  2:3281692/146‑1 (MQ=25)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              >  1:2629270/1‑150 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              >  2:3181118/1‑150 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              >  1:4222046/1‑150 (MQ=255)
                                                                                                          tttttGTCCCCCAGCCCCGGGGGGTGCTGGGGGGGGGGGCTATGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                             <  1:1902193/149‑1 (MQ=255)
                                                                                                           ttttGTTCGCCAGCGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                            <  2:4446638/150‑1 (MQ=255)
                                                                                                              tGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAATGGCCTg                                         >  2:440638/1‑150 (MQ=255)
                                                                                                              tGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                                         >  2:4017367/1‑150 (MQ=255)
                                                                                                                ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTAATGAAAGGAGCGCAACAAGATGATGTATAAACATATAACAACATCTTAAAAAAATGCCTGAc                                       >  1:1971902/1‑150 (MQ=25)
                                                                                                                 tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACa                                      >  2:731542/1‑150 (MQ=255)
                                                                                                                 tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACa                                      >  2:132480/1‑150 (MQ=255)
                                                                                                                 tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTAATGAAAGGAGCGCAACAAGATGATGAATCAACATATAAAAACATCTTAAAAAAAGGCCTGaaa                                      >  2:132829/1‑148 (MQ=255)
                                                                                                                 gcGCCCGCGCCGGGGGGTGCTGGGGGGGGCGGCTATGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACa                                      <  2:354786/149‑1 (MQ=255)
                                                                                                                 ccGCCAGCGCCCGGGGGGGCTGGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACa                                      <  2:3456374/149‑1 (MQ=255)
                                                                                                                  cGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTATCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                                     >  2:3449931/1‑150 (MQ=255)
                                                                                                                       gagCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAGAGGCCTGACATTACGc                                >  2:4563293/1‑150 (MQ=255)
                                                                                                                        cgCCGGGGGGTGCTGGGGGGGGCGGTTTTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGTGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                               <  2:1741244/149‑1 (MQ=255)
                                                                                                                         gggggggggggcgggggggggcggGGTTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:4021620/146‑1 (MQ=255)
                                                                                                                          cggggggggGCTGGGGGGGGGGGCGTTGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                             <  1:4180968/147‑1 (MQ=38)
                                                                                                                            cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGACATTACGCCAGGc                           >  1:4087950/1‑150 (MQ=255)
                                                                                                                            cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                           >  2:33005/1‑150 (MQ=255)
                                                                                                                             ggggggGCTGGGGGGGGGGGCTTTGGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          <  2:2183641/150‑1 (MQ=38)
                                                                                                                             gggggTGCTGGGGGGGGGGGATTTTGCCCGCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          <  2:1135596/150‑1 (MQ=255)
                                                                                                                              ggggTGGGGGGGGGGGCGGCTATGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCt                         <  2:481998/150‑1 (MQ=255)
                                                                                                                               gggTGCTGGGGGCGGGGGGTTTGGGCCGCCCCTTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                        <  2:2219286/150‑1 (MQ=38)
                                                                                                                                gggggggggggCGGGGGGTTTGGCCCCCCCCGCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                       <  1:613467/144‑1 (MQ=255)
                                                                                                                                ggTGCTGGGGGGGGGGGGTTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCGACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                       <  1:613533/150‑1 (MQ=255)
                                                                                                                                 ggggcggggggggggggggTTGCCCCCCCCACCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                       <  1:1797818/145‑1 (MQ=255)
                                                                                                                                  tggggggggCGGCGGTGTTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     <  1:2420291/146‑1 (MQ=255)
                                                                                                                                  gtgcggggggggggggGTTTGGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                      <  2:210101/146‑1 (MQ=255)
                                                                                                                                    cAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                   >  1:1593813/1‑150 (MQ=255)
                                                                                                                                      ggggggggggGCGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                 <  1:358068/150‑1 (MQ=39)
                                                                                                                                      gggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                 >  2:3357200/1‑150 (MQ=255)
                                                                                                                                            ggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGACAGGCCTTCTGCGttaattaa           >  2:2394631/1‑150 (MQ=255)
                                                                                                                                              cggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc         >  1:3588013/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  2:358344/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  1:669404/1‑150 (MQ=255)
                                                                                                                                                   aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt    >  2:569661/1‑150 (MQ=255)
                                                                                                                                                    ttGGCCGCCCCCTGCGCGCTCCTTGCGCCANTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGACATTACGCCAGGCCTTCAGCGTTAATTAATCACGCtt   >  1:4278901/1‑150 (MQ=255)
                                                                                                                                                     tGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  >  2:3963841/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  CP009273/910522‑910820

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