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.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  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)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  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)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                        >  1:3204200/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                        >  2:4089520/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                        >  2:2392399/1‑150 (MQ=32)
                  gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCt                                                                                                                                    >  2:3811184/1‑150 (MQ=32)
                  gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCt                                                                                                                                    >  2:397444/1‑150 (MQ=32)
                    gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGcc                                                                                                                                  >  1:6248951/1‑150 (MQ=17)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                               >  1:3605162/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAGCGCGCTCCATg                                                                                                                               >  1:1226817/1‑150 (MQ=12)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                              >  1:4107665/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                              >  2:4369230/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                              >  1:926957/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAg                                                                                                                         >  1:2446947/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                        >  1:6357085/1‑150 (MQ=32)
                               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg                                                                                                                       >  1:2600197/1‑150 (MQ=32)
                               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCATGCGCCAGTg                                                                                                                       >  1:2798579/1‑150 (MQ=17)
                                 gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                     >  2:5288234/1‑150 (MQ=33)
                                  ttttACAGGGTTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGTATAAGTCGTGGGCCCCGCGCAGCGGCTTTTTTGTTCGCCTGCGCCGGGGGGTGCTGGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                    <  2:4717773/150‑1 (MQ=11)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGcc                                                                                                                    >  1:5490638/1‑149 (MQ=18)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                    >  2:5746932/1‑150 (MQ=34)
                                     tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAata                                                                                                                 >  1:5731114/1‑150 (MQ=37)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  1:3811481/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  1:2206113/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCACAGCCCGGGGGTGCCGGGGGCGGCGGCGATTTGGCGCCCCCAGCGCGCGCCGAGCGCCCGTGGCAataatt                                                                                                              >  1:2220244/1‑147 (MQ=11)
                                            aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                          >  1:5106592/1‑150 (MQ=255)
                                              gTTACCTCGCGCCCGACCGGGTAGCCGTAAGGTATATGTCGAGGGCACCGCGCAGGGGCGATTTTGTTCGCCTGTGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                        <  2:5452517/150‑1 (MQ=14)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTc                                                                                                    >  2:1528399/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCCCCAGTGGCAATATGTTGCTTAGCTc                                                                                                    >  2:869434/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  >  1:4883063/1‑150 (MQ=255)
                                                    ccgcgCCCGCCCCGGTAGCCTTAGGGTAAACGGGGTGGCCACGGCGCCGCGGCTTTTTTGTGCGCCAGCGCCGGGGGGGGCGGGGGGGGGGGGCTTTGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  <  2:976000/149‑1 (MQ=9)
                                                      gcgcCCGACCGGGTTGCCTTTCGGTATACGTCGTGGGCACCGCGCCGCGGGTATTTTGTTCGCCAGCGCCGGGGGGTGCGGGGGGGGGGGGCGATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGa                                                                                                <  2:5769907/150‑1 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCGGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCGTAGCTCAAGAAAGGAg                                                                                          >  2:2507882/1‑150 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAGCGCGCTCCGTGCGCCCGTGGCAATATGTTGCGTAGCTCATGTAAggcg                                                                                          >  1:5085546/1‑148 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  1:4445157/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  1:5132736/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAGCGCGCACCATGCGCCAGGGGCAAAATGTTTCTTAGCTCATGCAAGGAg                                                                                          >  1:1327503/1‑150 (MQ=14)
                                                             tcccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCCGGGGCAATATGTTGGTTAGGTCATGTAAGGGgc                                                                                         >  1:3663642/2‑150 (MQ=14)
                                                                 ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGcaa                                                                                     >  1:1681181/1‑150 (MQ=255)
                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                    >  2:5224922/1‑150 (MQ=255)
                                                                     gCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAg                                                                                 >  2:4202773/1‑150 (MQ=255)
                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                >  2:2886445/1‑150 (MQ=255)
                                                                       ggtctgggATAGGTCGCGGGCAGCGGGCCGGGGCGTTTTTGGGCGCCAGCGCCGGGGGGGGCGGGGGGGGGGGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                               <  2:585859/145‑1 (MQ=18)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGCTTGGCCGCCCCCAGCGCGCTCCATACCCCAGGGGCAAAAAGTTTATTATCTCATGAAAGGCGCGCAACAAgaa                                                                               >  1:5215877/1‑149 (MQ=17)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCAGTGGCAATAAGTTTCTTAGCTAATGAAAAGAGCGCAACAAgatg                                                                              >  1:5802827/1‑150 (MQ=17)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  2:6007528/1‑150 (MQ=255)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  2:4240559/1‑150 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTTGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                        >  2:3011256/1‑150 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGTGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                        <  2:2297498/150‑1 (MQ=255)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGCAAGGAGCGCAACAAGATGATGAAt                                                                        >  2:3245727/1‑150 (MQ=255)
                                                                                   gtggggggCCCCGCGCCGGGGCTTTTTGGTTCGCCAGCGCCGGGGGTTGCGGGGGGGGGGGGGTTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                   <  2:5067303/145‑1 (MQ=18)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGGGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACag                                                                  >  2:340015/1‑149 (MQ=255)
                                                                                         ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaaa                                                             >  2:2021014/1‑149 (MQ=255)
                                                                                           cACCGCGCAGCGGCTTTTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                           <  2:3060871/150‑1 (MQ=21)
                                                                                           cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                           >  2:749802/1‑150 (MQ=255)
                                                                                              cgggccgcggcgTTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGGGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        <  2:7222047/144‑1 (MQ=21)
                                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        >  2:2188360/1‑150 (MQ=255)
                                                                                               gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                       >  2:6681029/1‑150 (MQ=255)
                                                                                                     cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                 >  2:5643711/1‑150 (MQ=255)
                                                                                                        cgtttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                              <  2:5407641/147‑1 (MQ=255)
                                                                                                        cgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             <  2:1809294/150‑1 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             >  2:6937806/1‑150 (MQ=255)
                                                                                                          tttttGTCCGCCCGTGCGGGGGGGGGGGGGGGGGGGGGGCGTTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                            <  1:1838230/149‑1 (MQ=255)
                                                                                                          aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                            >  2:4092173/1‑150 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                  >  2:2938020/1‑150 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCGCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                  >  2:3332392/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                 >  2:6487255/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGCGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                 <  2:148229/150‑1 (MQ=255)
                                                                                                                            cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                          >  2:5185433/1‑150 (MQ=255)
                                                                                                                             gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATAACGCCAGGcc                         >  2:1340180/1‑150 (MQ=255)
                                                                                                                                gggggggggggcgggggcGTGGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                      <  1:6718614/144‑1 (MQ=255)
                                                                                                                                  tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                    >  1:2903900/1‑150 (MQ=255)
                                                                                                                                  tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                    >  1:2903806/1‑150 (MQ=255)
                                                                                                                                      ggggggggCGGGTCTTGCCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  1:6835904/150‑1 (MQ=255)
                                                                                                                                      gggggCGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  1:535255/150‑1 (MQ=255)
                                                                                                                                        gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa              >  2:5360222/1‑150 (MQ=255)
                                                                                                                                         ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat             >  1:2139526/1‑150 (MQ=255)
                                                                                                                                          ggggcggccAATCGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            >  2:414432/3‑150 (MQ=255)
                                                                                                                                          gcggcggctTGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            <  2:6370744/150‑1 (MQ=255)
                                                                                                                                             gcggcgCTTGGCCGCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat         <  1:6598237/150‑1 (MQ=255)
                                                                                                                                             gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat         >  2:6131916/1‑150 (MQ=255)
                                                                                                                                               ggcgtgtGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa       <  2:4923710/144‑1 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa       >  2:1190456/1‑150 (MQ=255)
                                                                                                                                                  taTTGGCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc    <  1:5948742/149‑1 (MQ=255)
                                                                                                                                                  gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc    >  2:5449888/1‑150 (MQ=255)
                                                                                                                                                    ttGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  2:6777117/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  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: