Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 11.3% 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,2260TG11.3% 546.8 / 5.9 287intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (163/46);  new base G (7/26);  total (202/85)
Fisher's exact test for biased strand distribution p-value = 3.94e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  NC_000913/915077‑915374
                                                                                                                                                     |                                                                                                                                                     
caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                                       >  2:1461326/1‑150 (MQ=32)
 aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                                      >  1:593418/1‑150 (MQ=32)
 aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                                      >  2:5464859/1‑150 (MQ=32)
       gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCt                                                                                                                                               >  1:2381504/1‑150 (MQ=32)
           gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCtt                                                                                                                                           >  2:1150410/1‑150 (MQ=32)
           gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCtt                                                                                                                                           >  2:2323877/1‑150 (MQ=32)
               aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTgcgc                                                                                                                                       >  1:314283/1‑150 (MQ=17)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                                    >  1:2084196/1‑150 (MQ=32)
                   cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCCGt                                                                                                                                   >  2:6964494/1‑150 (MQ=17)
                    gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCATGCGCCAGTg                                                                                                                                  >  1:6986922/1‑150 (MQ=17)
                      gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCATGCGCCAGTGGc                                                                                                                                >  2:6624387/1‑150 (MQ=18)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAAAATGt                                                                                                                         >  2:176194/1‑150 (MQ=255)
                              gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                                        >  2:1846878/1‑150 (MQ=40)
                              gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTG‑CGCGCTCCGTGCGCCAGGGGCCATATGtt                                                                                                                        >  2:3964788/1‑150 (MQ=255)
                                  tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAG‑CGCGCGCCATGCGCCAGTGGCAATATGTTTCtt                                                                                                                    >  2:3840598/1‑150 (MQ=14)
                                   gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                                   >  1:3161475/1‑150 (MQ=255)
                                      aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                                >  1:373308/1‑150 (MQ=255)
                                      aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                                >  1:2531314/1‑150 (MQ=255)
                                       ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                                               >  2:8578853/1‑150 (MQ=255)
                                        cgcgcccccgccccGGTCGCCTTAGGGTATAAGTGGGGGCCACCGCGCCGGGGCTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCTTTTGGCCGCCCCCTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                              <  2:423266/148‑1 (MQ=12)
                                        cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                              >  2:2559788/1‑150 (MQ=255)
                                         tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAG‑CGCGCGCCTTGCGCCAGGGGCAATATGTTGCTTAGCTCAt                                                                                                             >  1:2836073/1‑150 (MQ=255)
                                          cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                            >  2:3282813/1‑150 (MQ=255)
                                                 gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACAGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                                     >  2:7257035/1‑150 (MQ=255)
                                                     cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGca                                                                                                 >  2:6005176/1‑150 (MQ=255)
                                                       gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                               >  2:3008289/1‑150 (MQ=255)
                                                            cGTCAGGGATCTGTCGCGGGCACCGCGCTGCGGCGATTTTGTTCGCCTGAGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                          <  2:6520223/150‑1 (MQ=18)
                                                            cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCATAGCTCATGAAAGGAGCGCAACAAgat                                                                                          >  2:5895904/1‑150 (MQ=255)
                                                             gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                                         >  2:1024479/1‑150 (MQ=255)
                                                               atggtatatGTCGCGGGCACCGCGCCGGGGCTTTTTGGGTCGCCTGCGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                                       <  2:2859952/145‑1 (MQ=18)
                                                                  ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                                                                    >  2:676871/1‑150 (MQ=255)
                                                                   tatatGTCGCGGGCCCCGCGCCGCGGCGATTTGGTTCGCCAGCGCGGGGGGGTGCTGGGGGGGGGGGCTTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                                   <  2:4313730/149‑1 (MQ=18)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATTTTTCTTAGCTCATGAAAGGAGCGCAACCAGATGATGTAt                                                                                   >  1:7779751/1‑150 (MQ=18)
                                                                     tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCa                                                                                 >  2:1375487/1‑150 (MQ=255)
                                                                      aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAG‑CGCGCTCCATGCGCCCCGGGCAAGATTTTTCTTAGCTCAAGGAAGGGGCGCAACCAGATTATTaaaaaa                                                                                >  1:8486287/1‑146 (MQ=14)
                                                                        gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGGCGCCCCCAG‑CGCGCACCATGCGCCAGTGGGCATAATTTGCTTAGCGCAAGAAAGGTGCGCAACAAGATTATGAAacaaca                                                                              >  1:1161049/1‑150 (MQ=17)
                                                                        gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTTGCCGCCCCCCG‑CGCGCTCCATGCGCCCGTGGCAATATGTTTCTTAGCGCATTAAAGGAGCGCAACAAGAATTTTAATCAACa                                                                              >  1:5097073/1‑150 (MQ=17)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCCCCCCCCG‑CGCCCGCCGTGCGCCCCGGGCAATATGTTGCGTATCGCATGGAAGGAGCGCAAAAAGATGATGAATCAACaa                                                                             >  1:6079356/1‑149 (MQ=17)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                             >  2:7961341/1‑150 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                             >  2:7940702/1‑150 (MQ=255)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                           >  2:6441945/1‑150 (MQ=255)
                                                                             gggCACCGGGGTGGGGCTTTTTTTTTCGCCAGAGCCGGGGGGTGCTGGGGGGGGGGGGGTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                         <  2:503405/150‑1 (MQ=21)
                                                                             gggAACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTGGTGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                         >  2:8104463/1‑150 (MQ=255)
                                                                             ggccACCGCGCCCCGGCGATTTTGTTCGCCTGAGCCCGGGGTTGCAGGGGGCGCCGGCTATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                         <  2:7146639/147‑1 (MQ=21)
                                                                              ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaac                                                                        >  2:7402048/1‑150 (MQ=255)
                                                                                cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                                      >  2:2828358/1‑150 (MQ=255)
                                                                                cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                                      >  2:752534/1‑150 (MQ=255)
                                                                                         gcggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                             >  2:5400443/1‑150 (MQ=255)
                                                                                          cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                            >  2:5459531/1‑150 (MQ=255)
                                                                                          cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                            >  2:3533762/1‑150 (MQ=255)
                                                                                            gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCCCTCCTTGCGCCAGTGGCAATATGTTGCATAGCTCAAGAAACGAGCGCAACAAGAAGATGAATCAACAGATAACAACAACATaaaaaa                                                          >  2:5911499/1‑150 (MQ=21)
                                                                                              gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                        >  2:3464174/1‑150 (MQ=255)
                                                                                               aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                                       >  2:4541252/1‑150 (MQ=255)
                                                                                                        gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGCGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACCtt                                              >  2:7859951/1‑150 (MQ=255)
                                                                                                          ctgcgCCGGGGGGGGCTGGGGGGGGCGGGTTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                            <  2:2023773/146‑1 (MQ=37)
                                                                                                          ctgcgCCGGGGGGGGCGGGGGGGGGCGGCTTTTGGCCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                            <  2:4646699/146‑1 (MQ=255)
                                                                                                           agagCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                           >  2:8226895/1‑150 (MQ=255)
                                                                                                               cccGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTG‑CGCGCGCCTGGCGCCAGGGGCCATATGTTGCGTAGCTCCGGCAAGGCGCGCCACAAGATGATGAATCAACCGAGACCACCCTCTTAAAAAAAGGCCTGACCATACGcccg                                       >  2:6243304/1‑148 (MQ=16)
                                                                                                                ccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                      >  2:7722282/1‑150 (MQ=255)
                                                                                                                 cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                                     >  2:2152629/1‑150 (MQ=255)
                                                                                                                  gggggTGCA‑GGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCt                                   >  2:5808976/1‑150 (MQ=255)
                                                                                                                      gtgggggggggggggggGATGGCCCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                                <  2:1076393/145‑1 (MQ=255)
                                                                                                                      gTGCGGGGGGGGGGGGCGTTGGCCCCCCCCCTC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                                <  1:567984/150‑1 (MQ=255)
                                                                                                                      gTGCGGGGGGCGGCGGCG‑TTGGCCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttt                                <  2:312718/150‑3 (MQ=255)
                                                                                                                       tGCTGGGGTGGGGGGCGATGGGACGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                               <  1:176194/150‑1 (MQ=38)
                                                                                                                        gctgggggCGGCGGCGTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                              <  1:5108770/147‑1 (MQ=255)
                                                                                                                        gccgggggcgggggcCTGGGGCCGCCCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                              <  1:7949083/147‑1 (MQ=38)
                                                                                                                          aGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGCAAGGAGCGCAACAAGATGATGCAACAACATATAACAACATCTGAAAAAAAGGCCAGACATAACGACAGCCCTACTGCGt                            >  2:4784505/1‑150 (MQ=37)
                                                                                                                          aGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGACAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCgg                            >  1:5295411/1‑149 (MQ=255)
                                                                                                                            ggggCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                          >  2:6307210/1‑150 (MQ=255)
                                                                                                                            ggggCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAACTTAAAAAAATGCCTGACATTACGCCAGGCCTTCTGCGtta                          >  1:759606/1‑150 (MQ=255)
                                                                                                                            ggggCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                          >  1:6591368/1‑150 (MQ=255)
                                                                                                                             gggCGGCGGCGTTGGGCCCCCCCCGG‑CGCGCTCCTAGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                         <  1:266055/150‑1 (MQ=255)
                                                                                                                                  gcggcgTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    <  2:2723900/150‑1 (MQ=255)
                                                                                                                                  gcggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    >  2:5606861/1‑150 (MQ=255)
                                                                                                                                  gcggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                    <  1:5606861/150‑1 (MQ=255)
                                                                                                                                    ggcgTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGACCGCAACAAGATGATGAATCAACATATAACAACATCTTATAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  <  2:8469299/150‑1 (MQ=255)
                                                                                                                                    gccgaggggCCCCCCCCGC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                  <  1:397190/143‑1 (MQ=255)
                                                                                                                                     gcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc                 >  1:5555497/1‑150 (MQ=255)
                                                                                                                                       gtttGCCCCCCCCCCG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc               <  1:2950845/148‑1 (MQ=255)
                                                                                                                                        ctgcgcccccCCCGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt              <  2:7363419/146‑1 (MQ=255)
                                                                                                                                         tggccccccccccc‑cGCGCTCCTTGCGCCAGTGTCAATAAGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt             <  1:5827787/143‑1 (MQ=255)
                                                                                                                                           ggCCGCCCCGTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCACAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtttt           <  2:4924857/150‑1 (MQ=255)
                                                                                                                                           ggCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCACAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtttt           <  2:5140874/150‑1 (MQ=255)
                                                                                                                                            gccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          >  1:4348835/1‑150 (MQ=255)
                                                                                                                                            gccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          >  1:4830959/1‑150 (MQ=255)
                                                                                                                                            gccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          >  1:8337505/1‑150 (MQ=255)
                                                                                                                                            cccccccccGG‑CGCGCTCCTTGCGCCAGTGGCACTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAGGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc          <  1:6781106/146‑1 (MQ=255)
                                                                                                                                             ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAATAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTATTCACGCTTTTcc         >  2:3694672/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc         >  2:7687012/1‑150 (MQ=255)
                                                                                                                                              cgccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca        >  2:8361190/1‑150 (MQ=255)
                                                                                                                                              cccccccGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca        <  1:6150338/148‑1 (MQ=255)
                                                                                                                                               gccCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGACTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag       >  2:7770948/1‑150 (MQ=255)
                                                                                                                                                cccccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagc      >  1:1570686/1‑150 (MQ=255)
                                                                                                                                                cccccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagc      >  1:4113744/1‑150 (MQ=255)
                                                                                                                                                 ccccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca     <  1:7845011/150‑1 (MQ=255)
                                                                                                                                                 ccccGG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca     <  1:547880/150‑1 (MQ=255)
                                                                                                                                                  cccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag    >  1:6323636/1‑150 (MQ=255)
                                                                                                                                                  cccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag    >  2:6167196/1‑150 (MQ=255)
                                                                                                                                                   ccTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg   >  2:3034367/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGG  >  NC_000913/915077‑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: ‑