Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,213 A→T 6.8% intergenic (‑356/+139) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2130AT 6.8% 503.9 / 5.0 213intergenic (‑356/+139)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base A (140/40);  new base T (3/13);  total (160/53)
Fisher's exact test for biased strand distribution p-value = 3.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCT  >  NC_000913/915064‑915362
                                                                                                                                                     |                                                                                                                                                     
ccGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGa                                                                                                                                                       <  1:3054299/150‑1 (MQ=32)
  ccctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGGCAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGagt                                                                                                                                                     >  2:4161386/5‑148 (MQ=11)
    ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                   >  1:1338594/1‑150 (MQ=32)
    ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                   <  2:4736706/150‑1 (MQ=32)
     tcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                  >  2:5780712/2‑150 (MQ=17)
     gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                  >  2:4696129/1‑150 (MQ=32)
      cAGCCCACATGGATGTGGGCTGCGGGCGCGTTTTACAGGGATTTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATTAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGTGCGGCGGCGATTGgcc                                                                                                                                                 <  1:3240240/150‑1 (MQ=255)
      cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc                                                                                                                                                 <  1:1570228/150‑1 (MQ=32)
      cAGCCCACAAGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTTgcc                                                                                                                                                 >  1:702932/1‑150 (MQ=12)
          ccACAGGGATGGGGGCTGAGGGCGCGTTTTACAGGGATGTTACCGCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                             <  2:2267191/150‑1 (MQ=12)
           taacgcagaaggcctggcgtaatgtcaggcctttttttaagatgttgttatatgttgattcatcatcttgttgcgctcctttcatgagctaagcaacatattgccactggcgcaaggagcgcGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                            >  1:630594‑M1/123‑150 (MQ=255)
           cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                            >  1:3065602/1‑150 (MQ=32)
            acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGTATAATTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGACAGAGCCCGGGGGTGCATGGGGCGGCGGCGATAGGCCGccccc                                                                                                                                           >  1:567682/1‑150 (MQ=11)
            acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                           >  1:1485157/1‑150 (MQ=32)
            acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                           >  1:1485076/1‑150 (MQ=32)
              cgcagaaggcctggcgtaatgtcaggcctttttttaagatgttgttatttgttgattcatcatcttgttgcgctccttacatgagctaagcaacatattgccactggcgcatggagcgtGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                         >  1:908868‑M1/120‑150 (MQ=255)
              cgcagaaggcctggcgtaaggtcaggcctttttttaagaggttgttatatgttgattcctcatcttgttgcgctcctttcatgagctaagcaacatattgccactggcgcaaggagcgcGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                         >  1:4649856‑M1/120‑150 (MQ=255)
               tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgc                                                                                                                                        <  2:4675102/150‑1 (MQ=32)
                ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTgcg                                                                                                                                       >  1:2342932/3‑150 (MQ=17)
                  aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGgcgcg                                                                                                                                     >  2:1230725/1‑150 (MQ=17)
                   ggggggCTGCGGGCGCGTTTTACAGGGATGTGACCTCGCGCCCGACCCGGTAGCGGTACGGGATAAGTCGCGGGCACCGCGCAGCGGCGTTTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGGGATTGGCCGCCCCCTgcgcgc                                                                                                                                    <  2:5178282/147‑1 (MQ=11)
                      gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTcc                                                                                                                                 >  1:5156831/1‑150 (MQ=32)
                      gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTcc                                                                                                                                 >  1:4043402/1‑150 (MQ=32)
                      gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTcc                                                                                                                                 >  1:1191177/1‑150 (MQ=17)
                         cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                              >  2:3530183/1‑150 (MQ=32)
                         cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                              >  2:4669779/1‑150 (MQ=32)
                         cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                              >  2:4669723/1‑150 (MQ=17)
                                 gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg                                                                                                                      >  2:3352655/1‑150 (MQ=32)
                                 gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg                                                                                                                      >  2:3672181/1‑150 (MQ=32)
                                  ctttttACTGGGATGTTACCTCGCGCCCTACCGGTTAGCGTTATGGTATATGTCTTGGGCACGGCGACGCGGGGTTTTTGTTGGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTgg                                                                                                                     <  2:451370/148‑1 (MQ=2)
                                    ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                   >  2:1386771/1‑150 (MQ=34)
                                    ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                   >  2:2583769/1‑150 (MQ=34)
                                    aattACACGGATCTTACCTCGCGCCCGAACCGGTAGCCGGAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTGTGTTCGCCAGAGCCCGGGGGTGCAGCGGGCGGCGGCGATTCGCCCCCACCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                   >  2:5280497/3‑150 (MQ=11)
                                     tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                                  >  1:5818897/1‑150 (MQ=35)
                                       tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCaaaa                                                                                                                >  2:1614077/1‑148 (MQ=25)
                                        aCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAatat                                                                                                               >  2:5258257/1‑150 (MQ=38)
                                        aCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCAGTGGCAatat                                                                                                               >  1:5089329/1‑150 (MQ=255)
                                          aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                             >  2:4950799/1‑150 (MQ=39)
                                          aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                             >  2:2605459/1‑150 (MQ=39)
                                              aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                         >  2:2374377/1‑150 (MQ=255)
                                              aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                         >  1:2435849/1‑150 (MQ=255)
                                                   cccTCGCGCCGGACCCGGTTGCCTTATGGGATACGTCGTGGGCACCGCGCCGGGGGTTTTTTTTTGGCCAGAGCGGGGGGGTGCTGGGGGGGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                    <  2:422056/149‑1 (MQ=9)
                                                    ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCATTCGCCAGTGGCAAAATGTTGCTTTGCTc                                                                                                   >  2:2543401/1‑150 (MQ=14)
                                                     cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                  >  2:194642/1‑150 (MQ=255)
                                                     cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                  >  2:539568/1‑150 (MQ=255)
                                                       cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCCCCACTGGCAATATGTTGCTTAGCGCATg                                                                                                >  2:1030309/1‑150 (MQ=14)
                                                        gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGa                                                                                               >  1:3789109/1‑150 (MQ=255)
                                                          gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaaa                                                                                             >  1:1353115/1‑150 (MQ=255)
                                                           cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCCGGGGTGCACGGGGCGGCGGCCCATGGCCGCCCCCAGCCCCCCCCGTGCGCCCGGGGCCATATGTTGCTTAGCTCATGAAAg                                                                                            >  1:202233/1‑150 (MQ=17)
                                                            ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAgg                                                                                           >  1:5176823/1‑150 (MQ=255)
                                                            ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                           >  1:3685519/1‑150 (MQ=255)
                                                            ccGAACCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGGGGGTGGCCGCCCCCTGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGGTCATGattgg                                                                                           >  1:870286/1‑146 (MQ=17)
                                                              gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTGGGCCGCCCCCAGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCTTAGCGCCTGAAAggcg                                                                                         >  1:2738013/1‑148 (MQ=17)
                                                                 ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCCCCAGTGGCAATAAGTTGCTTAGCTCAAGAAAGGAgcgc                                                                                      >  2:270221/1‑150 (MQ=255)
                                                                   ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGcaa                                                                                    >  1:3791181/1‑150 (MQ=255)
                                                                      cgcggTATGGGATCCGACGCGGGCCCCGCCCCGCGGCGATTTGGTGGCCCAGCGCCGGGGGGGGCCGGGGGGGGCGGCGTTGGCCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaacaa                                                                                 <  2:3611864/146‑1 (MQ=255)
                                                                         cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCATTCGCCAGTGGCAATATGTTGCATAGCTCAAGAAAGGAGCGCAACAAgat                                                                              >  1:3454108/1‑150 (MQ=17)
                                                                          gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCCCCCCCAGCGCCCACCATTCCCCACTGGCAAAATTTTTCATAACACAAGAAACGAGCGcaacaaaaat                                                                             >  1:3917530/1‑146 (MQ=14)
                                                                              gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCACATGAAAGGAGCGCAACAAGAtgatga                                                                         >  2:4811914/1‑150 (MQ=255)
                                                                                gaaaaGCCGAGCCCACCGCCCAGCGGCGATTTTGTACGCCAGACCCCGGGGGTCCAGGGGGCGCCGGCGATTGCCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                       <  2:4344606/147‑1 (MQ=18)
                                                                                gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGGTGGCCGCCCCCAGCGCGCGCCATGCGCCCGTGGCAATATGTTGCGTAGCGCATGTAAGGAGCGCAACAAGAAGATGAAt                                                                       >  1:3530332/1‑150 (MQ=18)
                                                                                 aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAATATTTTGCTTATCTCAAGAAAGGAGCGCAACAAGATGATGAAta                                                                      >  1:3201331/1‑149 (MQ=18)
                                                                                 aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATc                                                                      >  2:5696240/1‑150 (MQ=255)
                                                                                 aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCCGTGGCAATATGTTGCTTATCTCATGTAAGGAGCGCAACAAGATGATGAATc                                                                      >  1:3630359/1‑150 (MQ=18)
                                                                                  tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCa                                                                     >  2:1490426/1‑150 (MQ=255)
                                                                                      tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCCCCCCCAGCGCGCTCCATGCGCCCCTGGCAAAATGTTGCTTATCTCAAGAAAGGAGCGCACCAAGAATATTAATCAACaa                                                                 >  1:4510885/1‑149 (MQ=17)
                                                                                         aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCCGCGGCGGGTGGCCCCCCCCCGCGCGCGCCGTGCGCCCGGGGCAAAAAGTTTCGTAGCTCATGGAAGGAGCGCACCAAGATGATGCATCAACAgata                                                              >  1:337394/1‑150 (MQ=17)
                                                                                              aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGGGTGGCCCCCCCCCGCGCCCGCCATGCGCCCCGGGCCAAAAGTTGCTTAGCGCCTGACAGGAGCGCCACCACATGGTGAATCAACAAATaacaac                                                         >  1:4345666/1‑150 (MQ=12)
                                                                                                 gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                      >  2:4062726/1‑150 (MQ=255)
                                                                                                    ccgcggcgTTTTTGTGCGCCTGTGCCGGGGGGTGCAGGGGGCGGCGGCGTTGGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTa                                                   <  2:1866126/148‑1 (MQ=25)
                                                                                                      agcaacatattgccactggcgcaaggagcgcGCAGGGGGCGGCGGCGATTGGCCAcccc                                                                                                                                            >  2:1462781‑M1/32‑59 (MQ=255)
                                                                                                            aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                           >  2:606737/1‑150 (MQ=255)
                                                                                                              tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                         <  2:3833261/150‑1 (MQ=255)
                                                                                                                 gTTCGCCCGCGCCGGGGGGTGGGGGGGGCGGCGGCGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                      <  2:3980495/150‑1 (MQ=37)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                  >  2:762468/1‑150 (MQ=255)
                                                                                                                      ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGACATTa                                 >  2:2080258/1‑150 (MQ=255)
                                                                                                                       cAGAGCCCGGGGGTGCAGGGGGCGGCGGGTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                <  2:3339336/150‑1 (MQ=255)
                                                                                                                            cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                           >  2:5582859/1‑150 (MQ=255)
                                                                                                                              ggggggTGGTGGGGGGGGGGGTGTTTGGCCGCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                         <  2:1600630/149‑1 (MQ=255)
                                                                                                                              ggggggTGCTGGGGGGGGCGGCGTTTGGCCGCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                         <  2:1820272/149‑1 (MQ=255)
                                                                                                                              cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATAAGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                         >  2:2042125/1‑150 (MQ=255)
                                                                                                                               gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTATCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                        >  2:1159213/1‑150 (MQ=255)
                                                                                                                               gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                        >  2:2431006/1‑150 (MQ=255)
                                                                                                                                  ggggctggggggggCGGGGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                     <  2:5123390/144‑1 (MQ=255)
                                                                                                                                     cgcgcaGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                  <  2:1039643‑M1/144‑1 (MQ=255)
                                                                                                                                       AGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                <  1:3148491‑M1/150‑1 (MQ=255)
                                                                                                                                         GGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta              <  2:4115513‑M1/150‑1 (MQ=255)
                                                                                                                                          GGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTATAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa             <  2:3000983‑M1/150‑1 (MQ=255)
                                                                                                                                          GGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa             <  2:485131‑M1/150‑1 (MQ=255)
                                                                                                                                           GGCGGCGGGGTTTGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat            >  2:870286‑M1/1‑150 (MQ=255)
                                                                                                                                               GCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCCGGCCTTCTGCGTtaattaat        <  1:2132762‑M1/150‑1 (MQ=255)
                                                                                                                                               GCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCgttaattaag        <  1:3584497‑M1/150‑2 (MQ=255)
                                                                                                                                               GCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat        <  1:5537060‑M1/150‑1 (MQ=255)
                                                                                                                                               GCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat        <  2:800919‑M1/150‑1 (MQ=255)
                                                                                                                                                 ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      >  2:3251307/1‑150 (MQ=255)
                                                                                                                                                 ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      >  2:4327786/1‑150 (MQ=255)
                                                                                                                                                 ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      >  2:3795111/1‑150 (MQ=255)
                                                                                                                                                    gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc   >  2:2307562/1‑150 (MQ=255)
                                                                                                                                                     aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt  >  2:2394514/1‑150 (MQ=255)
                                                                                                                                                     aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt  >  2:2530138/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCT  >  NC_000913/915064‑915362

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