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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6700AT18.9% 115.5 / 5.9 70V64E (GTA→GAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (12/41);  new base T (10/4);  total (22/48)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTC  >  NC_000913/361439‑361889
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑Tccccaga                                                                                                                                                                                                                              >  1:118795/1‑235 (MQ=255)
ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGCTGGGCGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAAACGAAAAAGAAACAAACAGGCGGT‑TTT‑ATCCacac                                                                                                                                                                                                                              >  1:379530/1‑236 (MQ=255)
          tgTTAATGCCAATCGTTATCGGAAAAGAGTTCATCTCGGCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTccc                                                                                                                                                                                                                     <  1:975179/236‑1 (MQ=255)
              aTTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATTGATGTTGGAACCGTAAGAGAAATAGAAAGGCGGA‑TCTCACCAAGGG‑GTTTTtccacacc                                                                                                                                                                                                                >  2:894396/1‑234 (MQ=255)
                     aTCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCTTCCTATATGGCTTTTGGAACCTTAAGTGAAATAGCCATGCGGG‑GTT‑TCCCCGGC‑TTTTTccccccacccgggaaa                                                                                                                                                                                                        >  2:500311/1‑227 (MQ=255)
                       tgTTATCGGAAAAGGGTACATCTCGCCGTTTTTGCTCAATTCATGGTGTACAGTGTGTCCCGTAACGGAAAGAGTAAAGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACa                                                                                                                                                                                                   <  2:853359/238‑1 (MQ=255)
                                       tACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGCCGGG‑GCT‑TCCCCGCG‑GTGTTTCCCTC‑CCCTGGGA‑AAACATTT‑TTTAATAAggta                                                                                                                                                                                   >  2:691349/1‑237 (MQ=255)
                                                                     tgtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATTAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGa                                                                                                                                                   <  1:507610/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCATTACGTTATCGCCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGAt                                                                                                                                              <  1:919086/239‑1 (MQ=255)
                                                                          aGGGTGTCCCGTAAAGGAAAGGAGAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGGTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGAt                                                                                                                                              <  1:408838/239‑1 (MQ=255)
                                                                            ggTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCCCGC‑GGTTTCCCCCC‑CCGC‑GGC‑AAAAAATTCTTTTATCCAGCATTTTCTTTTTTTAACTTCTG‑TTGG‑ATGGGCGTGGTTa                                                                                                                                            >  2:931228/1‑239 (MQ=255)
                                                                               gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTATTT‑TCCCAGCG‑CTTTTCCCCT‑‑CCGT‑GGC‑CCACATTTCTTTTATATGACGTTCTCTTTTTTTAACtttttttcgtagcggggagtataac                                                                                                                                           >  2:367870/1‑216 (MQ=255)
                                                                               gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTA‑‑AGAAATAGACAGGCGGG‑TCT‑ACCCAGGCGGTTTTCCCCT‑‑CCGT‑GGCAAAACATTTTTTTAAT‑AGGATTTCTCTATTTTCAACTTTCG‑TTGG‑ATTGGAGttgatatact                                                                                                                                        >  1:922358/1‑229 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTTTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCCA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTc                                                                                                                                       <  1:847511/239‑1 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGCGAAATAGACAGGCGGT‑TTT‑ACCCACAC‑GCTTTCCCCCC‑CCTC‑GGC‑AAACAATTCTTTTTATAGGATTTTTCTTTTTTTTACTTTTgtgtgttaggggtggtaaaaaca                                                                                                                                         >  2:630355/1‑217 (MQ=255)
                                                                                 cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTc                                                                                                                                       <  2:894630/239‑1 (MQ=255)
                                                                                               agTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCCAGC‑GGTTTTCCCCA‑CCCC‑GGC‑AAACATTTCTTTTTTTAGGCTTTTTCTTTTTTTTCATTCTC‑TTGG‑ATGAGAGTGATTAAACTCCTAAA‑TTaaatttt                                                                                                                        >  2:746293/1‑234 (MQ=255)
                                                                                                 tAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATATACAGGCGGG‑TTT‑ACCCCCGC‑GGTTTCCCCTT‑CCCTGGGC‑AAAAATTT‑TTTTATTAGGCGTTTTCTTTTTTTTACTTCTG‑TTGG‑TTTGTAGTGTTTTAAATCATCAA‑ATAAA‑Gtttt                                                                                                                      >  2:782579/1‑238 (MQ=255)
                                                                                                       ttggTGCAATCTGTTCTTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGa                                                                                                               <  2:385819/237‑1 (MQ=255)
                                                                                                         ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GGTTTCCCCTT‑CCGT‑GGC‑AAACATATCTTTAATCAGGCTTTCTCTTTTTTTAAATTCTG‑GTGG‑CTGCGAGGTGTTAAACACATACC‑TTTAA‑GTTTTTCCCCgcaa                                                                                                             >  1:842942/1‑236 (MQ=255)
                                                                                                                tCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TTT‑AACCAAGC‑GTTTTTCCCCT‑CCGG‑GGC‑AAACATTTTTTTAATTAGGCTTTTCTTTTTTTTAACTTCTG‑ATGGTATGCGGGGGTTTAAACCCATA‑A‑TTAAC‑CATTTTCCCCAAAGGtttttt                                                                                                      >  1:299954/1‑239 (MQ=255)
                                                                                                                     aCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TTT‑AACCAGGCGGTTTTCCCCAA‑CCG‑‑GGG‑AAACAATTCTTTTATCAGGCTTTCTTTTTTTTTACCTTCTT‑ATTG‑GTGCGCGGCGTTAAAATCAACAA‑TTAAA‑GTTTTTCCACCAAGACTTTTTTTCgg                                                                                                 >  1:550415/1‑239 (MQ=255)
                                                                                                                       gTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTa                                                                                               <  1:88016/239‑1 (MQ=255)
                                                                                                                         tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCCAGC‑GTTTTTCCCCC‑CCCT‑GGC‑CAAAATTTCTTTATTACAGCTTTTTCTTTTTTTTACTTCTT‑TTGG‑ATTGGCGTGGTTTAAATCAAACA‑TTTAA‑TTTTTTCCCCCAAATCttttttttcttaa                                                                                              >  2:88016/1‑231 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TTCCAAGC‑GTTTTTCCCCA‑CCGC‑GGC‑AAACATTTCTTTTATTCAGCTTTCCTTTTTTTTCACTTTTT‑ATTG‑ATTCGCGTGATTAAAATCCTACATTTAAC‑G‑TTTCCCCCCCAGTCCTTTTTTCGGTAAACCtt                                                                                         >  1:358087/1‑238 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCCGGC‑GTGTTTCCCTA‑CCGT‑GGC‑AAAAATTTCTTTTATCACGCTTTCTCTTTTTTTCACTTCTT‑CTGT‑AGTCGTGTGTTTAAAATCATACACTTAAC‑‑TTTTTCCACCAACGCTTTTTTTCGGTTAACCTTc                                                                                        >  2:666559/1‑239 (MQ=255)
                                                                                                                              ccGATTGTTACCGTGTAGTCATAGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGA‑TCT‑ACCCTCGCGGTTTCCCCCTC‑CCGC‑GCA‑AAACTTTATTTTTAT‑AGGCCTTTTCTTTTTTTCACTATTG‑ATGG‑ACGCAGGGGGTTATACTCAAACA‑TTTACTTTTTTTCCCCACAGTCTTTTTTcgttaaacatt                                                                                          >  2:504585/1‑228 (MQ=255)
                                                                                                                                 tttttttccgtgttggtaacGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                     <  2:550415/220‑1 (MQ=255)
                                                                                                                                 atggTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCCGGA‑GGTTTCCCCCACCCGT‑GCC‑AAA‑AATTCTTTTATCAGGCTTTCTCTTTTTTTCACTTCTGAATGG‑ATGAG‑GGGATTTAAATCCTACATTTACC‑GTTTTCCCCCCAAGTCTTTTTTCGGGAAACCTTTcc                                                                                       >  1:579098/4‑237 (MQ=255)
                                                                                                                                 aTTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAAGC‑GTTTTTCCCTA‑CCCG‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTTAACTTCTG‑TTTG‑TGGCCACGGATTAAACTCCTAAG‑TTAAC‑ATTTTTCCCCCAAGGCTTTTTTTCCGTTAGCCTTCgac                                                                                     >  1:428813/1‑237 (MQ=255)
                                                                                                                                   tttttccGTGTAGTCAATGACAATTGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTGCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  2:505559/234‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:746293/239‑1 (MQ=255)
                                                                                                                                   tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:329366/239‑1 (MQ=255)
                                                                                                                                      tACCGTGGAGTCATCGACAATGGTTAAATTGAAATTGGCATAAAAATGGCGGCCTATATGGATGGTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACatat                                                                                <  1:500311/239‑1 (MQ=255)
                                                                                                                                        ccGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGcaca                                                                                   <  1:782575/234‑1 (MQ=255)
                                                                                                                                                      aaaaTGGTTAAATTGAAATTTGCATAAAAATTGGGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACCCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  2:413753/237‑1 (MQ=255)
                                                                                                                                                       gattGGTTAAATTGAAATTTGCAAAAAAATTGCGGCGTATAAGGATGTTGGAACCGTAAGGGAAATAGTCAGGCGGG‑TGT‑ACCCTGGC‑GTTTGCGTCTT‑CCTG‑GGC‑AAAAATTTCGTTAATCAGGCTTTCTCTTTTTTCAACGTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTACCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCTCATATCGGTAAATAGCTTGcc                                                                <  2:507874/235‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACtt                                                                                                                                                                  <  1:210597/144‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:894396/238‑1 (MQ=255)
                                                                                                                                                       cAATGGTTAAATTGAAATTGGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGcc                                                                <  1:1083305/238‑1 (MQ=255)
                                                                                                                                                           tgTTAAATTGAAATTTGCATAAAAATTGCGGCCTTTATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGctc                                                           <  1:596445/238‑1 (MQ=255)
                                                                                                                                                                    gAAATTTGCATAAAAATTGCGGCCTATAGCGATGTTGGAACCGTAAGAGAAAAAGTATGGCGGT‑TCT‑ACCCATGC‑GTTTTCCCCTA‑CCGT‑TGC‑AAACATGTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCttt                                                  <  1:504585/239‑1 (MQ=255)
                                                                                                                                                                        attGAATAAAAATTTGGGACTATAGGGATGTTGGAACCGTA‑‑AGAAATGTACAGGCGGT‑TCT‑CCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:922358/238‑1 (MQ=255)
                                                                                                                                                                          ttcaaagaaaTTGGGGCCTATATGGATGTTTGACCCGTAAGAGAAATAGACAGGCGGT‑GCT‑ACCCAGGC‑GTTTTCCCCTT‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGGCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:118795/232‑1 (MQ=255)
                                                                                                                                                                          tgactAAAAATTGCGGCCTATTTGGTTGTTGGAACCGTAAGAGAAATAGACCGGCGGT‑TCT‑CCCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  2:75292/235‑1 (MQ=255)
                                                                                                                                                                          tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                            <  1:533688/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGGGGGCTTTATTGATGTTGGACCCCGAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  2:569451/239‑1 (MQ=255)
                                                                                                                                                                               aaaaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                       <  1:931228/239‑1 (MQ=255)
                                                                                                                                                                                  aaTTGCGGCCTATTTGGATGTTTGAACCGTAAGAGAAATAGACAGGCTGT‑TTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATg                                    <  2:229719/239‑1 (MQ=255)
                                                                                                                                                                                         gCCTATATGGATGTTGGAACCGTAAGAGAAAAAGTCAGGCGGT‑TTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAt                              <  2:381648/238‑1 (MQ=255)
                                                                                                                                                                                             atgtGGATTTTGGAACCGTAAGGGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:560783/236‑1 (MQ=255)
                                                                                                                                                                                             atatGGATGTTTGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑CACCAGGC‑GTGTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTGGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                         <  1:409652/239‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑cc                                                                                                                                                                                                               <  2:415234/62‑1 (MQ=255)
                                                                                                                                                                                              tatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTc                                                 <  1:861215/214‑1 (MQ=255)
                                                                                                                                                                                               tgggttttttggaCCGTAAGGGAAATAGCCAGGCTGT‑TTT‑TACCATGC‑GTTTTCCCCGA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTTTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  2:379530/227‑1 (MQ=255)
                                                                                                                                                                                               tggggtttttgAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑TCCCTGGC‑TTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCCGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  2:985772/229‑1 (MQ=255)
                                                                                                                                                                                               atGGATGTTGGAACCGTAAGAGAAATAGACCGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTT‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                       <  1:813402/239‑1 (MQ=255)
                                                                                                                                                                                                 ggATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                     <  1:889057/239‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCTTAAGAGAAATAGTAAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCACCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACGTCTG‑ATTG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCTGTAAGCCTTCGCTCATATCGGTAAATAGCTTGTCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTcc                   <  2:470274/236‑1 (MQ=255)
                                                                                                                                                                                                    ttttGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:579098/237‑1 (MQ=255)
                                                                                                                                                                                                     ttttgAACCGTTAGGGAAATATACCGGTGGT‑TTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  2:474976/234‑1 (MQ=255)
                                                                                                                                                                                                     gTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                  <  1:674025/238‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGGAAAATAGACAGGCGGT‑TCC‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTa               <  2:435365/239‑1 (MQ=255)
                                                                                                                                                                                                       tGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGtt                <  1:933058/238‑1 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGT‑TTT‑TCCCAGGC‑GGTTTTCCCTA‑CCCT‑GGC‑caaaa                                                                                                                                                                                                   >  1:415234/1‑60 (MQ=255)
                                                                                                                                                                                                        ggAACCGTAAGAGAAATAGACAGGCGGT‑TTT‑ACCCAGGC‑GTTTTCCCCTA‑CCGG‑GGC‑AAACATTTCTTTTATCAGGCTTTCTCTTTTTTTAACTTCTG‑ATGG‑ATGCGAGGGATTAAACTCATACA‑TTAAC‑GTTTTTCCCCGAAGTACTTTTTTCGGTAAGCCTTTCCACATATCGGTAAATAGCTTTCCTGCTCTTTTTTTTTCTGTCCTTTGCATTTTCCATGCGAGCACTCCGtttt              >  1:591705/1‑237 (MQ=255)
                                                                                                                                                                                                            ccGTAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGata          <  1:64759/239‑1 (MQ=255)
                                                                                                                                                                                                               tAAGAGAAATAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGt       <  1:25054/239‑1 (MQ=255)
                                                                                                                                                                                                                     aaaTAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑GGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑TTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca                                                                                     <  1:792572/155‑1 (MQ=255)
                                                                                                                                                                                                                      aaTAGACAGGCGGT‑TCT‑ACCCAGGC‑GTTTTCCCCTA‑CCGT‑TGC‑AAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTG‑ATGG‑ATGCGAGTGATTAAACTCATACA‑TTAAC‑GTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTc  <  1:371709/237‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                   
CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGT‑TCT‑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: ‑