Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,149,958 T→C 34.7% E343G (GAG→GGG)  yegI ← protein kinase‑related putative non‑specific DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,149,9580TC34.7% 30.8 / 3.3 25E343G (GAG→GGG) yegIprotein kinase‑related putative non‑specific DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (3/13);  new base C (7/2);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 8.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAT  >  NC_000913/2149739‑2150163
                                                                                                                                                                                                                             |                                                                                                                                                                                                               
cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCCACTGTCGCCCTCACCACCGCCCCAACTTTCGCCAGCCCCAAATTCCCG‑CCGGCGGTTATAACCCCCTCGccc                                                                                                                                                                                                         >  2:486493/1‑227 (MQ=255)
         gcgtGAAGGCCGAAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACGGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                                                                                                                                                                                <  2:20122/225‑1 (MQ=255)
          ccTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCCACTGGCGCCCTCACCATCGCCCCGAACTTTGCCCGACCAAAATTTCCC‑CCGGGGGGAATGACCCCTTccccccgttcacta                                                                                                                                                                                                >  2:633893/1‑217 (MQ=255)
            taaaGGCCGAAAAGCAGCAGCAAAAATGCTGAAAGGGGGTTCTAAAGCAAAAGGATGATTCCGCGCCGAAACAGGCCTAA‑AGGAAGGGGCCTGCCAGGCGGTTGGAAATGATCGGTTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                                                                                                                                                                              <  2:137902/225‑1 (MQ=255)
                             cagcGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCCACTGACCCCATCACCATCGCCCCGACTTTCGCCAGCCACAAATTCCCC‑CCCGGGGGAAGGACCCCCTTCCCGCGACCAATAAAATAAATAaaacacttata                                                                                                                                                                            >  1:659438/1‑221 (MQ=255)
                                             aGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATTACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCC‑CCGGTGGGAATGACCTCTTCCCCGGGATCAATAAAAAAGGAAAAACCtgtctcttttacaaatttccga                                                                                                                                                            >  2:9306/1‑207 (MQ=255)
                                                         aGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACCTACGCCATCACCAACCCCCAGACTTTCACCCGCACAAAATTTCCG‑CCGGGGGGAATGACCCCTTCGCCCGGATCAATAAAAAAGATAAAACCCTTATTTTTCAGCCCAcacccacggcagacaagg                                                                                                                                                >  1:655252/1‑214 (MQ=255)
                                                                tgatgTATTCGCGCCGTAACAGGCCTAA‑GGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTGGCAATGCTGGCGGTGCTACGGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAg                                                                                                                                        <  2:86614/228‑1 (MQ=255)
                                                                tgaggTTTTCGCGCCTTAACGGGCGTAA‑AGGAAGCGGCCGGCCAGTCGGTTGGAAAGGATCGGGTAATGGCAGTTGCAATGGTGGCGGTGCTACGGACGCCACCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAg                                                                                                                                        <  2:127331/228‑1 (MQ=255)
                                                                 tgagtgttcagagcGTAAAAGGCCCAAGAGGGAGCGCCCTGCCAGGGGGTGGAAAATGACGGGGTAATGGCATTTGCGATGTGGGCGGTGCTCCTGACCCCATCACCCCCGCCAGGTCTTCCGCCAGCAAAAATTTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAg                                                                                                                                        <  1:372019/215‑1 (MQ=255)
                                                                 ggttttgtggcgcCGTAACAGGCCTAA‑AGGAAGCGGCCGGCCAGTCGTTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGGGGTGCTACTGACGCCATCACCATCGCCCAGACTTCCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGc                                                                                                                                       <  2:167344/219‑1 (MQ=255)
                                                                             ccGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCCTCGCCCAGACCTTTGCCCGCACAAAAATCCCG‑CCCGGGGGAATGACCTCTTCGCCCAGATCAATAAAAAAAGAAAAACCCTGTTTTACCAGCGCACACCACCGGGCGA‑CCATCAAAAAAGCGGGGGAAACaac                                                                                                                           >  1:565792/1‑226 (MQ=255)
                                                                                tAACAGGCCAAA‑AGGAAGCGGCCTGCCGGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                                                                                                                                                                              <  2:354793/159‑1 (MQ=255)
                                                                                               gacggggCCTGCCGGTCGGTTGGAAATGATCGGGAAAGGGCAGTTGAAAGGCTGGGGGTGCAACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATAc                                                                                                          <  1:216940/223‑1 (MQ=255)
                                                                                                 gcgcgccctgcgtggggTTGGAAATGACGGGGTAATGGCGGTGGCAAGGCTGTGGGGCCTACGGCCCCATTCACCATCCCCCAGACTGTCGCAGGAAAAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACAt                                                                                                        <  1:633893/213‑1 (MQ=255)
                                                                                                 aGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGCGCTCCTGACCCCCCCCCCATCCCCCCAACATTCGCCCGCACAAAATTCCCC‑CCGGGGGGAATGACCCCCTCGCCGCGGTCAACAAAAAAAAAAAAACCCCGCTTATCCCGCCCCCCCCACCGGCAATCCAGACAAA‑GATCCGGGGAAACAAGCAATGcccacaacac                                                                                                            >  1:136864/1‑215 (MQ=255)
                                                                                                                            tCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGAAGCCATCACCATCCCCCAGACGTTCGCCCGCACAAAAATTCCC‑CCGGTGGTAATGACCCCCTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCACCGCCCACCCCAGGCCGG‑CAGTCAAAAAAAAGGGGGAAAAATTAATTGCCGAAAACAAACATTTT‑TTTTAcatgtgggggagaaaagcc                                                                            >  1:423868/1‑210 (MQ=255)
                                                                                                                              ggtgctgtgcggttacagtcgcggCTGTGCTATGGCCGCTATCCCCCTCGCCCAGACTTCCGCCAACAAAAATTCCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCa                                                                           <  1:14099/205‑1 (MQ=255)
                                                                                                                                        aGTTGCAAGGCTGGCGTTGCTACTGCCGCCATCACCATCGCCCAGACTGTCGCCAGCCAAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTAcc                                                                <  2:229356/228‑1 (MQ=255)
                                                                                                                                                          gCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGGGGTAATGACCCCTTCGCCGAGGTCAAAAAAAAAAATAACACCTTGATTATCCCACGCCCACCACGGGCAGA‑CAGGCAAAAAAGCGGGGGAAAAATGCCTTGCTGAAACGGTAACTTTCTTTTAATTTGTGGCTAAGAGAGTCC‑TCGCGGCACCCCACCCTTTAGCCCTCgg                                               >  1:375351/1‑227 (MQ=255)
                                                                                                                                                                                         tCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCTTGaagt               >  2:349581/1‑226 (MQ=255)
                                                                                                                                                                                                    aaaTTCCCG‑CCGGTGGTAATGACCCCTTCGCCGAGATCAATAAAATAGATAACACCGTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAAc    <  2:136864/228‑1 (MQ=255)
                                                                                                                                                                                                     aaatCCCGCCCGGTGGTAATGACCTCTTAGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAAc    <  2:655252/225‑1 (MQ=255)
                                                                                                                                                                                                      aTTCCCG‑GCGGTGGTAATGACCCTTGCGGCGAGATCAATAAAATAGATAACCACTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACGGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAt  <  2:295944/228‑1 (MQ=255)
                                                                                                                                                                                                      aTTCCCG‑GCGGGTGTATAGCCCTTCTCGCCGGGAACAATAAAATAGATAACCCCTGGATTAACCAGGGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACAAGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAt  <  1:486493/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                               
CGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAA‑AGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCG‑CCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGT‑CAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTC‑TTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAT  >  NC_000913/2149739‑2150163

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