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,0160CT96.3% 79.1 / NA 27E350K (GAA→AAA) yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (13/13);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCT  >  NC_000913/4488773‑4489234
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          
gctCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCCCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATc                                                                                                                                                                                                                      >  2:138817/1‑250 (MQ=255)
               gAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCCCCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGcc                                                                                                                                                                                                      <  1:225702/251‑1 (MQ=255)
                       ttttGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATgg                                                                                                                                                                                               <  1:294778/250‑1 (MQ=255)
                          tGCAACATCTCATAGGCGGATTCCCGGTCCACCTCACCCTCATATTTGCCATACACCGTGGAGTGATGAATCAAGCCATTACGCTCATCTTCCGCCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCTCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTGCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGttt                                                                                                                                                                                            <  1:265663/250‑1 (MQ=255)
                                        ggCGGATTCCCGGTCCACCTCATCCTCATATGTGCCAAACCCCGGAGAGTGATTAAACAAGCCATTACGCTCACCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCa                                                                                                                                                                             <  2:71735/251‑1 (MQ=255)
                                          cGGATTCCCGGACCACCTCATCATAATATTTGCCATACACCGGAGAGTGAGTAATCACGACATTACGCTCATCTTCCGTCACCGGCCCCATCCACGAACACAGCGCGCTCACCATCACACGCTCCACCACAGAAGGGATTCCTTTCGCATCCAGGAAATAGTTAAACGACTCGCAGGTGCACAGTTCCTGAATCGCCTCTTTGGTATCAAACGCCGGATTGGCCCGAATGGATTGCGCAGCTGCCTTCAc                                                                                                                                                                            >  2:184774/1‑250 (MQ=255)
                                           ggtttCCCGGCCCACCCCATCCTCATAGTCGCCAAACCCCGTAGAGTGATTAATCAAGCCATTACGCTCATCTCCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTg                                                                                                                                                                          <  2:77158/248‑1 (MQ=255)
                                            gATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGc                                                                                                                                                                         <  1:138817/251‑1 (MQ=255)
                                                        aCCTCATCCTCATATTTGCCATACACCGGAGAGAGATTAATCAAGCCATTACGCCCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGCTCACCATCGCACGCTCCACCACAGCAGGGATTCCTTTCGCATCCAGAAAAGAGATCCACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCttt                                                                                                                                                             >  1:274308/1‑251 (MQ=255)
                                                                                         tGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAAGCGCGATCACCATCGCAAGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCAGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCACAAAGCATGTTGAAcgcg                                                                                                                            >  1:188083/1‑251 (MQ=255)
                                                                                               aTCAAGCCACTCCACTCATCTTCCGTCACCGGCCCACTCCTCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGACATCATCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAAGCGATTAcc                                                                                                                      <  2:19625/251‑1 (MQ=255)
                                                                                                                     ccGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCGCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATcc                                                                                                <  1:10916/251‑1 (MQ=255)
                                                                                                                       gTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCCAAGCATGTTGAACGCGATTACCGGGCTGCCCGAGCACATTATCCgg                                                                                              >  2:214611/1‑251 (MQ=255)
                                                                                                                          cccGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAat                                                                                           <  1:217812/250‑1 (MQ=255)
                                                                                                                                cccATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATc                                                                                        >  1:234033/1‑248 (MQ=255)
                                                                                                                                      cgagAACAAGGCGCGATCACCATCTCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGCATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCCTGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAACGCCCGCAAAGCATGTTGAACGCGATTACCGCGCTGCCCGAGCACATTATCCGGAATATCAGACGGGtt                                                                               >  2:218757/4‑251 (MQ=255)
                                                                                                                                         gAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACCTTATCCGGAATATCAGACGGGTTTTg                                                                            >  1:202688/1‑251 (MQ=255)
                                                                                                                                          aaCAAGGCGTGAGCACCATCGCACTCTCCACCACAGAAGGGCTTACTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAAGCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGc                                                                           <  2:97511/251‑1 (MQ=255)
                                                                                                                                                  gcgATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGa                                                                    >  1:253731/1‑250 (MQ=255)
                                                                                                                                                   cgATCACCATCGCCCGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGATATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCTGCATGGTTTGCGCCGCTGCGTTCACTGCTTTCTGATCTTTTTGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAAc                                                                  <  1:125227/251‑1 (MQ=255)
                                                                                                                                                       cACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGa                                                              <  2:295078/251‑1 (MQ=255)
                                                                                                                                                                 cGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGcc                                                    >  1:171884/1‑251 (MQ=255)
                                                                                                                                                                       cccacAGAAGGGCTTCCTGTCGCATCCAGAAAAGAGATCACCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGGTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGa                                              <  2:222630/250‑1 (MQ=255)
                                                                                                                                                                            aGAATGGCTTTCTTTCGAATTCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAAGCGCCTTTTTGGTAACAAACGCCGGTTTTGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAACGCCCGCAAAGCAGGTTTAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTACGAAACGAACCAGACTCCCACGCCTTTTGAGTGaa                                         <  2:253952/251‑1 (MQ=255)
                                                                                                                                                                                   cTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCt                                  >  1:70055/1‑251 (MQ=255)
                                                                                                                                                                                                            gcAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGAAAAAGCCCGCAAAGCATGTTTAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGACACGAACCAGACGCCCACGCCTTTTGAGCGAATCAGCCTTATCACCTGCTCACTCTTATCcagc         >  2:231506/2‑251 (MQ=255)
                                                                                                                                                                                                                   cTCGCCGGTGCCCAGTTACTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCAGGAATATCAGACGGGTTTTGCGAAACGCACCAGACGCACACGCATTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTCTCcagcagaaact  >  2:282855/1‑246 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          
GCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCT  >  NC_000913/4488773‑4489234

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