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% 300.9 / 2.8 119intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (79/27);  new base G (2/11);  total (81/38)
Fisher's exact test for biased strand distribution p-value = 5.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915066‑915363
                                                                                                                                                      |                                                                                                                                                     
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt‑t                                                                                                                                                       <  2:1887485/150‑1 (MQ=17)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt‑t                                                                                                                                                       <  2:112981/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                                      <  1:110212/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                                      <  2:3234682/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                                      <  2:1910374/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                                      <  2:1845883/150‑1 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgc                                                                                                                                                    <  1:217068/150‑1 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTCACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgc                                                                                                                                                    <  2:3577259/150‑1 (MQ=17)
     aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccg                                                                                                                                                  >  1:2139012/1‑150 (MQ=32)
      gCCCACCTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 <  2:185200/150‑1 (MQ=17)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 >  1:372529/1‑150 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 <  1:3249302/150‑1 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 >  1:2892370/1‑150 (MQ=32)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 >  2:3322037/1‑150 (MQ=17)
      gCCCACAGGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGGTACCGCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                                 <  1:1247452/150‑1 (MQ=255)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                                >  1:2887995/1‑150 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                                >  1:2888309/1‑150 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                                >  1:1598147/1‑150 (MQ=32)
         cacaGGGTTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGcccc                                                                                                                                              <  1:2092602/150‑1 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑Tgcg                                                                                                                                        >  2:53672/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTAAAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAA‑TGGCCGCCCCC‑Tgcg                                                                                                                                        >  2:53595/1‑150 (MQ=17)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTccat                                                                                                                                >  2:3012395/1‑148 (MQ=17)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTccgt                                                                                                                                >  1:3756487/1‑148 (MQ=17)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCtt                                                                                                                                >  2:1455653/1‑150 (MQ=17)
                                 tttttACTGGGTTGTTACCTCGCGCCCGACCGGGTAGCGTTATGGTATATGTCGTGGCCACCGCGCTGCGGCTTTTTTTTACGCAAGAGCCCGGGGGTGCGGGGGGCGGCGGCTTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGc                                                                                                                     <  2:579695/149‑1 (MQ=2)
                                 gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGc                                                                                                                     >  2:2579567/1‑150 (MQ=33)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                    >  1:2446567/1‑150 (MQ=34)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCATGCGCCAGTGGCa                                                                                                                    >  1:2739283/1‑150 (MQ=21)
                                   tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                                   >  2:775605/1‑150 (MQ=35)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  1:253428/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTACATGCGCCAGTGGCAATATGt                                                                                                              >  1:1897433/1‑150 (MQ=255)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGCGTT‑TTGCCGCCCCC‑TGGGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                            >  2:1373214/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCGGTTCGCGTTATGGTCAAAGTCGGGGGCACGGCGCCGCGGCTTTTTTTTTCGCCAGGGCCGGGGGGTGCCGGGGGGGTGGGCGTT‑TGGCCGCCCCC‑TCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                       <  2:1647926/150‑1 (MQ=9)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTTTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGGGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGGCAATAAGTTGCATAg                                                                                                       >  1:3402687/1‑150 (MQ=14)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TTGCCGCCCCC‑AGCCCGCACCCTGCGCCCGTGGCAATAATTTGCTTAg                                                                                                       >  1:3402647/1‑150 (MQ=14)
                                                tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                                      >  2:3735356/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                                    >  2:728090/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                   >  2:3228785/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                   >  2:371296/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                   >  2:371253/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                   >  2:371290/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCag                                                                                                  >  1:3879567/1‑149 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCAt                                                                                                  >  1:3879605/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTATCTCAt                                                                                                  >  1:3540203/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCGTGCGCCCGGGGCAATATGTTGCTTAGCTCAt                                                                                                  >  1:1812638/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  >  1:1812619/1‑150 (MQ=255)
                                                    ccgcgCCCCACCCGGTTGCGGTACGGTTTTCGTGGGGGCCACCGCGCTGGGGGTTTTTTTTTCGCCAGCGCCGGGGGGGGCTGGGGGCGGGGGGTTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                  <  2:3805790/149‑1 (MQ=12)
                                                     cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATg                                                                                                 >  2:1101644/1‑150 (MQ=255)
                                                     cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                 >  2:776079/1‑150 (MQ=255)
                                                         cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTTGCGCCTGAAAg                                                                                             >  1:3679929/1‑150 (MQ=14)
                                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                           >  1:2284225/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  2:3817314/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                          >  2:3050960/1‑150 (MQ=255)
                                                                 ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACCGGAAAGGAGCGcaa                                                                                     >  2:3697319/1‑150 (MQ=255)
                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                    >  2:3581451/1‑150 (MQ=255)
                                                                      ccttatgggATACGTCGAGGGCATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCGGGGGGTGCTGGGGGGGGGGGCGTT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                <  2:1334016/144‑1 (MQ=18)
                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCCAT‑TGGCCCCCCCC‑TGCGCGCGCCTTTCGCCACTGGCAATAATTTTCTTATCTCATTCAAGGAGCGCAACAAgaat                                                                              >  1:146865/1‑148 (MQ=17)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  2:357468/1‑150 (MQ=255)
                                                                         tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCGCCTTGCGCCACTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                             >  1:3816009/1‑150 (MQ=255)
                                                                           aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGGCGCCCCC‑AGCGCGCGCCATGCGCCAGTGGCAAAATTTTTCTTAGCTCAAGAAAGGAGCGCAACAAgatgatt                                                                           >  1:3776378/1‑149 (MQ=17)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                                                         >  2:3529613/1‑150 (MQ=255)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCG‑CCCC‑TGCGCGCTCCTTGCCCCAGTGGCAAAAAGTTTCTTATCTCATGAAAGGAGCGCAAAAAGAATATGAAt                                                                        >  2:2835342/1‑150 (MQ=17)
                                                                                 acgccgtgggCACCGCGCCGCGGCTTTTTTTTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGTT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                     <  2:1997124/143‑1 (MQ=18)
                                                                                 aaGTCGTGGGCACCGCGCAGCGGCGTTTTTGTTCGCCAGAGCCGGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                     <  2:2281277/150‑1 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑AGCGCCCGCCCTGCGCCCCTGGCAAAATGTTTCATATCTCAGGGAAGGAGCGCAACAAGAAGAATAATAAACaa                                                                  >  1:3529964/1‑149 (MQ=14)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                  >  2:2613928/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCCCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACat                                                                  >  2:3015331/1‑150 (MQ=255)
                                                                                            aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGAGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCCACAAGATGATGAATCAACAAATaacaac                                                          >  2:813836/1‑150 (MQ=255)
                                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        >  2:1735789/1‑150 (MQ=255)
                                                                                               gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                       >  2:807614/1‑150 (MQ=255)
                                                                                                   cgcggcgTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                                   <  2:3956706/149‑1 (MQ=255)
                                                                                                    ggggcttttttGTTCGCCTGAGCGGGGGGGGGGGGGGGGGGGCGGCTTT‑GGCCCGCCCCC‑CGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                  <  1:3289126/143‑1 (MQ=255)
                                                                                                       gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACATATAACAACATCTTaaaaaa                                               >  2:928315/1‑150 (MQ=255)
                                                                                                        tttttttGTTCGCCAGTGCGGGGGGGTGCGGGGGGGGGGGGCTTT‑GGCCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                              <  2:2759414/147‑1 (MQ=255)
                                                                                                        cttttttGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCGAT‑TGGCGGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                              <  2:2130842/147‑1 (MQ=25)
                                                                                                         ttttttGTCCGCCGGCGCGGGGGGGGGGGGGGGGCGGGGGGGTG‑GGCCCCCCCCC‑CCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             <  1:25731/148‑1 (MQ=255)
                                                                                                         gtttttGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             <  2:2972112/148‑1 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                             >  2:122581/1‑150 (MQ=255)
                                                                                                          tttttGTTCGCCTGTGCCGGGGGGTGCTGGGGGGGGGGGGTTT‑TGGCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                            <  2:3491647/149‑1 (MQ=25)
                                                                                                           ttttttTCGCCTGCGCCGGGGGGTGCGGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                           <  2:3920946/150‑1 (MQ=25)
                                                                                                           ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGTGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                           <  2:883642/150‑1 (MQ=255)
                                                                                                             ttGTTCGCCAGAGCCCGGGGGGGCCGGGGGGGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                                         <  2:3946399/150‑1 (MQ=255)
                                                                                                                  cGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                                    >  2:651283/1‑150 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAAATCTTAAAAAAAGGCCTGAAATTa                                  >  2:1469379/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑GGGCCGCCCCC‑TGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCGTAGCTCATGACAGGAGCGCAACAAGATGAAGCAGCAACATATCACAACATCTTCCAAAAAGGCCTGAccgtcc                                 >  2:3519878/1‑145 (MQ=18)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCC‑TGCGCGCGCCTTGCGCCCGGGGCAAGATGGTGCGTCGCGCCGGAAAGGAGCGCACCAAGATGAGGACTCAACCTAGAACACCCGCTTACAAAAAGGCCTGACagtcc                                 >  2:3519876/1‑146 (MQ=12)
                                                                                                                        ttgcggggggTGTGGGGGGGGGGGGGGTTGTGGCCCCCCCG‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                               <  2:3402647/145‑1 (MQ=255)
                                                                                                                           ccGGGTGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                           >  2:2610340/1‑150 (MQ=255)
                                                                                                                                   gggggggggggggGCGGT‑GGGCCCCCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATGTAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCttttt                    <  2:1159784/147‑4 (MQ=255)
                                                                                                                                      ggggggggCGGCGTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  2:2887995/150‑1 (MQ=255)
                                                                                                                                      gggggCGGCGGCTAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                <  2:1598147/150‑1 (MQ=255)
                                                                                                                                       ggggCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta               >  1:3333389/1‑150 (MQ=255)
                                                                                                                                        gggCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa              >  1:3331629/1‑150 (MQ=255)
                                                                                                                                          ggggcggccAA‑TCGCCGCCGCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            >  2:251787/3‑150 (MQ=255)
                                                                                                                                          gcgggggcgAT‑TGGCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt            <  2:1288777/150‑1 (MQ=255)
                                                                                                                                                 gttt‑tGGCCGCCCCG‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg     <  2:3816009/147‑1 (MQ=255)
                                                                                                                                                    t‑tGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  2:2567022/1‑150 (MQ=255)
                                                                                                                                                    t‑tGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  2:1611465/1‑150 (MQ=255)
                                                                                                                                                    t‑tGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  1:1463559/1‑150 (MQ=255)
                                                                                                                                                    t‑tGGCCGCCCAC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAATGGCCTGACATTACTCCAGGCCTTCTGCGTTAATTAATCACGCtt  >  2:1920430/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915066‑915363

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