Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 6.2% 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
*NC_000913915,2150TG 6.2% 723.9 / 12.3 252intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (186/46);  new base G (3/17);  total (189/63)
Fisher's exact test for biased strand distribution p-value = 5.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915066‑915363
                                                                                                                                                     |                                                                                                                                                     
ggggCAGCCGACATGGATGTGGGCTGAGGGCGCGTTTTTCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCTCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       <  2:7191257/150‑1 (MQ=12)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGCGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       <  1:3140309/150‑1 (MQ=17)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      >  1:4262082/1‑150 (MQ=32)
  ggCAGCCCCCATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAGGGGATAAGTCGATGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTga                                                                                                                                                     <  2:2701451/150‑2 (MQ=12)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                     <  1:478734/150‑1 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                    >  1:4923710/1‑150 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGgccgc                                                                                                                                                 >  2:3165206/1‑150 (MQ=17)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                 >  1:7168458/1‑150 (MQ=255)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  1:4808351/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  1:4039653/1‑150 (MQ=32)
         caccTGGCTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTACGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGcccc                                                                                                                                              <  2:449256/150‑1 (MQ=255)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              <  2:3898464/150‑1 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              <  2:299641/150‑1 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              <  1:4392644/150‑1 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              <  2:2221114/150‑1 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTGACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              <  2:104327/150‑1 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                        >  1:3204200/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                        >  2:4089520/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                        >  2:2392399/1‑150 (MQ=32)
                  gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCt                                                                                                                                    >  2:3811184/1‑150 (MQ=32)
                  gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCt                                                                                                                                    >  2:397444/1‑150 (MQ=32)
                    gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCGcc                                                                                                                                  >  1:6248951/1‑150 (MQ=17)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTg                                                                                                                               >  1:3605162/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAG‑CGCGCTCCATg                                                                                                                               >  1:1226817/1‑150 (MQ=12)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTgc                                                                                                                              >  1:4107665/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTgc                                                                                                                              >  2:4369230/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTgc                                                                                                                              >  1:926957/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTGCGCCAg                                                                                                                         >  1:2446947/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGt                                                                                                                        >  1:6357085/1‑150 (MQ=32)
                               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTg                                                                                                                       >  1:2600197/1‑150 (MQ=32)
                               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGG‑CGCGCTCCATGCGCCAGTg                                                                                                                       >  1:2798579/1‑150 (MQ=17)
                                 gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGc                                                                                                                     >  2:5288234/1‑150 (MQ=33)
                                  ttttACAGGGTTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGTATAAGTCGTGGGCCCCGCGCAGCGGCTTTTTTGTTCGCCTGCGCCGGGGGGTGCTGGGGGCGGGGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCa                                                                                                                    <  2:4717773/150‑1 (MQ=11)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCCGTGGcc                                                                                                                    >  1:5490638/1‑149 (MQ=18)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCa                                                                                                                    >  2:5746932/1‑150 (MQ=34)
                                     tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAata                                                                                                                 >  1:5731114/1‑150 (MQ=37)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  1:3811481/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  1:2206113/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCACAGCCCGGGGGTGCCGGGGGCGGCGGCGATTTGGCGCCCCCAG‑CGCGCGCCGAGCGCCCGTGGCAataatt                                                                                                              >  1:2220244/1‑147 (MQ=11)
                                            aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                          >  1:5106592/1‑150 (MQ=255)
                                              gTTACCTCGCGCCCGACCGGGTAGCCGTAAGGTATATGTCGAGGGCACCGCGCAGGGGCGATTTTGTTCGCCTGTGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                        <  2:5452517/150‑1 (MQ=14)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTc                                                                                                    >  2:1528399/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCCCCAGTGGCAATATGTTGCTTAGCTc                                                                                                    >  2:869434/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAG‑CGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  >  1:4883063/1‑150 (MQ=255)
                                                    ccgcgCCCGCCCCGGTAGCCTTAGGGTAAACGGGGTGGCCACGGCGCCGCGGCTTTTTTGTGCGCCAGCGCCGGGGGGGGCGGGGGGGGGGGGCTTTGGGCCCCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  <  2:976000/149‑1 (MQ=9)
                                                      gcgcCCGACCGGGTTGCCTTTCGGTATACGTCGTGGGCACCGCGCCGCGGGTATTTTGTTCGCCAGCGCCGGGGGGTGCGGGGGGGGGGGGCGATGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGa                                                                                                <  2:5769907/150‑1 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCGG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCGTAGCTCAAGAAAGGAg                                                                                          >  2:2507882/1‑150 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAG‑CGCGCTCCGTGCGCCCGTGGCAATATGTTGCGTAGCTCATGTAAggcg                                                                                          >  1:5085546/1‑148 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  1:4445157/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGG‑CGCGCTCCGTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  1:5132736/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAG‑CGCGCACCATGCGCCAGGGGCAAAATGTTTCTTAGCTCATGCAAGGAg                                                                                          >  1:1327503/1‑150 (MQ=14)
                                                             tcccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCCGGGGCAATATGTTGGTTAGGTCATGTAAGGGgc                                                                                         >  1:3663642/2‑150 (MQ=14)
                                                                 ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGcaa                                                                                     >  1:1681181/1‑150 (MQ=255)
                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                    >  2:5224922/1‑150 (MQ=255)
                                                                     gCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAg                                                                                 >  2:4202773/1‑150 (MQ=255)
                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                >  2:2886445/1‑150 (MQ=255)
                                                                       ggtctgggATAGGTCGCGGGCAGCGGGCCGGGGCGTTTTTGGGCGCCAGCGCCGGGGGGGGCGGGGGGGGGGGGCGTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                               <  2:585859/145‑1 (MQ=18)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGCTTGGCCGCCCCCAG‑CGCGCTCCATACCCCAGGGGCAAAAAGTTTATTATCTCATGAAAGGCGCGCAACAAgaa                                                                               >  1:5215877/1‑149 (MQ=17)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCGCCATGCGCCAGTGGCAATAAGTTTCTTAGCTAATGAAAAGAGCGCAACAAgatg                                                                              >  1:5802827/1‑150 (MQ=17)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  2:6007528/1‑150 (MQ=255)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  2:4240559/1‑150 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTTGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                        >  2:3011256/1‑150 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGTGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                        <  2:2297498/150‑1 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGCAAGGAGCGCAACAAGATGATGAAt                                                                        >  2:3245727/1‑150 (MQ=255)
                                                                                   gtggggggCCCCGCGCCGGGGCTTTTTGGTTCGCCAGCGCCGGGGGTTGCGGGGGGGGGGGGGTTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                   <  2:5067303/145‑1 (MQ=18)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGGGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACag                                                                  >  2:340015/1‑149 (MQ=255)
                                                                                         ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaaa                                                             >  2:2021014/1‑149 (MQ=255)
                                                                                           cACCGCGCAGCGGCTTTTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGGGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                           <  2:3060871/150‑1 (MQ=21)
                                                                                           cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                           >  2:749802/1‑150 (MQ=255)
                                                                                              cgggccgcggcgTTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGGGTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        <  2:7222047/144‑1 (MQ=21)
                                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        >  2:2188360/1‑150 (MQ=255)
                                                                                               gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                       >  2:6681029/1‑150 (MQ=255)
                                                                                                     cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                 >  2:5643711/1‑150 (MQ=255)
                                                                                                        cgtttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                              <  2:5407641/147‑1 (MQ=255)
                                                                                                        cgATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             <  2:1809294/150‑1 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             >  2:6937806/1‑150 (MQ=255)
                                                                                                          tttttGTCCGCCCGTGCGGGGGGGGGGGGGGGGGGGGGGCGTTGGCCCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                            <  1:1838230/149‑1 (MQ=255)
                                                                                                          aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                            >  2:4092173/1‑150 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                  >  2:2938020/1‑150 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCGCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                  >  2:3332392/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                 >  2:6487255/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGCGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                 <  2:148229/150‑1 (MQ=255)
                                                                                                                            cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                          >  2:5185433/1‑150 (MQ=255)
                                                                                                                             gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATAACGCCAGGcc                         >  2:1340180/1‑150 (MQ=255)
                                                                                                                                gggggggggggcgggggcG‑TGGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                      <  1:6718614/144‑1 (MQ=255)
                                                                                                                                  tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                    >  1:2903900/1‑150 (MQ=255)
                                                                                                                                  tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                    >  1:2903806/1‑150 (MQ=255)
                                                                                                                                      ggggggggCGGGTCTTGCCCCCCCCCCG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  1:6835904/150‑1 (MQ=255)
                                                                                                                                      gggggCGGCGGGGTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  1:535255/150‑1 (MQ=255)
                                                                                                                                        gggCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa              >  2:5360222/1‑150 (MQ=255)
                                                                                                                                         ggcggcggcGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat             >  1:2139526/1‑150 (MQ=255)
                                                                                                                                          ggggcggccAATCGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            >  2:414432/3‑150 (MQ=255)
                                                                                                                                          gcggcggctTGTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            <  2:6370744/150‑1 (MQ=255)
                                                                                                                                             gcggcgCTTGGCCGCCCCCTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat         <  1:6598237/150‑1 (MQ=255)
                                                                                                                                             gcggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat         >  2:6131916/1‑150 (MQ=255)
                                                                                                                                               ggcgtgtGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa       <  2:4923710/144‑1 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa       >  2:1190456/1‑150 (MQ=255)
                                                                                                                                                  taTTGGCCGCCCCCGC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc    <  1:5948742/149‑1 (MQ=255)
                                                                                                                                                  gATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc    >  2:5449888/1‑150 (MQ=255)
                                                                                                                                                    ttGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  2:6777117/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915066‑915363

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