Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 343,467 A→C H195P (CAC→CCC)  yahK → broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,4670AC95.5% 151.4 / ‑3.7 44H195P (CAC→CCC) yahKbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (21/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACG  >  NC_000913/343373‑343560
                                                                                              |                                                                                             
tATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGa                                                                                          >  2:301365/1‑100 (MQ=255)
 aTTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGAt                                                                                         <  2:208638/100‑1 (MQ=255)
       cgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGc                                                                                   <  2:808005/100‑1 (MQ=255)
            cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACAtg                                                                              >  1:553565/1‑100 (MQ=255)
             gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACAtgt                                                                             >  2:456918/1‑100 (MQ=255)
               cACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGCGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGTGGCGCATgtgg                                                                           <  2:827094/100‑1 (MQ=255)
                  tGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGtggtgg                                                                        <  1:118904/100‑1 (MQ=255)
                    gCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCa                                                                      <  2:274017/100‑1 (MQ=255)
                     cAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCAt                                                                     >  1:1119666/1‑100 (MQ=255)
                     cAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCAt                                                                     >  2:137985/1‑100 (MQ=255)
                        gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTa                                                                  <  2:339613/100‑1 (MQ=255)
                          cgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTAcc                                                                <  1:389606/100‑1 (MQ=255)
                                       aaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACACGTGGTGGCATTTACCACTTCTGAGGCaa                                                   >  2:1136115/1‑100 (MQ=255)
                                        aaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCaaa                                                  <  2:1539392/100‑1 (MQ=255)
                                           tGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTACGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAcg                                               <  1:292875/100‑1 (MQ=255)
                                              gCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGa                                            <  1:456920/100‑1 (MQ=255)
                                              gCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGa                                            <  2:590377/100‑1 (MQ=255)
                                                 tGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGc                                         >  1:840288/1‑100 (MQ=255)
                                                   gTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCgg                                       >  1:1004170/1‑100 (MQ=255)
                                                     cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCa                                     >  1:1371394/1‑100 (MQ=255)
                                                      ggCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaa                                    >  1:1490692/1‑100 (MQ=255)
                                                       gCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaaa                                   >  2:607002/1‑100 (MQ=255)
                                                        cATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaaaa                                  >  1:1125152/1‑100 (MQ=255)
                                                        cATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaaaa                                  >  2:679033/1‑100 (MQ=255)
                                                          tCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGc                                >  2:365084/1‑100 (MQ=255)
                                                               ggTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTgg                           <  2:477262/100‑1 (MQ=255)
                                                                  cTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAACAGCCCTggggg                        >  2:713943/1‑100 (MQ=255)
                                                                  cTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTggggg                        >  1:494812/1‑100 (MQ=255)
                                                                    gggACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                      <  1:454227/100‑1 (MQ=255)
                                                                    gggACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                      <  1:454215/100‑1 (MQ=255)
                                                                    gggACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                      <  1:254849/100‑1 (MQ=255)
                                                                        cATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                  <  2:847385/100‑1 (MQ=255)
                                                                        cATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                  <  2:893714/100‑1 (MQ=255)
                                                                           atGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAg               <  1:607005/100‑1 (MQ=255)
                                                                           atGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAg               <  1:1339752/100‑1 (MQ=255)
                                                                            tGGGGATTAAGCTGGCTCACGGGGTGGGGGGACATGTGGTGGCGTTTACAACTTCTGAGGCAAAACACGAAGCGGGAAAAGGCCGGGGGGCCGAtcgag               >  2:151235/1‑95 (MQ=255)
                                                                            tGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgt              >  2:479746/1‑100 (MQ=255)
                                                                               ggATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgt           >  1:1274127/1‑100 (MQ=255)
                                                                                gATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt          <  2:1180096/100‑1 (MQ=255)
                                                                                gATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgt           <  2:1175213/99‑1 (MQ=255)
                                                                                 aTTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTa         >  1:1320032/1‑100 (MQ=255)
                                                                                     aGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTc     >  1:420690/1‑100 (MQ=255)
                                                                                       cTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCAc   <  2:175104/100‑1 (MQ=255)
                                                                                        tGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACg  >  2:329593/1‑100 (MQ=255)
                                                                                              |                                                                                             
TATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACG  >  NC_000913/343373‑343560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: