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% 142.4 / NA 40T16A (ACC→GCC) mhpFacetaldehyde‑CoA dehydrogenase II, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/20);  total (20/20)

ATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372842‑373099
                                                                                                                            |                                                                                                                                     
aTCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGAtt                                                                                                                         >  1:378039/1‑139 (MQ=255)
           gCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCg                                                                                                              >  2:1074619/1‑139 (MQ=255)
             ccATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTc                                                                                                            <  1:223819/139‑1 (MQ=255)
                   tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:266470/139‑1 (MQ=255)
                     aaGGAATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    >  2:322992/1‑139 (MQ=255)
                       ggCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCa                                                                                                  >  1:389882/1‑139 (MQ=255)
                          aTAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               <  2:538783/107‑1 (MQ=255)
                          aTAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               >  1:538783/1‑107 (MQ=255)
                           tAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTg                                                                                              >  2:478614/1‑139 (MQ=255)
                                       ctgcGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  1:428520/139‑1 (MQ=255)
                                       ctgcGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  1:1074619/139‑1 (MQ=255)
                                       ctgcGACATTTTCAAGCGCAGCCCCAAAAGGAAGCCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    <  2:1088236/121‑1 (MQ=255)
                                       ctgcGACATTTTCAAGCGCAGCCCCAAAAGGAAGCCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    >  1:1088236/1‑121 (MQ=255)
                                            aCATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCa                                                                                                           <  1:38608/109‑1 (MQ=255)
                                            aCATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCa                                                                                                           >  2:38608/1‑109 (MQ=255)
                                              aTTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  1:561395/1‑112 (MQ=255)
                                              aTTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:561395/112‑1 (MQ=255)
                                                 ttCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCa                                                                        >  1:742307/1‑139 (MQ=255)
                                                         cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATgg                                                                              >  2:1094282/1‑125 (MQ=255)
                                                         cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATgg                                                                              <  1:1094282/125‑1 (MQ=255)
                                                                aaaaGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCg                                                         <  2:429788/139‑1 (MQ=255)
                                                                      aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTc                                                   >  1:237185/1‑139 (MQ=255)
                                                                      aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTc                                                   >  1:1088939/1‑139 (MQ=255)
                                                                        gTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTg                                                 >  2:1049420/1‑139 (MQ=255)
                                                                           tGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcg                                              <  2:378039/139‑1 (MQ=255)
                                                                               aTGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcgcgcg                                          >  2:871897/1‑139 (MQ=255)
                                                                                  aGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCa                                       <  1:871897/139‑1 (MQ=255)
                                                                                  aGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCa                                       <  1:1049420/139‑1 (MQ=255)
                                                                                   gTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAg                                      >  1:977061/1‑139 (MQ=255)
                                                                                            aaGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATggg                             >  2:103231/1‑139 (MQ=255)
                                                                                                   cATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGcc                      >  2:874069/1‑139 (MQ=255)
                                                                                                     ttATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  2:389882/139‑1 (MQ=255)
                                                                                                     ttATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  2:330718/139‑1 (MQ=255)
                                                                                                              cTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGt                                                                                                       <  2:755415/47‑1 (MQ=255)
                                                                                                              cTGCCAACATTGGTGCCGATGTGATGATTAAAATTTTGCGTCACGGt                                                                                                       >  1:755415/1‑47 (MQ=255)
                                                                                                                      aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCGTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGc                       >  2:797832/1‑119 (MQ=255)
                                                                                                                      aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGc                       <  1:797832/119‑1 (MQ=255)
                                                                                                                      aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCGCCCATGAAGGGGTGa   <  2:742307/139‑1 (MQ=255)
                                                                                                                       ttGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAt  <  1:322992/139‑1 (MQ=255)
                                                                                                                       ttGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAt  <  2:977061/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                     
ATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372842‑373099

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