Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,845,026 G→A I482I (ATC→ATT)  hycE ← formate hydrogenlyase subunit HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,845,0260GA91.7% 25.7 / NA 12I482I (ATC→ATT) hycEformate hydrogenlyase subunit HycE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/7);  minor base C (1/0);  total (5/7)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CACATCGACCACGGTCATGCGGTCGGTACAGGAGTAGCAAGGGTCGAGGCTACCGATAATCAGCGGCGCATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCTGCTCGCTGTGGACTTCCATTGGCAGCAGGCCATAGCCGACAAACGGGTGATCGGCGCGGGTATCACG  >  NC_000913/2844843‑2845271
                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
cacaTCGACCACGGTCATGCGGTCGGTACAGGAGTAGCAAGGGTCGAGGCTACCGATAATCAGCGGCGCATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAatat                                                                                                                                                                                                                                                   >  1:18204/1‑188 (MQ=255)
                       cggtACAGGAGTAGCAAGGGTCGAGGCTACCGATAATCAGCGGCGCATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGGGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGc                                                                                                                                                             <  1:95283/251‑1 (MQ=255)
                                     cAAGGGTCGAGGCTACCGATAATCAGCGGCGCATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGAcc                                                                                                                                               >  1:109858/1‑251 (MQ=255)
                                          gTCGAGGCTACCGAGAATCAGCGGCGCATCGGAAACGAAGTTGCCGCGCAGCATGTAGAGCAGGGTAGGACAGTTAGCATAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCCCAGGTCAGGATCCAGTGCATATCATCGCCGCGCGGAGCTTCGGCAAAGCCCGGCGCAACGCGGTGCGGAATGACGGTAAAGCCTTCCACCATCAGGTGGCAACCCAGACGGTTATCCAGACCGt                                                                                                                                             >  1:169812/1‑248 (MQ=255)
                                                                                                      cagtgtctgCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGAGAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATcc                                                                                                                                                    <  1:139070/173‑1 (MQ=255)
                                                                                                        ggggcGGCCAGGTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGGTATCCAGan                                                                                                                                                <  1:64877/179‑2 (MQ=255)
                                                                                                             ggCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTGCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGGAGGTAAAGCCGTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCAt                                                                                                                                   <  1:86695/191‑1 (MQ=255)
                                                                                                               ccAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAACGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATc                                                                                                                                                     <  1:82298/171‑1 (MQ=255)
                                                                                                                     cgcgTAGGTCGCGGCACGGCCGTGCCAGCGGTACAGATTCTGGTTGTCGCCGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCTGCTCGc                                                               >  1:73387/1‑251 (MQ=255)
                                                                                                                                                                      ccGGTCATGCTCCAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTGATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCTGCTCGCTGTGGACTTCCATTGGcagcag                                         <  1:8008/224‑1 (MQ=255)
                                                                                                                                                                               cNCCNGTGAANATCANCGCCGCGCGGCGCTTNGGCAAAGCCCAGCGCAAAGNNGNGCGGANTGTAg                                                                                                                                                                                              <  2:91947/66‑1 (MQ=255)
                                                                                                                                                                                  cAGTGAATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCTGCTCGCTGTGGACTTCCATTGGCAGCAGGCCATAGCCGACAAACGGGTGATCGGCGCGGGTATCACg  >  1:161522/1‑251 (MQ=255)
                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CACATCGACCACGGTCATGCGGTCGGTACAGGAGTAGCAAGGGTCGAGGCTACCGATAATCAGCGGCGCATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCCTGCTCGCTGTGGACTTCCATTGGCAGCAGGCCATAGCCGACAAACGGGTGATCGGCGCGGGTATCACG  >  NC_000913/2844843‑2845271

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