Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,138,560 C→T G159D (GGC→GAC)  gldA ← L‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,5600CT96.0% 68.3 / ‑3.5 26G159D (GGC→GAC) gldAL‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/5);  total (20/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

CTCAAATAGGTGTTCGCTTCAATCACGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGCCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCGAGAATTGCGCCACACTGCGCAGTCTCCGCGAT  >  NC_000913/4138314‑4138801
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
cTCAAATAGGTGTTCGCTTCAATCACGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGc                                                                                                                                                                                                                                               >  1:9652/1‑251 (MQ=255)
                  tCAATCACGCGCTCCAGCGCCGGCGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACCGACGTGCAGGCGCGTCAGCGACGATTTTGGTGTCgcc                                                                                                                                                                                                                             >  2:137321/1‑249 (MQ=255)
                   cAATCACGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGAc                                                                                                                                                                                                                             >  1:111388/1‑250 (MQ=255)
                        aCGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGAc                                                                                                                                                                                                                       >  1:264300/1‑251 (MQ=255)
                        aCGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGGGCGTCAGCGACGATTGTGGTGTCGACAATGAc                                                                                                                                                                                                                       >  2:160063/1‑251 (MQ=255)
                          gcgcTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGAcc                                                                                                                                                                                                                      >  1:66967/1‑250 (MQ=255)
                                cAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCg                                                                                                                                                                                                               >  1:216925/1‑251 (MQ=255)
                                             aCTACATGCTGTTAGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGAAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCACGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCACACCAGGTTGCCAGCGCCTCGCCGATACCCGCCGCTAACAGACGTGAAGGTGCGTCAGCGCCGATTTTGGTGTCGACAATGACCATTATCGGGTTATTTGGCa                                                                                                                                                                                                   >  2:181085/1‑250 (MQ=255)
                                                                    cctCGCTTTTTAGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGACAGAGCAAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAAAAGCAGATAGCGGTCAAACTCAc                                                                                                                                                                           >  2:275112/3‑251 (MQ=255)
                                                                                                            aGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGAGGTCAAACTCACCCTCATCGGTGTAGATAACAGAAAATGCGCTGCACGGTGc                                                                                                                                   >  2:87575/1‑251 (MQ=255)
                                                                                                             gTTCCGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGCGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACACTGACCCTATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAACCTCAACCTCATCGTTGTAGATAACAGACAATGCGCTGCACGGTGcc                                                                                                                                  >  1:251723/1‑250 (MQ=255)
                                                                                                                  gCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGt                                                                                                                             >  2:247772/1‑251 (MQ=255)
                                                                                                                                  ccTGGGTGCACTTGCCTCCCGCCATGGTGGTCGCGCCGATACGAGCGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGCTACCCGCCGCTAACCGCCGTGCAGGCGCGTCAGCGACGATTTTGGTGCCGACAATGACCATATTCGGGTTATTTCGCACCAGCAGATCGCGGTCAACCTCACCCTCATCGTTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGAGATAGCCggg                                                                                                             >  2:86100/1‑250 (MQ=255)
                                                                                                                                                                   cgccgcTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATAGTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGCTCGCTACCGGAACACCCATGAAATGAGcc                                                                            >  2:79596/1‑251 (MQ=255)
                                                                                                                                                                       gcTACGAGAGCAGGCACGCGCTTCAAACCAGGCTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGAGTCAGCGACGATTTTGATGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGCAGATAACAGACAATGCGCTGCACGGGGCATCGGAAGAGGCGATAGAAAGTGCGGTCGCTACCGGAACACCCagtaaatgagccgga                                                                         >  1:252788/1‑236 (MQ=255)
                                                                                                                                                                             agagCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTgg                                                                  >  2:24138/1‑251 (MQ=255)
                                                                                                                                                                              gagCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCCAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGc                                                                 >  1:133145/1‑251 (MQ=255)
                                                                                                                                                                                                      gtgtccagcagctctcaaatccccGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCAcc                                         <  2:31538/230‑1 (MQ=255)
                                                                                                                                                                                                      gtggcgagcccatgaccGATCCCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCAcc                                         <  2:264300/236‑1 (MQ=255)
                                                                                                                                                                                                             gggccgctatgttctccgccgCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCg                                  <  2:7888/236‑1 (MQ=255)
                                                                                                                                                                                                                    gCCGATACCCGCCGCTAACAGACGTGCAGGTGCGCCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAAcagca                                                                                                                                                                                             >  2:257032/1‑89 (MQ=255)
                                                                                                                                                                                                                        aTACCCGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCAACCGATACCGAGAATTGCGcc                       >  2:154897/1‑251 (MQ=255)
                                                                                                                                                                                                                           cccGCCGCTAACAGACGTGCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCGAGAATTGCGCcac                     >  2:84181/1‑250 (MQ=255)
                                                                                                                                                                                                                               ccgcTAACAGACGTGCAGGTGCGTCAGCGACCATTTTAGTGACGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGAAGATAACAGACCATGCGCTGCACGGTGCATCCGTAGAGGCGATAGTCTGTGCGATCGCTACCGGAACGCCCATGAAATGTGCCAGTGCTTTGGCAGTATTGAGGGTTTTTCCGCCGCCGATCCCGAGAAATTCGCcaccct                  <  2:130130/249‑4 (MQ=255)
                                                                                                                                                                                                                                    aaCAGATGTGCAGGTGCGGCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCGAGAATTGCGCCACACTGCGCa             <  1:162303/249‑1 (MQ=255)
                                                                                                                                                                                                                                             gCAGGTGCGTCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCGAGAATTGCGCCACACTGCGCAGTCTCCGCGAt  <  2:133145/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
CTCAAATAGGTGTTCGCTTCAATCACGCGCTCCAGCGCCGGAGTCACTACATGCTGTTCGGCAGCAAGCATCGCTTTTTCGCCTTCTTCCAGCAGGGTGTTGTAGCACAGTTCAGCCAGTGCCAGCGCAGCCTGGGTGCACTTGCCGCCCGCCATGGTGGTCGCGCCGCTACGAGAGCAGGCACGCGCTTCAAACCAGGTTGCCAGCGCATCGCCGATACCCGCCGCTAACAGACGTGCAGGTGCGCCAGCGACGATTTTGGTGTCGACAATGACCATATTCGGGTTATTTGGCAACAGCAGATAGCGGTCAAACTCACCCTCATCGGTGTAGATAACAGACAATGCGCTGCACGGTGCATCGGTAGAGGCGATAGTCGGTGCGATCGCTACCGGAACACCCATGAAATGTGCCAGTGCTTTGGCAGTATCGAGGGTTTTTCCGCCACCGATACCGAGAATTGCGCCACACTGCGCAGTCTCCGCGAT  >  NC_000913/4138314‑4138801

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