Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 841,944 C→T G38D (GGT→GAT)  mcbA ← DUF1471 domain‑containing protein McbA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913841,9440CT100.0% 18.1 / NA 10G38D (GGT→GAT) mcbADUF1471 domain‑containing protein McbA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/5);  total (5/5)

TTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACG  >  NC_000913/841734‑842161
                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
ttttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTGTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTATCGGCGGGCCGTAGCTGACCg                                                                                                                                                                                                       <  1:62836/231‑1 (MQ=255)
                                                                                                             cGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTATCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCGCAGGAc                                                                                                                                            <  1:90674/181‑1 (MQ=255)
                                                                                                                                   cgcgCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTTCCGAAACGGTATCGTCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCgg                                                                                                                                                                                   <  1:123443/120‑1 (MQ=255)
                                                                                                                                                             ccAGTTTAGCCTCCAGATCGCTTAGGTTCTATGCGCCAGTTGCCGAAACGGTGTCGGCGGGGCGTAGCTGACCGCTGACCAGGTTACTCATCTGTTGCGCGGCATAAGCCGGGAGAGAGATTCCGCTCAGGACGGTTGCAACCAGTATTGTGAGGCACTGTTTCAGAATCACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACg                                      <  1:111210/235‑1 (MQ=255)
                                                                                                                                                             ccAGTNTCTGCNCNNGGTCGCTTAGTCTTCATGCNCCAGTTGCCGAAACGNTATCNGCGGGnccna                                                                                                                                                                                                               >  2:111210/1‑63 (MQ=255)
                                                                                                                                                                       ctccAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAACCGGTATCGGCGGGGCGTAGCTGACCGCTGTCAAGGTTTCTCATCGGtttcgg                                                                                                                                                                             >  1:228171/1‑86 (MQ=255)
                                                                                                                                                                              tCGCTTATGTTTGATGCGCCAGTTGCCGAACCGGTATCGTCTGGCCGTATCTGCCCGCTGGCCCTGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGTACGGTGGCAATCAGTAGTGTGAGTCACTTTTTCATAATTACATACTCATAAATACAACGGGTAGTGTCTGATTATTGTACGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGAtgt                  <  1:161888/238‑1 (MQ=255)
                                                                                                                                                                                 cTTAGGTTTGATGCGCCAGTTGCCGAAACGGTATCGGCGGGCCGTAGCTGACCGCTGTCCCGGTTACTCATCGGTTGCGCGGCATAAGCCGTGTGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACg  >  1:15282/1‑251 (MQ=255)
                                                                                                                                                                                         tGATGCGCCAGTTGCCGAAACGGAATCGGCGGGCCGTAGATGACCGCTGCCCAGGCTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGCCGGTGGCAATCCGTAGTGTGAGGCACTTTTTTATAATGACATCCTCTTAAATACAACGGGTAGTGCCTGATAAGTGTAGGTAGCCTGCGCGATGCAAATGACg                                      >  1:226301/1‑207 (MQ=255)
                                                                                                                                                                                               gcCAATTGCCGAAACGGTATCGGCGGGCTGTAGCTGACCGCCGTCCCGGGTACTCACCGGTTGCCCGGCATAAGCCGTGAGAGAGAttnncccc                                                                                                                                                 >  2:12160/1‑88 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
TTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTAGGTTTGATGCGCCAGTTGCCGAAACGGTACCGGCGGGCCGTAGCTGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACG  >  NC_000913/841734‑842161

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