Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 372,966 A→G T16A (ACC→GCC)  mhpF → acetaldehyde‑CoA dehydrogenase II, NAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913372,9660AG100.0% 76.8 / NA 22T16A (ACC→GCC) mhpFacetaldehyde‑CoA dehydrogenase II, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/11);  total (11/11)

ATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372857‑373099
                                                                                                             |                                                                                                                                     
atagTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGCTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCAc                                                                                                          >  2:539273/1‑139 (MQ=255)
    tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:174179/139‑1 (MQ=255)
    tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:519218/139‑1 (MQ=255)
    tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:799819/139‑1 (MQ=255)
      aaGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    <  2:196274/139‑1 (MQ=255)
        ggCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               >  2:91113/1‑110 (MQ=255)
        ggCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               <  1:91113/110‑1 (MQ=255)
             aGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTgg                                                                                             >  1:214754/1‑139 (MQ=255)
                        ctgcGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  2:644040/139‑1 (MQ=255)
                          gcGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGAt                                                                                >  2:234006/1‑139 (MQ=255)
                                                       aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATgg                                                                                       <  1:65343/103‑1 (MQ=255)
                                                       aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATgg                                                                                       >  2:65343/1‑103 (MQ=255)
                                                       aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTc                                                   >  2:14041/1‑139 (MQ=255)
                                                                     tAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGa                                     >  2:650892/1‑139 (MQ=255)
                                                                        gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:387300/71‑1 (MQ=255)
                                                                        gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:387300/1‑71 (MQ=255)
                                                                           tAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  1:484018/88‑1 (MQ=255)
                                                                           tAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  >  2:484018/1‑88 (MQ=255)
                                                                                                       aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGga                                                                                            <  2:144669/50‑1 (MQ=255)
                                                                                                       aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGga                                                                                            >  1:144669/1‑50 (MQ=255)
                                                                                                       aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGa   >  1:176686/1‑139 (MQ=255)
                                                                                                        ttGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAt  <  2:176686/139‑1 (MQ=255)
                                                                                                             |                                                                                                                                     
ATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372857‑373099

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