Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,116,753 C→T R157H (CGT→CAT)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,116,7530CT95.0% 52.3 / ‑5.3 21R157H (CGT→CAT) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/11);  total (8/13)
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.95e-01

GTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTG  >  NC_000913/4116516‑4116996
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
gTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGc                                                                                                                                                                                                                                                 <  1:352279/242‑1 (MQ=255)
   ccgccTTTGCCGCCAATGTTAGTATGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTcgcg                                                                                                                                                                                                                                        >  1:364269/1‑248 (MQ=255)
                 aaTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATGTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCGTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCa                                                                                                                                                                                                                                                <  1:249537/226‑1 (MQ=255)
                                          gaagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGATGATCCACTTCACTTTGGTGcc                                                                                                                                                                                              <  1:340753/251‑1 (MQ=255)
                                                                        cgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGCTGAACAACATGGTAAGAGCGGCGTTGGAGTAATCGGTCACATGGACACGGCCCTGAGTCTTTTTCCAGCTAAGCCACGCATCAACAGTACCAATGAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGCCGAGGATCCACTTCACTTTGGTGTCAGAAAAGTACGGGTCAATCACCAGACCgg                                                                                                                                                                 >  1:441864/1‑250 (MQ=255)
                                                                                         ttCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGAtgtg                                                                                                                                                                                                                                                  <  1:300495/152‑1 (MQ=255)
                                                                                                                   gggTATGGATGTTGAACAACATGGTACGAGAGGAGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcact                                                                                                                                                                                                           >  1:126664/1‑165 (MQ=255)
                                                                                                                       aTGGATGTTGAACAACATTGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAGTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATTGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTTAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTcn                                                                                                                  <  1:83411/250‑2 (MQ=255)
                                                                                                                       aTGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACGTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGc                                                                                                                 <  1:395380/251‑1 (MQ=255)
                                                                                                                                         ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGAtgtg                                                                                                                                                                                                                                                  >  1:340225/1‑104 (MQ=255)
                                                                                                                                               agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgcn                                                                                                                                                       <  1:129562/189‑2 (MQ=255)
                                                                                                                                                      ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAGTTTCCATATATGCCACGTATCAACCGTACCAAACAGCAATTGACCACGATGGGCACGCTCGCGAGAGCCTTCCACAGGGTCGAGGATCCACTGCACTTTGGTGCCAGAAAAg                                                                                                                                                                                       <  1:180193/150‑1 (MQ=255)
                                                                                                                                                                  gTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCTACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCGCCAGACCGGTATTGCGGCGGATATAATCTTCTAAACCGGCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGcn                                                                           <  1:55790/246‑2 (MQ=255)
                                                                                                                                                                                gCCCTGAGTCATTGTCCAGATAAGCCACGTATTAACCGGACCAAACATCAATTCACCACGAGGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                                                                                                                                                                                  <  1:57576/129‑1 (MQ=255)
                                                                                                                                                                                 cncTGAGTCATTTTCCAGCTAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTAGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTatag                                                                     >  1:349242/1‑237 (MQ=255)
                                                                                                                                                                                       gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                                                                                                                                                                                >  1:25322/1‑124 (MQ=255)
                                                                                                                                                                                                       gCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGagag                                                                                                                                                                                                                                    >  1:444596/1‑56 (MQ=255)
                                                                                                                                                                                                       gCCACGTATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGCGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGACAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTc                                              >  1:369547/1‑238 (MQ=255)
                                                                                                                                                                                                             tATCAACCGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCTAGAGCCTTCCACATGTTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCg                                                                                                                                                                  <  1:318733/116‑1 (MQ=255)
                                                                                                                                                                                                                  nccGTACNNAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                                                                                                                                                                                <  2:25322/96‑1 (MQ=255)
                                                                                                                                                                                                                   ccGTACCAAACAGCAATTCACCACGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGTTTTCACGCTgg                     <  1:375337/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCNAAACAGCAATTCACCTCGACATGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCGTGTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAATCCg                                                                                                                                    <  2:263779/136‑1 (MQ=255)
                                                                                                                                                                                                                                      accacGATGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTg  >  1:168792/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
GTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTG  >  NC_000913/4116516‑4116996

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