Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,489,016 C→T E350K (GAA→AAA)  yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,0160CT93.3% 36.3 / NA 15E350K (GAA→AAA) yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/7);  minor base G (0/1);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGC  >  NC_000913/4488788‑4489241
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
gAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTgn                                                                                                                                                                                                               <  1:102115/249‑2 (MQ=255)
     ccTTTTTGACACCTCTCATAGGCGGGTTCCCGGTAAACCTCATCCTCATATTTGCACTACACCGGAGAGTGATTACTCAAGCCATTAGGCTCATCTTCCGTCACCGGCCCTATCCGCGCACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAACAGAGATCAACGCCTCGCCGGTTCCCAGTTCCTGCATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCAt                                                                                                                                                                                                        >  1:295785/1‑251 (MQ=255)
       tttttGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTAGTATCAAACGCCTGATTGGCCCGCATg                                                                                                                                                                                                       <  1:77740/250‑1 (MQ=255)
                                        cACCTCATCCTCATATTTGCCATACCCCGGAGAGCGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCACCATCCGCGAACAAGGCGAGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCCCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCtt                                                                                                                                                                     >  1:170870/1‑251 (MQ=255)
                                          cctcatcctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGTCGCGATCACCATCGAAAGCTCCACCACAGAGGGGCTTCCTTTCGCATACAGAAAAGAGATCAGCGCCTCGCCGGTGCCCAGTTCCTGAATCGACTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCtttt                                                                                                                                                                   >  1:244148/1‑251 (MQ=255)
                                              atcctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGTCGGTGCCCAGTTCCTGAAGCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCg                                                                                                                                                               <  1:237973/251‑1 (MQ=255)
                                                                  ncGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCAGTCACCGGCCCCATCCGCGAACAAGGAGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCAGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGCTCTTTTGGCGtcaa                                                                                                                                                           >  1:37307/2‑232 (MQ=255)
                                                                                                              ggCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCTGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCGTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCg                                                                                                      <  1:141543/244‑1 (MQ=255)
                                                                                                                          gAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTg                                                                                   >  1:281589/1‑251 (MQ=255)
                                                                                                                                   gcgATCACCATCGCACGCTCCACCACAGTAGGCCTTCGTTTCGGATCTACAAAAGAGATCAACTCCTCGCCCGTGCCCATTTCCTGAATCGCGTTAGGGGTATCAAACGCCGGATTGGCCCGCATGGTTTgcg                                                                                                                                                                                                <  1:237217/133‑1 (MQ=255)
                                                                                                                                                                   gCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGGAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAg                                                    <  1:158321/241‑1 (MQ=255)
                                                                                                                                                                                    gAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAATCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCTAGCTCATTATCCGGAATATCAGTCGGGTTTTGCGAAACGAACCAGACGCCCACTCCTTTTGAGCGAATAAGCCTTATCCCCTGc                               <  1:342354/245‑1 (MQ=255)
                                                                                                                                                                                        nganaTCAANGCCTNGCCGGTGCCCAGTTCNTGAAGCGCCTTTTTGGTACNNANCGCCGNATTgg                                                                                                                                                                                                               <  2:153762/64‑1 (MQ=255)
                                                                                                                                                                                                           gTGCCCAGTTCCTGATTCGCCTTTTTGGTATCGATCGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGc  >  1:119936/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCCTTTTTGGTATCAAACGCAGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAACGCATGTTGAACGCGATTACCGCGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGcc                                          >  1:290464/1‑193 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
GAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGC  >  NC_000913/4488788‑4489241

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