Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,5910A.42.5% 147.0 / 122.3 73intergenic (+18/‑16)metW/metVtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (38/4);  new base . (29/2);  total (67/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

AAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947478‑2947740
                                                                                                                    |                                                                                                                                                     
aatacgccaattACGCTCACGCGGGCGGGAGCAGCCTGGTCGCTCTGCGGCATCATAACCCGAAGGGCGTCGGGTCAAAACCGCCCCCCGCAACCAA‑TA‑AAA‑TATGCTTCACTACAAGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                                  <  8:473033/140‑1 (MQ=2)
                                              tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTCATCAAGTAAAAGCAGTACGGTGAACCGGGCTCGAGAAGCCTGGAAGCTCGACGGGCTCATAACCCCAACCTCGTCCGGTCACATCCGcccc                                                                    >  3:738715/1‑151 (MQ=14)
                                                    tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGc                                                                                                                                                  <  1:1132609/67‑1 (MQ=255)
                                                    tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGc                                                                                                                                                  >  5:1130921/1‑67 (MQ=255)
                                                             gAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TG‑AAA‑TTTTATGAAGTAAAAg                                                                                                                                                   >  4:218054/1‑57 (MQ=255)
                                                             gAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAT‑GTACATGCTGTAACAGCAGTAGTTGGAAGGGGGCTCGAGAAGAATGGTAGCTAGGAGGGCTGATTAACCGTACGGGGACGGTACCACTCCGCCCCCCGCCACaaaagatt                                                     >  2:57121/1‑146 (MQ=2)
                                                             gAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAg                                                                                                                                                   <  2:217787/57‑1 (MQ=255)
                                                                              aTCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCCGCCCCCGCAACAAATTATTGAACACCCAAAAGGgtg                                    >  4:1081329/1‑151 (MQ=255)
                                                                                   gCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTACCGGGTGttttt                               >  1:1051831/1‑151 (MQ=255)
                                                                                   gCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTCGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttt                               >  7:885/1‑151 (MQ=255)
                                                                                     ccccGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  5:1408643/1‑151 (MQ=255)
                                                                                     ccccGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  5:939928/1‑151 (MQ=255)
                                                                                     ccccGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  6:1313707/1‑151 (MQ=255)
                                                                                     ccccGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  2:189050/1‑151 (MQ=255)
                                                                                       ccGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGt                           >  6:1372138/1‑151 (MQ=255)
                                                                                        cGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGtt                          >  6:916093/1‑151 (MQ=255)
                                                                                         gCAACCAATTA‑AAA‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGCAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATAATATAACACCATCACGTGTCTGTtttttttt                         >  6:1279071/1‑151 (MQ=11)
                                                                                         gCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTGATTTGttt                         >  8:1345369/1‑151 (MQ=255)
                                                                                          cAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc                        >  8:1287025/1‑151 (MQ=255)
                                                                                           aaCCAA‑TTAAAA‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCTAAAGAGGGttttttttttacg                       >  3:531855/1‑148 (MQ=11)
                                                                                           aaCCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  1:541067/1‑151 (MQ=255)
                                                                                            cccAA‑TC‑AAA‑TGTTGTGCAGGAAAAGCAGTTCGTTGTCGGGGCGGTGACCAGCCTTGTAGCTCGTCGCGGTCATAACCGGACTGGAGTCGGTTCAAATCCGGCGCGCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTg                      <  2:580895/150‑1 (MQ=255)
                                                                                            aCCAA‑TT‑AAATTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtt                                  <  8:1037316/139‑1 (MQ=255)
                                                                                            aCCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGGCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTg                      >  2:344204/1‑151 (MQ=255)
                                                                                             ccAA‑TC‑AAA‑TTTGATGAAGTAAAAAAAGTACCGTGACGCGGGGAGGAGAAGCCTGGTACCTCGTCGGGCTCATAAACATAAGGTCGTAGTTTCAAATCAGGCCCACGAAACCAATTATTGAACACACTAACTGGTGTTATTTTGTTTCAgg                     >  3:1002687/1‑151 (MQ=255)
                                                                                             accaattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTAACGGGTGtt                                  >  3:1041639/9‑139 (MQ=11)
                                                                                             accaattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACAACCTAAATAGTGTTTTTTTTTTTCTg                      >  4:702035/9‑151 (MQ=11)
                                                                                              ccaattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTATCAACCTGAACTGTGTTTTTTTGTATCTgg                     >  1:1227865/8‑151 (MQ=1)
                                                                                               caattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCTAACGGCGCtttttttttaaggggt                    >  1:609904/7‑142 (MQ=11)
                                                                                               caattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCTCACGGGGGTGTTGTTGTTTCTggg                    >  8:787610/7‑150 (MQ=11)
                                                                                               caattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTATCAGCCTAACGGCTGTGTTGTTGGTTCgggg                    >  2:355969/7‑150 (MQ=11)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGCGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  4:600499/1‑151 (MQ=255)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTGAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtg                   >  4:848647/1‑150 (MQ=255)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtg                   >  8:1272554/1‑150 (MQ=255)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  2:87505/1‑151 (MQ=255)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  8:508681/1‑151 (MQ=255)
                                                                                               aa‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  4:936549/1‑150 (MQ=255)
                                                                                                aattaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCTAAAGGGTGTTTTTTTGTTTCTGGtg                   >  7:1456370/6‑150 (MQ=14)
                                                                                                a‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtct                  >  5:1343698/1‑151 (MQ=255)
                                                                                                 attaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTTACGTGGGTTTTTTTTGTTCTGGgct                  >  7:442737/5‑151 (MQ=11)
                                                                                                 attaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACCTAAAGGGTGTTTTTTTGTTTCTGGtct                  >  3:337638/5‑151 (MQ=14)
                                                                                                 attaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTTAACACGATAAAGCGGGTTTTTTTGGTTCGGGgct                  >  6:767449/5‑151 (MQ=1)
                                                                                                 attaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCTAAAGGGGCTGTATTTGTTTCTGGtct                  >  5:23912/5‑151 (MQ=14)
                                                                                                  ttaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTAAAGGGTGTTTTTTTGTTTCTGGtctc                 >  1:858186/4‑151 (MQ=255)
                                                                                                  ttaaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCAAACGGGTGTGTTTTTGTTTCTGGtctc                 >  2:881089/4‑151 (MQ=14)
                                                                                                  tC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc                 >  6:1123205/1‑151 (MQ=255)
                                                                                                  tC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc                 >  8:1393522/1‑151 (MQ=255)
                                                                                                   taaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCAGCCTATCGGGTGTTTTTTTGTTTGTGGTCTcc                >  2:511079/3‑151 (MQ=14)
                                                                                                    aaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTAATGAACACCCTGACGGGTGTTTTTgtgtgggggggctcgc               >  5:666836/2‑139 (MQ=14)
                                                                                                    aaaa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCAAACGGGGCTTTTTTTGTATCAGGTCTcc                >  7:1287150/2‑150 (MQ=14)
                                                                                                      aa‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  8:637426/1‑150 (MQ=255)
                                                                                                      aa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTAAAGTGTGTTTTTTTGTATGGGGTCTCCCa              >  2:658395/1‑150 (MQ=14)
                                                                                                      aa‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGCTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCTAACCGGTGTTTTGTTGTTTCTGGTCTCCCAt             >  6:440963/1‑151 (MQ=17)
                                                                                                       a‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  7:817855/1‑151 (MQ=255)
                                                                                                       a‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTTTTTGTTTCTGGTCTCCCa              >  7:1331868/1‑150 (MQ=255)
                                                                                                       a‑TTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCATAACGGGGGTTTTTTTGTTTCGGGTCTCCCaa             >  3:514754/1‑149 (MQ=17)
                                                                                                       a‑TTTGATGAAGT‑AAAGCAGTACCGTGACCCGGGGTGGAGCAGCCTGGTAGCACGTCGCCCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACCAAAACAGTCAATTTTTGTTTCAGCTCTCCCaa             >  3:712739/1‑149 (MQ=11)
                                                                                                         tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCAAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  7:863586/1‑150 (MQ=255)
                                                                                                         tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACAACCTAAAGGGTGTTTTTTTGTTTCTGGTCTCCCaaa            >  3:334984/1‑148 (MQ=17)
                                                                                                          ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  6:1299167/1‑150 (MQ=255)
                                                                                                           tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  2:478005/1‑150 (MQ=255)
                                                                                                           tgatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTCACGGGTGTTTTGTTGTTTCTGGTCTCCCATaaa          >  8:23800/1‑150 (MQ=255)
                                                                                                            gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaa          >  6:1097906/1‑150 (MQ=255)
                                                                                                             atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAGAACACGCTATAGCGTCTTTATTTGTTACGGGGCTCCCCTaaaaac       >  5:897234/1‑150 (MQ=11)
                                                                                                              tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  6:267636/1‑150 (MQ=255)
                                                                                                              tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACATAAACCGTCTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  1:245988/1‑149 (MQ=17)
                                                                                                              tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATAGAACACCCTCACGGCGGTTTTTTTGTTTCTGGGCTCCCATaaaaacg      >  6:1418056/1‑149 (MQ=17)
                                                                                                               gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  6:309760/1‑149 (MQ=255)
                                                                                                               gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   <  2:297766/137‑1 (MQ=255)
                                                                                                               gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGATCACCCTAACAAGTGTTTTTTTGtaacaggtc                   >  4:298152/1‑129 (MQ=14)
                                                                                                                gacgaaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAg      <  5:857429/146‑1 (MQ=255)
                                                                                                                aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  7:569570/1‑151 (MQ=255)
                                                                                                                aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:625562/1‑151 (MQ=255)
                                                                                                                 aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCCGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  8:387455/1‑151 (MQ=255)
                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  2:1142929/1‑150 (MQ=255)
                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  3:551025/1‑150 (MQ=255)
                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  4:693377/1‑150 (MQ=255)
                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  4:1424478/1‑150 (MQ=255)
                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTGTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  7:541868/1‑150 (MQ=255)
                                                                                                                    |                                                                                                                                                     
AAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAA‑TTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947478‑2947740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑