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

CGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAT  >  NC_000913/2149739‑2150163
                                                                                                                                                                                                                             |                                                                                                                                                                                                               
cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCCACTGTCGCCCTCACCACCGCCCCAACTTTCGCCAGCCCCAAATTCCCGCCGGCGGTTATAACCCCCTCGccc                                                                                                                                                                                                         >  2:486493/1‑227 (MQ=255)
         gcgtGAAGGCCGAAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACGGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                                                                                                                                                                                <  2:20122/225‑1 (MQ=255)
          ccTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCCACTGGCGCCCTCACCATCGCCCCGAACTTTGCCCGACCAAAATTTCCCCCGGGGGGAATGACCCCTTccccccgttcacta                                                                                                                                                                                                >  2:633893/1‑217 (MQ=255)
            taaaGGCCGAAAAGCAGCAGCAAAAATGCTGAAAGGGGGTTCTAAAGCAAAAGGATGATTCCGCGCCGAAACAGGCCTAAAGGAAGGGGCCTGCCAGGCGGTTGGAAATGATCGGTTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                                                                                                                                                                              <  2:137902/225‑1 (MQ=255)
                             cagcGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCCACTGACCCCATCACCATCGCCCCGACTTTCGCCAGCCACAAATTCCCCCCCGGGGGAAGGACCCCCTTCCCGCGACCAATAAAATAAATAaaacacttata                                                                                                                                                                            >  1:659438/1‑221 (MQ=255)
                                             aGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATTACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCCCCGGTGGGAATGACCTCTTCCCCGGGATCAATAAAAAAGGAAAAACCtgtctcttttacaaatttccga                                                                                                                                                            >  2:9306/1‑207 (MQ=255)
                                                         aGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACCTACGCCATCACCAACCCCCAGACTTTCACCCGCACAAAATTTCCGCCGGGGGGAATGACCCCTTCGCCCGGATCAATAAAAAAGATAAAACCCTTATTTTTCAGCCCAcacccacggcagacaagg                                                                                                                                                >  1:655252/1‑214 (MQ=255)
                                                                tgatgTATTCGCGCCGTAACAGGCCTAAGGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTGGCAATGCTGGCGGTGCTACGGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                                                                                        <  2:86614/228‑1 (MQ=255)
                                                                tgaggTTTTCGCGCCTTAACGGGCGTAAAGGAAGCGGCCGGCCAGTCGGTTGGAAAGGATCGGGTAATGGCAGTTGCAATGGTGGCGGTGCTACGGACGCCACCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                                                                                        <  2:127331/228‑1 (MQ=255)
                                                                 tgagtgttcagagcGTAAAAGGCCCAAGAGGGAGCGCCCTGCCAGGGGGTGGAAAATGACGGGGTAATGGCATTTGCGATGTGGGCGGTGCTCCTGACCCCATCACCCCCGCCAGGTCTTCCGCCAGCAAAAATTTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAg                                                                                                                                        <  1:372019/215‑1 (MQ=255)
                                                                 ggttttgtggcgcCGTAACAGGCCTAAAGGAAGCGGCCGGCCAGTCGTTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGGGGTGCTACTGACGCCATCACCATCGCCCAGACTTCCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGc                                                                                                                                       <  2:167344/219‑1 (MQ=255)
                                                                             ccGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCCTCGCCCAGACCTTTGCCCGCACAAAAATCCCGCCCGGGGGAATGACCTCTTCGCCCAGATCAATAAAAAAAGAAAAACCCTGTTTTACCAGCGCACACCACCGGGCGACCATCAAAAAAGCGGGGGAAACaac                                                                                                                           >  1:565792/1‑226 (MQ=255)
                                                                                tAACAGGCCAAAAGGAAGCGGCCTGCCGGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                                                                                                                                                                              <  2:354793/159‑1 (MQ=255)
                                                                                               gacggggCCTGCCGGTCGGTTGGAAATGATCGGGAAAGGGCAGTTGAAAGGCTGGGGGTGCAACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATAc                                                                                                          <  1:216940/223‑1 (MQ=255)
                                                                                                 gcgcgccctgcgtggggTTGGAAATGACGGGGTAATGGCGGTGGCAAGGCTGTGGGGCCTACGGCCCCATTCACCATCCCCCAGACTGTCGCAGGAAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACAt                                                                                                        <  1:633893/213‑1 (MQ=255)
                                                                                                 aGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGCGCTCCTGACCCCCCCCCCATCCCCCCAACATTCGCCCGCACAAAATTCCCCCCGGGGGGAATGACCCCCTCGCCGCGGTCAACAAAAAAAAAAAAACCCCGCTTATCCCGCCCCCCCCACCGGCAATCCAGACAAAGATCCGGGGAAACAAGCAATGcccacaacac                                                                                                            >  1:136864/1‑215 (MQ=255)
                                                                                                                            tCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGAAGCCATCACCATCCCCCAGACGTTCGCCCGCACAAAAATTCCCCCGGTGGTAATGACCCCCTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCACCGCCCACCCCAGGCCGGCAGTCAAAAAAAAGGGGGAAAAATTAATTGCCGAAAACAAACATTTTTTTTAcatgtgggggagaaaagcc                                                                            >  1:423868/1‑210 (MQ=255)
                                                                                                                              ggtgctgtgcggttacagtcgcggCTGTGCTATGGCCGCTATCCCCCTCGCCCAGACTTCCGCCAACAAAAATTCCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCa                                                                           <  1:14099/205‑1 (MQ=255)
                                                                                                                                        aGTTGCAAGGCTGGCGTTGCTACTGCCGCCATCACCATCGCCCAGACTGTCGCCAGCCAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTAcc                                                                <  2:229356/228‑1 (MQ=255)
                                                                                                                                                          gCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGTAATGACCCCTTCGCCGAGGTCAAAAAAAAAAATAACACCTTGATTATCCCACGCCCACCACGGGCAGACAGGCAAAAAAGCGGGGGAAAAATGCCTTGCTGAAACGGTAACTTTCTTTTAATTTGTGGCTAAGAGAGTCCTCGCGGCACCCCACCCTTTAGCCCTCgg                                               >  1:375351/1‑227 (MQ=255)
                                                                                                                                                                                         tCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCTTGaagt               >  2:349581/1‑226 (MQ=255)
                                                                                                                                                                                                    aaaTTCCCGCCGGTGGTAATGACCCCTTCGCCGAGATCAATAAAATAGATAACACCGTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAAc    <  2:136864/228‑1 (MQ=255)
                                                                                                                                                                                                     aaatCCCGCCCGGTGGTAATGACCTCTTAGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAAc    <  2:655252/225‑1 (MQ=255)
                                                                                                                                                                                                      aTTCCCGGCGGTGGTAATGACCCTTGCGGCGAGATCAATAAAATAGATAACCACTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACGGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAt  <  2:295944/228‑1 (MQ=255)
                                                                                                                                                                                                      aTTCCCGGCGGGTGTATAGCCCTTCTCGCCGGGAACAATAAAATAGATAACCCCTGGATTAACCAGGGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACAAGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAt  <  1:486493/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                               
CGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCGTGAATGCCTGCTGAAACAT  >  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: