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.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                      |                                                                                                                                                       
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                         <  1:1287226/150‑1 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                         <  1:40935/150‑1 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       <  1:5790438/150‑1 (MQ=32)
  ctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  1:1198870/3‑150 (MQ=12)
  ctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACATAGCGCCCGACCCGGAAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑AGCCAGATCCCGGGGGTGCAGGGGGCGGCCGCGATTgg                                                                                                                                                       >  1:4562051/3‑150 (MQ=11)
  ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  2:1023937/3‑150 (MQ=12)
  ctcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                       >  2:1181293/3‑150 (MQ=12)
   gCAGCCCCCATGGTTGTGGGCTGCGGGCGCTTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑c                                                                                                                                                     <  1:5567290/150‑1 (MQ=255)
   acAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCA‑GGGGCGGCGGCGATTGg‑cc                                                                                                                                                    <  1:2343031/149‑1 (MQ=17)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgc                                                                                                                                                  <  2:2346803/150‑1 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgc                                                                                                                                                  >  1:336642/1‑150 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGACTGg‑ccgc                                                                                                                                                  >  2:1076476/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgcc                                                                                                                                                 >  2:3012319/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccc                                                                                                                                                >  1:1584444/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccc                                                                                                                                                >  1:178654/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGA‑TTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑Ccccccc                                                                                                                                               >  1:2023595/1‑150 (MQ=12)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                               >  1:591525/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                               >  1:3876435/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                               >  1:486947/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                               <  2:274546/150‑1 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccccc                                                                                                                                              >  2:3594410/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccccc                                                                                                                                              >  1:3744190/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccccc                                                                                                                                              <  1:3813234/150‑1 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccccc                                                                                                                                              >  1:2545654/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccccc                                                                                                                                              >  1:4061281/1‑150 (MQ=32)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCTCGCAGCGGCGATTTTGTT‑CGCCCGAGCCCGGGGGTGCAGGGGGCCGCGGCGCTTTG‑ccccccccgg                                                                                                                                            >  1:5142089/1‑148 (MQ=11)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTg                                                                                                                                            >  2:6669946/1‑150 (MQ=32)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCTg                                                                                                                                            >  2:5038810/1‑150 (MQ=17)
               gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTg‑cgc                                                                                                                                        >  1:5325296/1‑150 (MQ=32)
                 tgtgggctgtgggcGAGTTTTACTGGGTTGTAACCTCGCGCCCTACGGGGTAGCCTTACGGGTTAAGTCGAGGGCACCGCGCAGCGGCTATTTTTTT‑GGCCTGTGCCGGGGGGTGCTGGGGGGGGCGGCGATGGG‑CCGCCCCGTg‑cgcgc                                                                                                                                      <  2:3634154/150‑1 (MQ=255)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCt                                                                                                                                  >  2:5585235/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCTTg                                                                                                                                >  1:5726181/1‑150 (MQ=17)
                            ggCGCTTTTTACCGGGATGTTCCCGCGCGCCCGACCGGGTCGCCTTACGGGATATGTGGCGGGCACCGCGCAGCGGCTTTTTTGTT‑CGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCGTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCa                                                                                                                           <  2:5535570/150‑1 (MQ=11)
                            ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCa                                                                                                                           >  1:208613/1‑150 (MQ=32)
                             gcggtttttACTGGTTTGTTACCTCGCGCCCGACCGGGTCGCCGTATGGTATATGTGGTGGGCACCGCGCCGGGGCTTTTTTGTT‑CGCCAGCGCCCGGGGGGGCCGGGGGCGGCGGCTATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                          <  2:3767465/145‑1 (MQ=2)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                          >  2:557280/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCGTGCGCCAg                                                                                                                          >  1:3036550/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAgg                                                                                                                         >  1:4528875/1‑149 (MQ=17)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCa                                                                                                                     >  1:3368240/1‑150 (MQ=21)
                                    ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAat                                                                                                                   >  2:4246375/1‑150 (MQ=35)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  1:1479260/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:3105123/1‑150 (MQ=255)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCCTGG‑CCGCCCCCTG‑CGCGCTCCGTGCGCCCGGGGCAATATgg                                                                                                               >  1:1456644/1‑149 (MQ=11)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCATGCGCCAGTGGCAATATGt                                                                                                               >  1:2796747/1‑150 (MQ=255)
                                         ggtaTGTTACCTCGCGCCCGACCGGTTAGCCGTAAGGGATATTTCGTGGGCACCGCGCAGCGGCTATTTTGTT‑CGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              <  2:4975864/147‑1 (MQ=14)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:4564756/1‑150 (MQ=40)
                                            aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGACAGGGGGTGCCGGGGGCGGCGGCGGTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGGGGCAATATGTTGCt                                                                                                           >  1:429602/1‑150 (MQ=14)
                                              gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                         >  2:2615126/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGTGGGGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:1797080/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  1:3473551/1‑150 (MQ=255)
                                                 aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                      >  1:3595478/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTGGG‑CCGCCCCCAG‑CGCGCTCCGTGCGCCAGTGGCACTATTTTTGTTAtctc                                                                                                     >  1:1706353/1‑150 (MQ=11)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                                     >  2:1896718/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    >  2:5278311/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    >  2:1940879/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGG‑CCCCCCCCAG‑CGCGCGCCATGCGCCCGGGGCAAGAAGTTGCGTAGcgccg                                                                                                   >  1:5603254/1‑146 (MQ=11)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                   >  2:6350489/1‑150 (MQ=255)
                                                     cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCACCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATg                                                                                                  >  1:6348202/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  2:682657/1‑150 (MQ=255)
                                                              ccggttAGCCGTATGGTATATGTCGAGGGCACCGCGCTGCGGCTTTTTTGTT‑CGCCTGCGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:6070022/145‑1 (MQ=17)
                                                              ccgggTAGCCTTATGGTATATGTGGTGGGCCCCGCGCGGGGGCTTTTTTGTT‑CGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCTTTTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:3744190/147‑1 (MQ=12)
                                                              cccGGTTGCCTTATGGTATACGTCGAGGGCACGGCGCTGCGGCGATTTTGTT‑CGCCTGAGCCCGGGGGTGCCGGGGGCGGCGGCTTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                         <  2:4389712/150‑1 (MQ=17)
                                                               ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGACCCCGGGGGTGCACGGCGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                        >  2:5840533/1‑150 (MQ=255)
                                                                      ccGTATGGGATAAGTCGAGGGCCCCGCGCCGGGGCGATTTTGTG‑CGCCAGAGCCGGGGGGTGCAGGGGGCGGGGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                 <  2:3866307/150‑1 (MQ=18)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                >  2:5424080/1‑150 (MQ=255)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCCCCCCCGG‑CGCGCGCCATGCGCCCGGGGCCAAATGTTGCTTAGCGCAAGAAAGGAGCGCAACAAcat                                                                                >  1:3696585/1‑150 (MQ=17)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTG‑CCGCCCCCTG‑CGCCCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                               >  1:4985133/1‑150 (MQ=255)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGG‑CCCCCCCCAT‑CGCGCACCAATCGCCAATGGCAAAAATTTTCTTATCTCAAGAAACGAGCGCAACAagaaga                                                                              >  2:6170324/1‑150 (MQ=14)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATgaagaa                                                                          >  2:3220985/1‑150 (MQ=255)
                                                                               aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGG‑CCGCCCCCAG‑CGCGCTCCGTGCGCCAGTGGCAATATTTTTGTTATCTCATGAAAGGAGCGCAACAAGATGTTGTAta                                                                        >  1:104591/1‑149 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTTG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:2309858/1‑150 (MQ=255)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCCGGGGTGCACGGGGCGGCGGCGAGTTG‑CCCCCCCCCG‑CGCGCTCCGTGCGCCCCTGGCAAAATGTTGCTTAGCGCGAGAAAGGGGCGCCACAAGATTATTGATCAACataa                                                                 >  1:871020/1‑149 (MQ=17)
                                                                                           cACGGCGCTGCGGGTTTTTTGTT‑GGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGGGTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                            <  2:4458032/150‑1 (MQ=21)
                                                                                             ccGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAacaaca                                                          >  2:5988971/1‑150 (MQ=255)
                                                                                                       gcgATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGTGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                >  2:3370168/1‑150 (MQ=255)
                                                                                                         ttttttGTT‑CGCCCGGGGGGGGGGGGGCGGGGGGGGGGGGGGGTGGG‑CCCCCCCCGC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  1:237937/148‑1 (MQ=255)
                                                                                                          tttttGTTGCGCCTGCGCCGGGGGGGGCCGGGGGGGGGGGCTTTGGC‑CCCCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  2:3711537/149‑1 (MQ=255)
                                                                                                            tttGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCCCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAAAAAGATGATGAATCAACATAAAACAACATCTTaaaaaaacgac                                           >  2:4707894/1‑146 (MQ=25)
                                                                                                                tt‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                       >  2:4393984/1‑150 (MQ=255)
                                                                                                                       agagCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                 >  2:3139002/1‑150 (MQ=255)
                                                                                                                          gcggggggggGCTGGGGGGGGGGGCGTGGGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:5871347/146‑1 (MQ=255)
                                                                                                                          gCCGGGGGGTGCTGGGGGCGGCGGGGTTGGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:5992586/150‑1 (MQ=255)
                                                                                                                          ccgggggggTGCCGGGGGCGGCGG‑GTTTGGCCCCCCCCCTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  1:3254207/146‑1 (MQ=255)
                                                                                                                           ggggggggTGCGGGGGGGGGGGGCTTTTGC‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                             <  1:6068684/147‑1 (MQ=38)
                                                                                                                              gggggTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          >  2:5669585/1‑150 (MQ=255)
                                                                                                                                   tGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     >  1:5170431/1‑150 (MQ=255)
                                                                                                                                    gctggggggggCGGGGTTTGC‑CCCCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                    <  2:1811499/147‑1 (MQ=38)
                                                                                                                                        ggggCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  2:4093966/1‑150 (MQ=255)
                                                                                                                                        ggggCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  2:6211384/1‑150 (MQ=255)
                                                                                                                                           ggggcggccAATCG‑CCGCCGCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt             >  1:6524532/3‑150 (MQ=255)
                                                                                                                                            ggggggcgTTGGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta            <  1:1096096/146‑1 (MQ=255)
                                                                                                                                             gggggcgATGGG‑CCGCCCCCCG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa           <  1:400584/147‑1 (MQ=255)
                                                                                                                                             ggcggccAATCG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTAGCTTAGAGCATGAAAGGAGCGCAACAAGAAGATGAATCAACATATAACAACATCTTAAAAAATGGCCTGACATTAAGCCATGACTTCTGAGttaattta           >  2:880532/1‑148 (MQ=37)
                                                                                                                                                ggggtttGG‑GCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc         <  2:1959892/145‑1 (MQ=255)
                                                                                                                                                 gcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc       >  1:2583142/1‑150 (MQ=255)
                                                                                                                                                 gcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc       >  1:2583084/1‑150 (MQ=255)
                                                                                                                                                  ggATGGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg      <  2:1553316/149‑1 (MQ=255)
                                                                                                                                                   gATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  2:2268255/1‑150 (MQ=255)
                                                                                                                                                   gATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  2:6561767/1‑150 (MQ=255)
                                                                                                                                                      tGG‑CCGCCCGATG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:4929043/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTT‑CGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  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: ‑