Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 5.8% 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 5.8% 1390.4 / 25.0 495intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (333/125);  new base G (3/34);  total (336/159)
Fisher's exact test for biased strand distribution p-value = 3.74e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                     |                                                                                                                                                     
ggggCAGCCCACATGGCTGTGGGCTGAGGGCGCGTTGTACAGGGATGTTACCTCGCGCCCGACCCGGTTGCCGTAAGGGATAAGTCGTGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       <  1:1367957/150‑1 (MQ=255)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      <  2:3241792/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTg                                                                                                                                                      >  1:8919786/1‑150 (MQ=17)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                     >  2:5261517/1‑150 (MQ=32)
   tcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                    >  1:4825752/2‑150 (MQ=17)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc                                                                                                                                                   <  2:14519770/150‑1 (MQ=32)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                 >  2:14191014/1‑150 (MQ=17)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                <  2:1685341/150‑1 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                <  1:11379136/150‑1 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                <  2:5537217/150‑1 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  1:12944965/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  1:14423946/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                              >  2:4319029/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                             >  1:6030136/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                         >  1:12543777/1‑150 (MQ=32)
               gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc                                                                                                                                        >  1:5128506/1‑150 (MQ=32)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCTGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  1:3112698/1‑150 (MQ=17)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  1:10460179/1‑150 (MQ=32)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                 >  1:2621765/1‑150 (MQ=17)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCtt                                                                                                                                 >  2:12132460/1‑150 (MQ=17)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                                >  2:8037529/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATgc                                                                                                                               >  2:1372007/1‑150 (MQ=17)
                           gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGcc                                                                                                                            >  1:8053176/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                     <  2:13571772/123‑1 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:620941/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:11358333/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:13087613/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:10767250/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  1:7795278/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                         <  2:3428476/150‑1 (MQ=17)
                                   tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                                    >  1:13348970/1‑150 (MQ=25)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCGGCGCGCTCCATGCGCCCCTGGCActattt                                                                                                               >  1:5299950/1‑148 (MQ=11)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  1:74829/1‑150 (MQ=39)
                                         ggtttGTTCCCGCGCGCCCGACCCGTTAGCCGTATGGTTTATTTGTTGGGCACCGGGCCGCGGCTTTTTTGTTCGCAAGAGCCGGGGGGTGCCGGGGGGGGGGCGGTTTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              <  2:8647722/146‑1 (MQ=9)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                             >  2:13274472/1‑150 (MQ=40)
                                             tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCtt                                                                                                          >  2:605043/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAg                                                                                                        >  2:787124/1‑150 (MQ=255)
                                                         cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAg                                                                                              >  2:12384068/1‑150 (MQ=255)
                                                           cGACCCGGTAGCCGTATGGTTTCAGTCGTGGGCACCGGGCAGCGGCTTTTTTGTTGGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGGGTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                            <  2:6550268/150‑1 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  1:13533731/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  1:9570090/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2790448/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGCCGCTTGGCCGCCCCCCGCCCGCGCCAGGCGCCCCTGGCAAAAATTTGCATAGCGCAAGTAAGGAgc                                                                                          >  1:3443397/1‑150 (MQ=11)
                                                                 gttAGCCTTATGGTTTAAGTCGAGGGCACCGCGCCGCGGCTTTTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaa                                                                                      <  2:12893323/148‑1 (MQ=17)
                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                     >  2:4794402/1‑150 (MQ=255)
                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                     >  2:10960139/1‑150 (MQ=255)
                                                                     gCCGTATGTTATATGTCGCGGCCACCGCCCAGGGGCTATTTTTTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGCCCGCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAg                                                                                  <  2:7745794/150‑1 (MQ=255)
                                                                      ccGTATGGGATAAGTCGAGGGCACCGCGCAGCGGCGTTTTTGTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                 <  2:9818332/150‑1 (MQ=18)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCCCTCCTTGCGCCAGTGGCAATAATTTTCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                >  1:10035661/1‑150 (MQ=17)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCCCGCGCCATGCGCCAGTGGCAAAAATTTTCTTAGCTCAATAAAGGATCCcaacaaaaa                                                                                >  1:15419195/1‑147 (MQ=14)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                               >  2:7305409/1‑150 (MQ=255)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGTGGCAAAAATTTTCTTAGCTAATGAAAGGGGCGcaacaacaag                                                                               >  1:5589703/1‑146 (MQ=14)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       >  1:11617291/1‑77 (MQ=34)
                                                                          aaGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                             >  2:836341/1‑150 (MQ=255)
                                                                          aaGGCATATGACGCGGGAACCGCGCCGCGGCGATTTTGTACGCCTGAGCCGGGGGGTGCCGGGGGCGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                             <  2:14049601/150‑1 (MQ=17)
                                                                            gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCGGCGCGCGCCATGCGCCAGTGGCAAAATGTTTCTTAGCGCAAGAAAGGAGCGCAACAACAagatga                                                                           >  1:7117239/1‑150 (MQ=17)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAATATGATGaa                                                                          >  2:5089502/1‑150 (MQ=255)
                                                                               aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGAGGAATc                                                                        >  2:414762/1‑150 (MQ=255)
                                                                                 aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCCCGCCATGCGCCAGTGGCAATATTTTGCTTAGCGCAAGAAAGGTGCGCCACAAAATGTTGAATaaa                                                                      >  1:3794649/1‑150 (MQ=17)
                                                                                   gTCGAGGGCACAGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCGCCCCCCGCGCGCGCCATGCGCCCCGGGCAAAATGTTGCTTAGCGCCAGCAAGGAGCGCAACAACCAGATGAATCAACa                                                                    >  1:9922541/1‑150 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTTGCCCCCCCCAGCGCGCTCCATGCGCCAGTGGCAAAATTTTGCTTAGCTCAAGAAAGGAGCGCAACAAAAATAAGAATCAACaa                                                                   >  1:10392507/1‑149 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTTGCCCCCCCCAGCGCGCACCATGCGCCACGGGCAAAAATTTTCTTATCTCATGAAAGGAGCGCAACAACAAGATGAATaaaaaa                                                                   >  1:2013842/1‑147 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:12234394/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:9161293/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:10177591/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGTGGCAATATGTTTCTTATCGCATGAAATGAGCGCAACAAGATTATGAAAAAACat                                                                   >  1:11058473/1‑150 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGCGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGAGCGCAACAAGATGAGGAATCAACat                                                                   >  2:14346417/1‑150 (MQ=255)
                                                                                            aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                           >  2:4556074/1‑150 (MQ=255)
                                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                         >  2:7003158/1‑150 (MQ=255)
                                                                                               gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                        >  2:228847/1‑150 (MQ=255)
                                                                                                    gcggctTTTTTTTTCGCCAGAGCGCGGGGGTGCTGGGGGGGGGGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                   <  2:12619674/150‑1 (MQ=25)
                                                                                                     cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                  >  2:11664541/1‑150 (MQ=255)
                                                                                                        gttttttGTTCGCCTGCGCGGGGGGGTGGTGGGGGGGGGGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  1:9142014/147‑1 (MQ=25)
                                                                                                        cttttttGTTCGCCTGAGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  2:11503154/147‑1 (MQ=255)
                                                                                                        cttttttGTGCGCCTGCGCGGGGGGGGGGTGGGGGGGGGGGGGCTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  2:9558396/147‑1 (MQ=255)
                                                                                                              tGTTCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                                         <  2:10392507/150‑1 (MQ=255)
                                                                                                               gtgccccAGAGCCGGGGGGGGCCGGGGGCGGCGGCTCTGGCCCGCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGGTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                        <  2:8841126/145‑1 (MQ=25)
                                                                                                                 tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGACa                                      >  2:4024142/1‑150 (MQ=255)
                                                                                                                   gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                    >  2:11090163/1‑150 (MQ=255)
                                                                                                                     cAGAGCCGGGGGGTGCTGGGGGGGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                  <  2:10460179/150‑1 (MQ=255)
                                                                                                                       gagCGCGGGGGTGGTGGGGGCGGCGGGGTGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                                <  2:11424308/150‑1 (MQ=37)
                                                                                                                           ccGGGGGTGCGGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                            <  2:15445871/150‑1 (MQ=255)
                                                                                                                               gggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                        <  2:9060472/150‑1 (MQ=255)
                                                                                                                                 gtggtgggggcgggggcGCTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                      <  1:12084764/145‑1 (MQ=255)
                                                                                                                                 gTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                   <  1:12416050/137‑1 (MQ=255)
                                                                                                                                  tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     >  1:11246011/1‑150 (MQ=255)
                                                                                                                                   gccggggggggCGGCGAGGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                    <  2:15358466/147‑1 (MQ=255)
                                                                                                                                    cAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                 >  2:3694645/1‑88 (MQ=255)
                                                                                                                                     aGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttt                      <  2:3870156/146‑3 (MQ=37)
                                                                                                                                     aGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttt                      <  1:14076551/146‑3 (MQ=37)
                                                                                                                                     aGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                  >  1:10764793/1‑150 (MQ=255)
                                                                                                                                       ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  2:2014661/1‑150 (MQ=255)
                                                                                                                                        gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa               >  1:5768570/1‑150 (MQ=255)
                                                                                                                                         ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat              >  1:6411797/1‑150 (MQ=255)
                                                                                                                                             gggggttttGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat          <  1:2368260/143‑1 (MQ=255)
                                                                                                                                               ggggATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        <  2:2445013/147‑1 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  1:43029/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  1:8595272/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGGTAATTAATCa        >  1:3964217/1‑150 (MQ=255)
                                                                                                                                                     tGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  >  1:14577504/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364

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: