Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,670 A→T 16.0% V64E (GTA→GAA)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6700AT16.0% 114.7 / 3.6 65V64E (GTA→GAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (11/40);  new base T (8/3);  total (19/46)
Fisher's exact test for biased strand distribution p-value = 2.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTC  >  NC_000913/361439‑361889
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑Tccccaga                                                                                                                                                                                                                              >  1:118795/1‑235 (MQ=255)
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGCTGGGCGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAAACGAAAAAGAAACAAACAGGCGGTTTT‑ATCCacac                                                                                                                                                                                                                              >  1:379530/1‑236 (MQ=255)
          tgTTAATGCCAATCGTTATCGGAAAAGAGTTCATCTCGGCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTccc                                                                                                                                                                                                                     <  1:975179/236‑1 (MQ=255)
              aTTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATTGATGTTGGAACCGTAAGAGAAATAGAAAGGCGGATCTCACCAAGGG‑GTTTTtccacacc                                                                                                                                                                                                                >  2:894396/1‑234 (MQ=255)
                       tgTTATCGGAAAAGGGTACATCTCGCCGTTTTTGCTCAATTCATGGTGTACAGTGTGTCCCGTAACGGAAAGAGTAAAGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACa                                                                                                                                                                                                   <  2:853359/238‑1 (MQ=255)
                                       tACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGCCGGGGCT‑TCCCCGCG‑GTGTTTCCCTC‑CCCTGGGA‑AAACATTT‑TTTAATAAggta                                                                                                                                                                                   >  2:691349/1‑237 (MQ=255)
                                                                     tgtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATTAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGa                                                                                                                                                   <  1:507610/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCATTACGTTATCGCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGAt                                                                                                                                              <  1:919086/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAAAGGAAAGGAGAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGGTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGAt                                                                                                                                              <  1:408838/239‑1 (MQ=255)
                                                                            ggTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCCCGC‑GGTTTCCCCCC‑CCGC‑GGC‑AAAAAATTCTTTTATCCAGCATTTTCTTTTTTTAACTTCTG‑TTGG‑ATGGGCGTGGTTa                                                                                                                                            >  2:931228/1‑239 (MQ=255)
                                                                               gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTA‑‑AGAAATAGACAGGCGGGTCT‑ACCCAGGCGGTTTTCCCCT‑‑CCGT‑GGCAAAACATTTTTTTAAT‑AGGATTTCTCTATTTTCAACTTTCG‑TTGG‑ATTGGAGttgatatact                                                                                                                                        >  1:922358/1‑229 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTTTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCCA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTc                                                                                                                                       <  1:847511/239‑1 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTc                                                                                                                                       <  2:894630/239‑1 (MQ=255)
                                                                                               agTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCCAGC‑GGTTTTCCCCA‑CCCC‑GGC‑AAACATTTCTTTTTTTAGGCTTTTTCTTTTTTTTCATTCTC‑TTGG‑ATGAGAGTGATTAAACTCCTAAA‑TTaaatttt                                                                                                                        >  2:746293/1‑234 (MQ=255)
                                                                                                 tAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATATACAGGCGGGTTT‑ACCCCCGC‑GGTTTCCCCTT‑CCCTGGGC‑AAAAATTT‑TTTTATTAGGCGTTTTCTTTTTTTTACTTCTG‑TTGG‑TTTGTAGTGTTTTAAATCATCAA‑ATAAA‑Gtttt                                                                                                                      >  2:782579/1‑238 (MQ=255)
                                                                                                       ttggTGCAATCTGTTCTTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGa                                                                                                               <  2:385819/237‑1 (MQ=255)
                                                                                                         ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GGTTTCCCCTT‑CCGT‑GGC‑AAACATATCTTTAATCAGGCTTTCTCTTTTTTTAAATTCTG‑GTGG‑CTGCGAGGTGTTAAACACATACC‑TTTAA‑GTTTTTCCCCgcaa                                                                                                             >  1:842942/1‑236 (MQ=255)
                                                                                                                tCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTT‑AACCAAGC‑GTTTTTCCCCT‑CCGG‑GGC‑AAACATTTTTTTAATTAGGCTTTTCTTTTTTTTAACTTCTG‑ATGGTATGCGGGGGTTTAAACCCATA‑A‑TTAAC‑CATTTTCCCCAAAGGtttttt                                                                                                      >  1:299954/1‑239 (MQ=255)
                                                                                                                     aCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTT‑AACCAGGCGGTTTTCCCCAA‑CCG‑‑GGG‑AAACAATTCTTTTATCAGGCTTTCTTTTTTTTTACCTTCTT‑ATTG‑GTGCGCGGCGTTAAAATCAACAA‑TTAAA‑GTTTTTCCACCAAGACTTTTTTTCgg                                                                                                 >  1:550415/1‑239 (MQ=255)
                                                                                                                       gTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTa                                                                                               <  1:88016/239‑1 (MQ=255)
                                                                                                                         tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCCAGC‑GTTTTTCCCCC‑CCCT‑GGC‑CAAAATTTCTTTATTACAGCTTTTTCTTTTTTTTACTTCTT‑TTGG‑ATTGGCGTGGTTTAAATCAAACA‑TTTAA‑TTTTTTCCCCCAAATCttttttttcttaa                                                                                              >  2:88016/1‑231 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TTCCAAGC‑GTTTTTCCCCA‑CCGC‑GGC‑AAACATTTCTTTTATTCAGCTTTCCTTTTTTTTCACTTTTT‑ATTG‑ATTCGCGTGATTAAAATCCTACATTTAAC‑G‑TTTCCCCCCCAGTCCTTTTTTCGGTAAACCtt                                                                                         >  1:358087/1‑238 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCCGGC‑GTGTTTCCCTA‑CCGT‑GGC‑AAAAATTTCTTTTATCACGCTTTCTCTTTTTTTCACTTCTT‑CTGT‑AGTCGTGTGTTTAAAATCATACACTTAAC‑‑TTTTTCCACCAACGCTTTTTTTCGGTTAACCTTc                                                                                        >  2:666559/1‑239 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATAGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGATCT‑ACCCTCGCGGTTTCCCCCTC‑CCGC‑GCA‑AAACTTTATTTTTAT‑AGGCCTTTTCTTTTTTTCACTATTG‑ATGG‑ACGCAGGGGGTTATACTCAAACA‑TTTACTTTTTTTCCCCACAGTCTTTTTTcgttaaacatt                                                                                          >  2:504585/1‑228 (MQ=255)
                                                                                                                                 atggTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCCGGA‑GGTTTCCCCCACCCGT‑GCC‑AAA‑AATTCTTTTATCAGGCTTTCTCTTTTTTTCACTTCTGAATGG‑ATGAG‑GGGATTTAAATCCTACATTTACC‑GTTTTCCCCCCAAGTCTTTTTTCGGGAAACCTTTcc                                                                                       >  1:579098/4‑237 (MQ=255)
                                                                                                                                 aTTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAAGC‑GTTTTTCCCTA‑CCCG‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTTAACTTCTG‑TTTG‑TGGCCACGGATTAAACTCCTAAG‑TTAAC‑ATTTTTCCCCCAAGGCTTTTTTTCCGTTAGCCTTCgac                                                                                     >  1:428813/1‑237 (MQ=255)
                                                                                                                                   tttttccGTGTAGTCAATGACAATTGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTGCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  2:505559/234‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:746293/239‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:329366/239‑1 (MQ=255)
                                                                                                                                      tACCGTGGAGTCATCGACAATGGTTAAATTGAAATTGGCATAAAAATGGCGGCCTATATGGATGGTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACatat                                                                                <  1:500311/239‑1 (MQ=255)
                                                                                                                                        ccGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:782575/234‑1 (MQ=255)
                                                                                                                                                      aaaaTGGTTAAATTGAAATTTGCATAAAAATTGGGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACCCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  2:413753/237‑1 (MQ=255)
                                                                                                                                                       gattGGTTAAATTGAAATTTGCAAAAAAATTGCGGCGTATAAGGATGTTGGAACCGTAAGGGAAATAGTCAGGCGGGTGT‑ACCCTGGC‑GTTTGCGTCTT‑CCTG‑GGC‑AAAAATTTCGTTAATCAGGCTTTCTCTTTTTTCAACGTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTACCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCTCATATCGGTAAATAGCTTGcc                                                                <  2:507874/235‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACtt                                                                                                                                                                  <  1:210597/144‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:894396/238‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTGGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:1083305/238‑1 (MQ=255)
                                                                                                                                                           tgTTAAATTGAAATTTGCATAAAAATTGCGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGctc                                                           <  1:596445/238‑1 (MQ=255)
                                                                                                                                                                    gAAATTTGCATAAAAATTGCGGCCTATAGCGATGTTGGAACCGTAAGAGAAAAAGTATGGCGGTTCT‑ACCCATGC‑GTTTTCCCCTA‑CCGT‑TGC‑AAACATGTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCttt                                                  <  1:504585/239‑1 (MQ=255)
                                                                                                                                                                        attGAATAAAAATTTGGGACTATAGGGATGTTGGAACCGTA‑‑AGAAATGTACAGGCGGTTCT‑CCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:922358/238‑1 (MQ=255)
                                                                                                                                                                          ttcaaagaaaTTGGGGCCTATATGGATGTTTGACCCGTAAGAGAAATAGACAGGCGGTGCT‑ACCCAGGC‑GTTTTCCCCTT‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGGCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:118795/232‑1 (MQ=255)
                                                                                                                                                                          tgactAAAAATTGCGGCCTATTTGGTTGTTGGAACCGTAAGAGAAATAGACCGGCGGTTCT‑CCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:75292/235‑1 (MQ=255)
                                                                                                                                                                          tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  1:533688/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGGGGGCTTTATTGATGTTGGACCCCGAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  2:569451/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  1:931228/239‑1 (MQ=255)
                                                                                                                                                                                  aaTTGCGGCCTATTTGGATGTTTGAACCGTAAGAGAAATAGACAGGCTGTTTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATg                                    <  2:229719/239‑1 (MQ=255)
                                                                                                                                                                                         gCCTATATGGATGTTGGAACCGTAAGAGAAAAAGTCAGGCGGTTTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAt                              <  2:381648/238‑1 (MQ=255)
                                                                                                                                                                                             atgtGGATTTTGGAACCGTAAGGGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:560783/236‑1 (MQ=255)
                                                                                                                                                                                             atatGGATGTTTGAACCGTAAGAGAAATAGACAGGCGGTTCT‑CACCAGGC‑GTGTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTGGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:409652/239‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑cc                                                                                                                                                                                                               <  2:415234/62‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTc                                                 <  1:861215/214‑1 (MQ=255)
                                                                                                                                                                                               tggggtttttgAACCGTAAGAGAAATAGACAGGCGGTTCT‑TCCCTGGC‑TTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCCGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  2:985772/229‑1 (MQ=255)
                                                                                                                                                                                               atGGATGTTGGAACCGTAAGAGAAATAGACCGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTT‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  1:813402/239‑1 (MQ=255)
                                                                                                                                                                                                 ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                     <  1:889057/239‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCTTAAGAGAAATAGTAAGGCGGTTCT‑ACCCAGGC‑GTTTTCACCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACGTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCTGTAAGCCTTCGCTCATATCGGTAAATAGCTTGTCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTcc                   <  2:470274/236‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:579098/237‑1 (MQ=255)
                                                                                                                                                                                                     ttttgAACCGTTAGGGAAATATACCGGTGGTTTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:474976/234‑1 (MQ=255)
                                                                                                                                                                                                     gTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  1:674025/238‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGGAAAATAGACAGGCGGTTCC‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTa               <  2:435365/239‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGtt                <  1:933058/238‑1 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGTTTT‑TCCCAGGC‑GGTTTTCCCTA‑CCCT‑GGC‑caaaa                                                                                                                                                                                                   >  1:415234/1‑60 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGTTTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGG‑GGC‑AAACATTTCTTTTATCAGGCTTTCTCTTTTTTTAACTTCTG‑ATGG‑ATGCGAGGGATTAAACTCATACA‑TTAAC‑GTTTTTCCCCGAAGTACTTTTTTCGGTAAGCCTTTCCACATATCGGTAAATAGCTTTCCTGCTCTTTTTTTTTCTGTCCTTTGCATTTTCCATGCGAGCACTCCGtttt              >  1:591705/1‑237 (MQ=255)
                                                                                                                                                                                                            ccGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGata          <  1:64759/239‑1 (MQ=255)
                                                                                                                                                                                                               tAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGt       <  1:25054/239‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑TTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                     <  1:792572/155‑1 (MQ=255)
                                                                                                                                                                                                                      aaTAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑TGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTc  <  1:371709/237‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTC  >  NC_000913/361439‑361889

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