Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,628,644 A→T C99S (TGT→AGT)  yfgF ← cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6440AT100.0% 65.9 / NA 18C99S (TGT→AGT) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (10/8);  total (10/8)

AGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAAT  >  NC_000913/2628425‑2628845
                                                                                                                                                                                                                           |                                                                                                                                                                                                         
aGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAACCTGTAACCGGCCCAACTGGGaatga                                                                                                                                                                                                   >  2:89110/1‑228 (MQ=255)
        tAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCCGGGAGTAAAGACCCTGTAACCGCACCAACTTAGAATGGTAATGATAATCAAATTaca                                                                                                                                                                                >  2:405471/1‑237 (MQ=255)
                                      aTGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                                                                                                                   >  1:44772/1‑238 (MQ=255)
                                                    aTTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                                                                                                    >  1:633990/1‑239 (MQ=255)
                                                                          ccaaaaatggcagaaAACGGAAAAAGTATAAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCTGGAAATTTAAACGGGTATCCCCATGTAAAACGTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                               <  2:493140/181‑1 (MQ=255)
                                                                             aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACACAAATCCCTAAGGCGATc                                                                                                             >  1:56117/1‑237 (MQ=255)
                                                                                    cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                                                                                                     >  1:493140/1‑190 (MQ=255)
                                                                                                aaaTTAGCAAAAACGGCGTTGCAGAAACCAAATCCTGCAAGAATATACGCTGGAAAATAAACCGGGTATCCCAATGAAAACCGTTTTTACGCGGGGAAAAAAAGCCCCTATAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTg                                                                                        <  2:44772/239‑1 (MQ=255)
                                                                                                    tATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaa                                                                                    >  1:733737/1‑239 (MQ=255)
                                                                                                           aaCGGCTTTGAAGGAACCCAAATCGGCAAGAATACAGGATGGAAAATTAACCGGTTACCACAATGTAAACCGTTGTTCCGCCGGGAAGAAAAGCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcat                                                                             <  2:56117/239‑1 (MQ=255)
                                                                                                                         aaaCACAAGCTGACAAAATATACGCTAGGAAATTACGCGAGTATCACCAGGTAAATGGTTGTACCGCGGGGGAATAAACACCTTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                 <  2:633990/237‑1 (MQ=255)
                                                                                                                                           tataCGCTGGGAAATTAACCGGGTATGACCATGTGAACCGTTGTAACGCGGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATaaa                                             <  1:20098/239‑1 (MQ=255)
                                                                                                                                                                 gTATCACCATGTAAAACGTTGTTACGCCGGGAAGAAAAGCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa                                                                                   <  2:38583/179‑1 (MQ=255)
                                                                                                                                                                    tctccAGGTGAATCGTGTTAACGCCGGGGAGGAAAGCCCCGTTATCCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                  <  2:504498/110‑1 (MQ=255)
                                                                                                                                                                           gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATg                                                                         >  1:38583/1‑179 (MQ=255)
                                                                                                                                                                             gAAACGTTGTTACGCCGGGAAGTAAAGCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa            <  1:405471/238‑1 (MQ=255)
                                                                                                                                                                              aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                                                                                                        >  1:504498/1‑113 (MQ=255)
                                                                                                                                                                                      ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCCGCAAAGGATTAAt  >  2:313213/1‑239 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                         
AGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAAT  >  NC_000913/2628425‑2628845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: