Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 9.9% intergenic (‑369/+126) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2260TG 9.9% 430.4 / 7.3 218intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (135/25);  new base G (13/13);  total (171/47)
Fisher's exact test for biased strand distribution p-value = 2.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  NC_000913/915081‑915374
                                                                                                                                                  |                                                                                                                                                     
gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTg‑cgc                                                                                                                                                  >  2:6174002/1‑150 (MQ=32)
     gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCTG‑CGCGCTcc                                                                                                                                             >  2:3387269/1‑150 (MQ=17)
       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCtt                                                                                                                                           >  2:4573366/1‑150 (MQ=32)
       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCtt                                                                                                                                           >  2:3776684/1‑150 (MQ=32)
        cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTg                                                                                                                                          >  1:2548298/1‑150 (MQ=32)
        cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCTTg                                                                                                                                          >  1:3952556/1‑150 (MQ=17)
         tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCGCCTTgc                                                                                                                                         >  1:1684144/1‑150 (MQ=17)
             ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCa                                                                                                                                     >  1:649831/1‑150 (MQ=17)
              gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                                    >  2:3354892/1‑150 (MQ=32)
              gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                                    >  1:3244632/1‑150 (MQ=32)
               cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCCGt                                                                                                                                   >  1:5975852/1‑150 (MQ=17)
                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGTACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCATGCGCCAGTg                                                                                                                                  >  1:5668236/1‑150 (MQ=12)
                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTg                                                                                                                                  >  2:3743109/1‑150 (MQ=32)
                   ttttCCGGGGTTGTTACCCCCCGCCCGACCGGGGAGCGGTCCGGTATAGGTGGCGGGCACCGCGCCGGGGGGTTTTTGGCCGCAACACCGGGGGGGGGCTGGGGGGGGGGGCGTTGGC‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCa                                                                                                                               <  2:728343/150‑1 (MQ=2)
                        cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCCGTGGCAATATg                                                                                                                          >  1:4030519/1‑150 (MQ=38)
                        cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCGTGCGCCAGTGGCAATATg                                                                                                                          >  1:2132890/1‑150 (MQ=255)
                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAatattt                                                                                                                         >  1:1293421/1‑148 (MQ=255)
                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCAATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                         >  2:6197812/1‑150 (MQ=255)
                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCAAATGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAAAATGt                                                                                                                         >  2:2902257/1‑150 (MQ=255)
                          gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                                        >  1:2220556/1‑150 (MQ=40)
                           gtttGTTCCCGCGCGCCCGCCCGGGTAGCCGTATGGTATATGTCGCGGGCACCGCGCCGCGGCGTTTTTGTGCGCCAGCGCCGGGGGGTGCCGGGGGGGGGGGCGTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                       <  2:3592630/147‑1 (MQ=14)
                           ggttGTTCCCGCCCGCCGGACCGGGTCGCCGAACGGTTAAAGTGGTGGGCACCGGGCCGGGGTTTTTTTGTTCGCTAGCGCCGGGGGGTGCTGGGGGGGGCGGCTTTGGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                       <  2:5839029/147‑1 (MQ=9)
                           ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCATGGGGCGGCGGCGAGTGG‑CCCCCCCCAG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTg                                                                                                                       >  1:4943768/1‑150 (MQ=14)
                           ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGGGGCAATATGTTg                                                                                                                       >  1:880363/1‑150 (MQ=255)
                           ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTGGG‑CCGCCCCCGG‑CGCGCGCCTTGCGCCAGTGGCAATATGTTg                                                                                                                       >  1:1316828/1‑150 (MQ=14)
                            gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGc                                                                                                                      >  2:566794/1‑150 (MQ=255)
                                ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCACTATGTTGCTTAg                                                                                                                  >  1:4157176/1‑150 (MQ=255)
                                  aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGGGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCGTGCGCCAGTGGGCAGATGTTGGttttct                                                                                                                >  1:1132445/1‑146 (MQ=11)
                                  aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCATGCGCCAGTGGCAATATGTTTCTTAGCt                                                                                                                >  1:3601414/1‑150 (MQ=255)
                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                              >  1:5237294/1‑150 (MQ=255)
                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                             >  2:1576555/1‑150 (MQ=255)
                                       gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTATCTCATGa                                                                                                           >  2:2642085/1‑150 (MQ=255)
                                       gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGa                                                                                                           >  1:4965747/1‑150 (MQ=255)
                                        cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaa                                                                                                          >  2:651046/1‑150 (MQ=255)
                                             gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCACGGGCAATATGTTGCTTAGCTCAAGAAAGGAg                                                                                                     >  2:587502/1‑150 (MQ=255)
                                             gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCGCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAg                                                                                                     >  1:2992744/1‑150 (MQ=255)
                                                ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCGCCGTGCGCCCCTGGCCAAAATTTGCTTATCTCAAGAAAGGAgcgc                                                                                                  >  1:1148194/1‑150 (MQ=14)
                                                 cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCACTGGCAATAAGTTGCTTAGCTCATGAAAGGAGCGca                                                                                                 >  1:5786167/1‑150 (MQ=255)
                                                   gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                               >  1:4326344/1‑150 (MQ=255)
                                                   gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                               >  2:2564632/1‑150 (MQ=255)
                                                    tAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGATCCCGGGGGTGCAGGGGGCGGCGGCGTGTGG‑CCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAATGAGCGcaaca                                                                                              >  2:5127884/1‑150 (MQ=17)
                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTG‑CCCCCCCCGG‑CGCGCACCGTGCGGCACTGGCAAAATTTTTCATAGCTCATGAAAGGAGCGCAACAAgaag                                                                                         >  1:4392256/1‑148 (MQ=17)
                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                                         >  2:2205612/1‑150 (MQ=255)
                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATTG‑CCGCCCCCAG‑CGCGCTCCTTGCCCCAGTGGCAAAATTTTGCATAGCACAAGAAAGGAGCGCAACAAtatg                                                                                         >  1:4391958/1‑150 (MQ=17)
                                                                aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATc                                                                                  >  2:1224912/1‑150 (MQ=255)
                                                                aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCCGTGGCCAAATTGTGCTTATGGCATGAAAGGGGCGCCACAAGATGATGAATc                                                                                  >  1:155231/1‑150 (MQ=17)
                                                                 tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCa                                                                                 >  2:2012165/1‑150 (MQ=255)
                                                                 tAAGTCGAGGGCACCGCGCAGCGGCGATTGTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCGTGCGCCAGGGGCACGATGTTGCTTAGCGCGTGAGCGGCGCGCCACACGATGATGACTCa                                                                                 >  2:5159073/1‑150 (MQ=14)
                                                                     tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                             >  2:3207194/1‑150 (MQ=255)
                                                                       gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                           >  2:4366983/1‑150 (MQ=255)
                                                                        tgggCACCGGGCTGCGGCTTTTTTGTTCGCCTGTGCCGGGGGGTGCTGGGGGCGGGGGGGTTTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAtata                                                                          <  2:3655166/149‑1 (MQ=21)
                                                                             aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                                     >  2:819078/1‑150 (MQ=255)
                                                                             aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACCAATaacaac                                                                     >  2:3989136/1‑150 (MQ=255)
                                                                                       ggtttttttGTTCGCCTGTGCGGGGGGGTGCTGGGGGGGGGGGTTTTGGG‑CGGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                           <  2:328574/145‑1 (MQ=25)
                                                                                       ggcgTTTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGCGATGGG‑CCGCCCCGTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAATATATAACAACATCTTaaaaa                                                           <  2:6566313/150‑1 (MQ=25)
                                                                                       ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                           >  2:2060565/1‑150 (MQ=255)
                                                                                        gtttttttGGGCGCCCGCGCCGGGGGGGGCGGGGGGGGGGGGCGCGGGC‑CCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                          <  1:1001273/146‑1 (MQ=255)
                                                                                        gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                          >  2:101465/1‑150 (MQ=255)
                                                                                         cttttttGTCCGCCTGCGCCGGGGGGTGGTGGGGGCGGGGGCTTTGGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                                         <  1:3978488/147‑1 (MQ=255)
                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                        >  2:665801/1‑150 (MQ=255)
                                                                                            ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACATATAACAACATCTTAAAAAAAGGc                                                      >  2:4863665/1‑150 (MQ=255)
                                                                                             ttggtgcGCCTGAGCCGGGGGGGGCTGGGGGGGGCGGCTTTGGG‑CCGCCCCCTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                                     <  2:4030519/144‑1 (MQ=255)
                                                                                                 ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGAc                                                 >  2:1373461/1‑150 (MQ=255)
                                                                                                 ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATAAGTTGCTTAGCTCAAGGAAGGAGCGCAACAAGATGATGTAACAACCTATAACAACAGCTTAAAAAAAGGCCTGAc                                                 >  2:5547708/1‑150 (MQ=25)
                                                                                                           cgggggggTGCTGGGGG‑GGGGGGGTTTGGCCCCCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                                       <  1:4838674/147‑1 (MQ=38)
                                                                                                            ccGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                      >  2:5458017/1‑150 (MQ=255)
                                                                                                              gggggTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAGGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                                    >  2:2803836/1‑150 (MQ=255)
                                                                                                                 cgcgCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                                 >  2:56297/4‑150 (MQ=255)
                                                                                                                   tGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                               >  1:739444/1‑150 (MQ=255)
                                                                                                                    gcggggggcgggggcGATGGC‑CCCCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                              <  1:4628776/147‑1 (MQ=38)
                                                                                                                    gCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                              >  2:2238436/1‑150 (MQ=255)
                                                                                                                      ggggggCGGCGGCTATGGG‑CCCCCCCGGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCACCAAGCTGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCttttttta                            <  2:4102299/149‑7 (MQ=255)
                                                                                                                         gggCGGGGGCGTGGGG‑CCGCCCCCGC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                         <  2:1316828/150‑1 (MQ=255)
                                                                                                                         gggCGGCGGCGATTGG‑TCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                         >  2:5376453/1‑150 (MQ=255)
                                                                                                                         gggCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCACGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                         >  1:3809809/1‑150 (MQ=255)
                                                                                                                           ggggcggccAATCG‑CCGCCGCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt                       >  1:4339917/3‑150 (MQ=255)
                                                                                                                           ggggcggccAATCG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTCAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt                       >  2:758353/3‑150 (MQ=255)
                                                                                                                           ggggcggccAATCG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt                       >  2:1450085/3‑150 (MQ=255)
                                                                                                                           ggggcggccAATCG‑CCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt                       >  1:3536655/3‑150 (MQ=39)
                                                                                                                           gcggcggcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt                       >  1:1386780/1‑150 (MQ=255)
                                                                                                                              ggggcgTTGGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    <  2:5288130/148‑1 (MQ=255)
                                                                                                                              gcggcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    >  2:1267472/1‑150 (MQ=255)
                                                                                                                              gcggcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTCACATTACGCCAGGCCTTCTGCGAtaattaat                    >  2:3920365/1‑150 (MQ=255)
                                                                                                                              gcggcgAGTGG‑CCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    <  2:4026823/150‑1 (MQ=255)
                                                                                                                                gggctttGG‑CCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc                   <  2:1132445/145‑1 (MQ=255)
                                                                                                                                 gcgATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc                 >  1:5624396/1‑150 (MQ=255)
                                                                                                                                   gATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCCGGCCTTCTGCGTTAATTAATCACGc               >  2:6090519/1‑150 (MQ=255)
                                                                                                                                   gATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               >  2:1832263/1‑150 (MQ=255)
                                                                                                                                   gATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               >  2:4315203/1‑150 (MQ=255)
                                                                                                                                   cttgggcccgcccgcTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAATAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                <  2:5117559/143‑1 (MQ=255)
                                                                                                                                   cgtgggcccgccCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                <  2:6116672/143‑1 (MQ=255)
                                                                                                                                    tgggg‑CCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt              <  1:4533863/147‑1 (MQ=255)
                                                                                                                                        g‑ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          >  1:5813138/1‑150 (MQ=255)
                                                                                                                                        g‑ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          >  1:1271641/1‑150 (MQ=255)
                                                                                                                                         cccgccgCCCC‑CGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          <  2:2335776/149‑1 (MQ=255)
                                                                                                                                          ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         >  1:2076122/1‑150 (MQ=255)
                                                                                                                                          ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         >  1:1190050/1‑150 (MQ=255)
                                                                                                                                          ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         >  2:1862704/1‑150 (MQ=255)
                                                                                                                                           cgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca        >  2:3716671/1‑150 (MQ=255)
                                                                                                                                                ccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg   >  2:5008052/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  NC_000913/915081‑915374

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: ‑