Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 910,682 T→C 26.9% intergenic (‑369/+126) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273910,6820TC26.9% 86.4 / 23.4 74intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (37/2);  new base C (7/7);  total (59/15)
Fisher's exact test for biased strand distribution p-value = 6.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  CP009273/910538‑910830
                                                                                                                                                |                                                                                                                                                      
aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgc                                                                                                                                                  >  2:3202688/1‑149 (MQ=32)
aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgc                                                                                                                                                  <  2:3855096/149‑1 (MQ=32)
 tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgcg                                                                                                                                                 >  2:1706588/1‑149 (MQ=32)
     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTTTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGCGGGTTGGCGCCCCCCg‑‑cccccccc                                                                                                                                             >  1:1450843/1‑142 (MQ=11)
         gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGGGGCGGTTTGCCGCCCCCCG‑‑CCCGCTCCATgc                                                                                                                                         >  1:1524690/1‑149 (MQ=11)
           gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCCCCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGGGGCGGCGGGTGGGCGCCCCCGG‑‑CGCGCGCCGTgcgc                                                                                                                                       >  1:1825067/1‑149 (MQ=11)
           gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGGGGGTGGGCGCCCCCCG‑‑CGCGCGCCATgcgc                                                                                                                                       >  1:1825101/1‑149 (MQ=11)
            ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGCGTGGCCCCCCCCCG‑‑CGCGCGCCGTGCGcc                                                                                                                                      >  1:1818733/1‑149 (MQ=11)
               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGt                                                                                                                                   >  2:334123/1‑149 (MQ=32)
               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCCGt                                                                                                                                   >  2:453863/1‑149 (MQ=17)
                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                                          >  2:2783809/1‑149 (MQ=255)
                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGGTTGGGCGCCCCCCG‑‑CGCGCGCCGTGCGCCCGGGGCCAAATGt                                                                                                                         >  1:465212/1‑149 (MQ=11)
                           gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGCTTGGGCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCCATATTTTg                                                                                                                       >  1:1989186/1‑149 (MQ=11)
                             tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGCGGCGGCGCGTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                                     >  2:3451667/1‑149 (MQ=14)
                                tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                  >  2:2730466/1‑149 (MQ=255)
                                 aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCACTGGCCATATGTTCCTTAGc                                                                                                                 >  2:2114091/1‑149 (MQ=255)
                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCt                                                                                                                >  2:3168877/1‑149 (MQ=255)
                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCCGCGATTTGCCCCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTATCt                                                                                                                >  2:1109860/1‑149 (MQ=14)
                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGGGGTTTGCCCCCCCCCG‑‑CGCGCGCCCTTCGCCACGTGCCATATGTTGGTTAtctc                                                                                                               >  1:3870962/1‑149 (MQ=11)
                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGATTTGCCGCCCCCAG‑‑CGCGCTCCGTGCGCCACTGGCCATATGTTGCTTATCTCa                                                                                                              >  1:306278/1‑149 (MQ=14)
                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGGCACTATGTTGCGTAGCTCa                                                                                                              >  2:3962610/1‑149 (MQ=255)
                                      gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGGGGGGGTTTGGCGCCCCCAG‑‑CGCGCGCCGTGCGCCCGTGGCAAAATTTTTGTTTGCTCAtt                                                                                                            >  1:1219881/1‑148 (MQ=11)
                                      gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                            >  2:2526383/1‑149 (MQ=255)
                                       cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCGTGCGCCCGTGGCAAAATGTTGCGTAGCTCATGa                                                                                                           >  2:856720/1‑149 (MQ=255)
                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGGTGGCCGCCCCCCG‑‑CGCCCGCCCTGCGCCCCGGGCAAAATGTTTCTTATGTCATGGAAgg                                                                                                       >  1:1478560/1‑149 (MQ=11)
                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAg                                                                                                     >  2:3607358/1‑149 (MQ=255)
                                                 ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTG‑‑CGCGCTCCTTTCGCCAGTGGCAAAATTTTGCTTATCTAAAGAAAGGTGCGca                                                                                                 >  1:2285484/1‑149 (MQ=14)
                                                 ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGGATGGCCGCCCCCAG‑‑CGCCCACCATGCGCCCCTGGCAAAATTTTGCTTATCTCATGAAACGAGCGca                                                                                                 >  1:3200464/1‑149 (MQ=14)
                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCCCTTTGGCGCCCCCGG‑‑CGCGCGCCTTGCGCCACTGGCAAAATGTTTCTTAGCACAATAAAGGAGCGcaacaaaa                                                                                           >  2:1426214/1‑147 (MQ=17)
                                                            gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCACTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGaagatg                                                                                      >  2:3772695/1‑149 (MQ=255)
                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCGCCCCCAG‑‑CGCGCTCCCTGCGCAAGTGGCAAAAATTTGCTTATCGCATGAAAGGAGCGCAACAACAagatga                                                                                     >  1:979097/1‑149 (MQ=17)
                                                                  aGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                                >  2:3437305/1‑149 (MQ=255)
                                                                   gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATaaaa                                                                               >  2:3934954/1‑148 (MQ=255)
                                                                        gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGGCCCCTG‑‑CGCGCGCCGTGCGCCACTGGCAATATGTTGGTTAGCTAATGAAAGGAGCGCAACCAGATGATGAATCAACAtata                                                                          >  2:223423/1‑149 (MQ=21)
                                                                         ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAAATaa                                                                         >  2:941573/1‑149 (MQ=255)
                                                                          gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCGG‑‑CGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATaac                                                                        >  2:3519203/1‑149 (MQ=255)
                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGGGGCGAATTGCCGCCCCCAG‑‑CGCGCTCCATTCGCCCCTGGCCAAATTTTTCTTATCTCATGAAAGGTGCGCAACAAAATGATTAAAAAACAAATaaaaaaa                                                                    >  1:2068988/1‑147 (MQ=9)
                                                                               gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTG‑‑CCCGCTCCTTGCGCCAGTGGCAATAAGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATTATTCATACACATATAACAAAAt                                                                   >  2:2166145/1‑149 (MQ=18)
                                                                                          aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCATGCGCCAGGGGCAAAATGTTGCTTAGCTCATGTAAGGAGCGCAACAAGATGTTGTATCAACATATAAAAAACTCTTAAAAAAAg                                                        >  1:3711740/1‑149 (MQ=21)
                                                                                           ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCGCCGTGCGCCAGTGGCAATATTTTTCTTAGCTCATTAAAGGTGCGCAACAAGAATATGAATCAAAAAATAAAAACAACTTAAAAAAAgg                                                       >  1:616105/1‑149 (MQ=18)
                                                                                           ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCCATTGGGCCCCCCCAG‑‑CGCGCTCCATGCGCCACTGGCAAAATTTTTCTTATCGCAATGAATGAGCGCAACAAGATGATGTATAAAAAAATAACAACAACTTaaaaaaaat                                                       >  1:3566635/1‑147 (MQ=12)
                                                                                            tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCATGGCCGCCCCCAG‑‑CGCGCGCCTTGCGCCCGTGGCCATATTTTGCTTAGCTCATGTAAGGAGCGCAAAAACATGATGAATCAACATATAACAACAACTTAAAAAAAGGc                                                      >  2:1423392/1‑149 (MQ=25)
                                                                                             ttggccgcccgcgcccgGGGGGGCGGGGGGCGCGGGCGCGGGGCCCCCCCCCG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                                     <  2:2657110/144‑1 (MQ=255)
                                                                                                   gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTTCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACATATAACAACATCATAAAAAAAGGCCTGACAt                                               >  2:2847179/1‑149 (MQ=255)
                                                                                                   gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACCACATCTTAAAAAAAGGCCAGACAt                                               >  2:2847170/1‑149 (MQ=255)
                                                                                                    cccgcgCCGGGGGGGGCGGGGGGCGGGGGCTAGGGCCCCCCCCCCC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                              <  2:2068988/144‑1 (MQ=37)
                                                                                                        gggggggggggGCGGGGGGGGGGGGCTTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                                          <  2:3314436/144‑1 (MQ=255)
                                                                                                         gCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCACCGTTCGCCAGTGGCAAAATTTTTCTTATCTAAAGAAACGAGCGCAACAAAATTATGAAAAAACAAATAACAAAAAAATAAAAAAAAGCCAGACAATaaccc                                         >  2:19927/1‑145 (MQ=12)
                                                                                                           gcGGGGGGGGGGGGGGGGGGGGGGCTTGGGCCCCCCCCGCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                                         <  2:3977542/148‑1 (MQ=255)
                                                                                                              ggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                                    >  2:2510622/1‑149 (MQ=255)
                                                                                                                gggggggggggCGGCGGCGTGGCGCCCCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                                  <  1:1733313/143‑1 (MQ=38)
                                                                                                                   gCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                               >  1:3560242/1‑149 (MQ=255)
                                                                                                                      gggggCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                            >  1:584283/1‑149 (MQ=255)
                                                                                                                       ggggggggggCGGTGGCCCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                            <  2:798193/149‑1 (MQ=255)
                                                                                                                       ggggCGGGCGCGAGGGGCCGCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                           <  2:3312138/149‑1 (MQ=255)
                                                                                                                             gcggcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa                     >  2:3478547/1‑149 (MQ=255)
                                                                                                                             gcggcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa                     >  2:3478553/1‑149 (MQ=255)
                                                                                                                                ggttggggCCCCCCCCCC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  <  2:3486162/143‑1 (MQ=255)
                                                                                                                                gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCCGGCCTTCTGCGTTAATTAATCa                  >  2:3882852/1‑149 (MQ=255)
                                                                                                                                gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  >  2:3882253/1‑149 (MQ=255)
                                                                                                                                gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  >  1:2694649/1‑149 (MQ=255)
                                                                                                                                cctcggggCCCCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  <  2:3484712/143‑1 (MQ=255)
                                                                                                                                   aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               >  2:727920/1‑149 (MQ=255)
                                                                                                                                   aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               >  2:2789821/1‑149 (MQ=255)
                                                                                                                                   aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               >  1:2926979/1‑149 (MQ=255)
                                                                                                                                    ttGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt              >  2:1561896/1‑149 (MQ=255)
                                                                                                                                    gcgccccccccccc‑‑cGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt              <  1:856720/142‑1 (MQ=255)
                                                                                                                                     ggcccgccCCCCC‑‑CGCGCTCCTTGCGCCAGTGTCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt             <  2:381593/146‑1 (MQ=255)
                                                                                                                                     ggcccccccccGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt             <  2:1456625/143‑1 (MQ=255)
                                                                                                                                         cgccCCCCG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTAAATTAATCACGCTTTTcc         <  2:3065553/149‑1 (MQ=255)
                                                                                                                                         ccccccccc‑‑cGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         <  2:1825067/147‑1 (MQ=255)
                                                                                                                                         ccccccccc‑‑cGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         <  2:1825101/147‑1 (MQ=255)
                                                                                                                                         cccccccTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTAAATTAATCACGCTTTTcc         <  2:3065597/147‑1 (MQ=255)
                                                                                                                                         cccccccGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTAAATTAATCACGCTTTTcc         <  2:3065584/147‑1 (MQ=255)
                                                                                                                                          cccccccc‑‑cGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca        <  2:235515/148‑1 (MQ=255)
                                                                                                                                           cccccTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag       >  2:3593990/1‑149 (MQ=255)
                                                                                                                                            ccccTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagc      >  2:2236492/1‑149 (MQ=255)
                                                                                                                                              ccTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag    >  2:367014/1‑149 (MQ=255)
                                                                                                                                                tG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAggg  >  2:1373982/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
ATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  CP009273/910538‑910830

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: