Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 8.6% 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 8.6% 502.7 / 17.5 205intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (138/46);  new base G (2/19);  total (140/65)
Fisher's exact test for biased strand distribution p-value = 5.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                      |                                                                                                                                                       
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                         <  1:1287226/150‑1 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                         <  1:40935/150‑1 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       <  1:5790438/150‑1 (MQ=32)
  ctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  1:1198870/3‑150 (MQ=12)
  ctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACATAGCGCCCGACCCGGAAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCAGATCCCGGGGGTGCAGGGGGCGGCCGCGATTgg                                                                                                                                                       >  1:4562051/3‑150 (MQ=11)
  ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  2:1023937/3‑150 (MQ=12)
  ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  2:1181293/3‑150 (MQ=12)
   gCAGCCCCCATGGTTGTGGGCTGCGGGCGCTTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                     <  1:5567290/150‑1 (MQ=255)
   acAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGCGGCGGCGATTGgcc                                                                                                                                                    <  1:2343031/149‑1 (MQ=17)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                  <  2:2346803/150‑1 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                  >  1:336642/1‑150 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGACTGgccgc                                                                                                                                                  >  2:1076476/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                 >  2:3012319/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                                >  1:1584444/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                                >  1:178654/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCcccccc                                                                                                                                               >  1:2023595/1‑150 (MQ=12)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                               >  1:591525/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                               >  1:3876435/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                               >  1:486947/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                               <  2:274546/150‑1 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  2:3594410/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  1:3744190/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              <  1:3813234/150‑1 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  1:2545654/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  1:4061281/1‑150 (MQ=32)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCTCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCCGCGGCGCTTTGccccccccgg                                                                                                                                            >  1:5142089/1‑148 (MQ=11)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                            >  2:6669946/1‑150 (MQ=32)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTg                                                                                                                                            >  2:5038810/1‑150 (MQ=17)
               gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc                                                                                                                                        >  1:5325296/1‑150 (MQ=32)
                 tgtgggctgtgggcGAGTTTTACTGGGTTGTAACCTCGCGCCCTACGGGGTAGCCTTACGGGTTAAGTCGAGGGCACCGCGCAGCGGCTATTTTTTTGGCCTGTGCCGGGGGGTGCTGGGGGGGGCGGCGATGGGCCGCCCCGTgcgcgc                                                                                                                                      <  2:3634154/150‑1 (MQ=255)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  2:5585235/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTg                                                                                                                                >  1:5726181/1‑150 (MQ=17)
                            ggCGCTTTTTACCGGGATGTTCCCGCGCGCCCGACCGGGTCGCCTTACGGGATATGTGGCGGGCACCGCGCAGCGGCTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCa                                                                                                                           <  2:5535570/150‑1 (MQ=11)
                            ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCa                                                                                                                           >  1:208613/1‑150 (MQ=32)
                             gcggtttttACTGGTTTGTTACCTCGCGCCCGACCGGGTCGCCGTATGGTATATGTGGTGGGCACCGCGCCGGGGCTTTTTTGTTCGCCAGCGCCCGGGGGGGCCGGGGGCGGCGGCTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          <  2:3767465/145‑1 (MQ=2)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:557280/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAg                                                                                                                          >  1:3036550/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAgg                                                                                                                         >  1:4528875/1‑149 (MQ=17)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCa                                                                                                                     >  1:3368240/1‑150 (MQ=21)
                                    ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAat                                                                                                                   >  2:4246375/1‑150 (MQ=35)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  1:1479260/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:3105123/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCCTGGCCGCCCCCTGCGCGCTCCGTGCGCCCGGGGCAATATgg                                                                                                               >  1:1456644/1‑149 (MQ=11)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCATGCGCCAGTGGCAATATGt                                                                                                               >  1:2796747/1‑150 (MQ=255)
                                         ggtaTGTTACCTCGCGCCCGACCGGTTAGCCGTAAGGGATATTTCGTGGGCACCGCGCAGCGGCTATTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              <  2:4975864/147‑1 (MQ=14)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:4564756/1‑150 (MQ=40)
                                            aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGACAGGGGGTGCCGGGGGCGGCGGCGGTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGCAATATGTTGCt                                                                                                           >  1:429602/1‑150 (MQ=14)
                                              gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                         >  2:2615126/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGTGGGGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:1797080/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  1:3473551/1‑150 (MQ=255)
                                                 aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                      >  1:3595478/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTGGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCACTATTTTTGTTAtctc                                                                                                     >  1:1706353/1‑150 (MQ=11)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                                     >  2:1896718/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    >  2:5278311/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    >  2:1940879/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCCCCCCCAGCGCGCGCCATGCGCCCGGGGCAAGAAGTTGCGTAGcgccg                                                                                                   >  1:5603254/1‑146 (MQ=11)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                   >  2:6350489/1‑150 (MQ=255)
                                                     cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCACCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATg                                                                                                  >  1:6348202/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  2:682657/1‑150 (MQ=255)
                                                              ccggttAGCCGTATGGTATATGTCGAGGGCACCGCGCTGCGGCTTTTTTGTTCGCCTGCGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:6070022/145‑1 (MQ=17)
                                                              ccgggTAGCCTTATGGTATATGTGGTGGGCCCCGCGCGGGGGCTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCTTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:3744190/147‑1 (MQ=12)
                                                              cccGGTTGCCTTATGGTATACGTCGAGGGCACGGCGCTGCGGCGATTTTGTTCGCCTGAGCCCGGGGGTGCCGGGGGCGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:4389712/150‑1 (MQ=17)
                                                               ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGACCCCGGGGGTGCACGGCGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                        >  2:5840533/1‑150 (MQ=255)
                                                                      ccGTATGGGATAAGTCGAGGGCCCCGCGCCGGGGCGATTTTGTGCGCCAGAGCCGGGGGGTGCAGGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                 <  2:3866307/150‑1 (MQ=18)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                >  2:5424080/1‑150 (MQ=255)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCCCCCCCGGCGCGCGCCATGCGCCCGGGGCCAAATGTTGCTTAGCGCAAGAAAGGAGCGCAACAAcat                                                                                >  1:3696585/1‑150 (MQ=17)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTGCGCCCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                               >  1:4985133/1‑150 (MQ=255)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGGCCCCCCCCATCGCGCACCAATCGCCAATGGCAAAAATTTTCTTATCTCAAGAAACGAGCGCAACAagaaga                                                                              >  2:6170324/1‑150 (MQ=14)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATgaagaa                                                                          >  2:3220985/1‑150 (MQ=255)
                                                                               aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATTTTTGTTATCTCATGAAAGGAGCGCAACAAGATGTTGTAta                                                                        >  1:104591/1‑149 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTTGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:2309858/1‑150 (MQ=255)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCCGGGGTGCACGGGGCGGCGGCGAGTTGCCCCCCCCCGCGCGCTCCGTGCGCCCCTGGCAAAATGTTGCTTAGCGCGAGAAAGGGGCGCCACAAGATTATTGATCAACataa                                                                 >  1:871020/1‑149 (MQ=17)
                                                                                           cACGGCGCTGCGGGTTTTTTGTTGGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                            <  2:4458032/150‑1 (MQ=21)
                                                                                             ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAacaaca                                                          >  2:5988971/1‑150 (MQ=255)
                                                                                                       gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGTGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                >  2:3370168/1‑150 (MQ=255)
                                                                                                         ttttttGTTCGCCCGGGGGGGGGGGGGCGGGGGGGGGGGGGGGTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  1:237937/148‑1 (MQ=255)
                                                                                                          tttttGTTGCGCCTGCGCCGGGGGGGGCCGGGGGGGGGGGCTTTGGCCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  2:3711537/149‑1 (MQ=255)
                                                                                                            tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAAAAAGATGATGAATCAACATAAAACAACATCTTaaaaaaacgac                                           >  2:4707894/1‑146 (MQ=25)
                                                                                                                ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                       >  2:4393984/1‑150 (MQ=255)
                                                                                                                       agagCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                 >  2:3139002/1‑150 (MQ=255)
                                                                                                                          gcggggggggGCTGGGGGGGGGGGCGTGGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:5871347/146‑1 (MQ=255)
                                                                                                                          gCCGGGGGGTGCTGGGGGCGGCGGGGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:5992586/150‑1 (MQ=255)
                                                                                                                          ccgggggggTGCCGGGGGCGGCGGGTTTGGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  1:3254207/146‑1 (MQ=255)
                                                                                                                           ggggggggTGCGGGGGGGGGGGGCTTTTGCCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                             <  1:6068684/147‑1 (MQ=38)
                                                                                                                              gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          >  2:5669585/1‑150 (MQ=255)
                                                                                                                                   tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     >  1:5170431/1‑150 (MQ=255)
                                                                                                                                    gctggggggggCGGGGTTTGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                    <  2:1811499/147‑1 (MQ=38)
                                                                                                                                        ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  2:4093966/1‑150 (MQ=255)
                                                                                                                                        ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  2:6211384/1‑150 (MQ=255)
                                                                                                                                           ggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt             >  1:6524532/3‑150 (MQ=255)
                                                                                                                                            ggggggcgTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta            <  1:1096096/146‑1 (MQ=255)
                                                                                                                                             gggggcgATGGGCCGCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa           <  1:400584/147‑1 (MQ=255)
                                                                                                                                             ggcggccAATCGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTAGCTTAGAGCATGAAAGGAGCGCAACAAGAAGATGAATCAACATATAACAACATCTTAAAAAATGGCCTGACATTAAGCCATGACTTCTGAGttaattta           >  2:880532/1‑148 (MQ=37)
                                                                                                                                                ggggtttGGGCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc         <  2:1959892/145‑1 (MQ=255)
                                                                                                                                                 gcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc       >  1:2583142/1‑150 (MQ=255)
                                                                                                                                                 gcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc       >  1:2583084/1‑150 (MQ=255)
                                                                                                                                                  ggATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg      <  2:1553316/149‑1 (MQ=255)
                                                                                                                                                   gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  2:2268255/1‑150 (MQ=255)
                                                                                                                                                   gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  2:6561767/1‑150 (MQ=255)
                                                                                                                                                      tGGCCGCCCGATGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:4929043/150‑1 (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: