Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,071,959 A→C 8.6% H271P (CAC→CCC)  yiaY → Fe‑containing alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,071,9590AC 8.6% 127.7 / 3.4 46H271P (CAC→CCC) yiaYFe‑containing alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base C (2/2);  total (24/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 = 9.36e-01

GCCTGATCAATGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCCGCGGTCTCA  >  NC_002947/3071821‑3072102
                                                                                                                                          |                                                                                                                                               
gcCTGATCAATGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGggc                                                                                                                                       <  5:216229/149‑1 (MQ=255)
      tCAATGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGCCGCCCAGGCTGtgcgacatg                                                                                                                                 >  3:92775/1‑144 (MQ=255)
      tCAATGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGCCGCCCAGGCTGtgcgacatg                                                                                                                                 >  3:88856/1‑144 (MQ=255)
        aaTGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGCCGCCCAGGCTGtgcgacatggc                                                                                                                               <  4:92775/149‑8 (MQ=255)
        aaTGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGCCGCCCAGGCTGtgcgacatggc                                                                                                                               <  4:88856/149‑8 (MQ=255)
         aTGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTCCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTc                                                                                                                              <  8:59951/149‑1 (MQ=255)
            gCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGAc                                                                                                                           >  2:230721/1‑149 (MQ=255)
            gCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGAc                                                                                                                           <  3:228495/149‑1 (MQ=255)
             cAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACt                                                                                                                          >  3:223771/1‑149 (MQ=255)
              aaCCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACtt                                                                                                                         >  3:75128/1‑149 (MQ=255)
                        agaTGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGc                                                                                                               >  8:153373/1‑149 (MQ=255)
                         gatgatCGCCAACGCTACTGATATTGGGCTGCGCGAGAAGATCATGCTGGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGc                                                                                                               <  4:75128/148‑1 (MQ=255)
                           tgatCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTtg                                                                                                            <  1:76786/149‑1 (MQ=255)
                           ggatCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGCTCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTtg                                                                                                            <  4:42929/148‑1 (MQ=255)
                            gatCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTtgt                                                                                                           >  3:118520/1‑149 (MQ=255)
                               cGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGc                                                                                                         >  7:185382/1‑148 (MQ=255)
                                      ccTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCtt                                                                                                             <  1:164672/137‑1 (MQ=255)
                                      ccTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCtt                                                                                                             >  2:164672/1‑137 (MQ=255)
                                          cTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCtg                                                                                             >  3:182557/1‑149 (MQ=255)
                                              tatTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCtggtg                                                                                          >  3:149968/1‑148 (MQ=255)
                                                ttGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGa                                                                                        <  1:33230/148‑1 (MQ=255)
                                                ttGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAg                                                                                       >  5:154530/1‑149 (MQ=255)
                                                                 aTCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCaa                                                                      <  4:118520/149‑1 (MQ=255)
                                                                 aTCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTAGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCaa                                                                      >  5:1574/1‑149 (MQ=255)
                                                                  tCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCaa                                                                      <  6:154530/148‑1 (MQ=255)
                                                                  tCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAAc                                                                     <  1:49310/149‑1 (MQ=255)
                                                                            cagcaTGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCa                                                                                     >  6:254090/1‑123 (MQ=255)
                                                                               caTGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTGGCACAGGCTGGGCGGCTTCCTCGAGTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGgcgc                                                        <  8:185382/149‑1 (MQ=255)
                                                                                  gCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTc                                                                        >  4:242083/1‑130 (MQ=255)
                                                                                  gCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTc                                                                        <  3:242083/130‑1 (MQ=255)
                                                                                     ggCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGc                                                  >  5:135734/1‑149 (MQ=255)
                                                                                             tGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCGCGCCATGTCGCACAGCCTGGGCGGCTTCCTc                                                                                                                              <  5:21529/65‑1 (MQ=255)
                                                                                             tGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCGCGCCATGTCGCACAGCCTGGGCGGCTTCCTc                                                                                                                              >  6:21529/1‑65 (MQ=255)
                                                                                               gCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAgg                                         <  4:149968/148‑1 (MQ=255)
                                                                                               gCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGt                                        >  7:254422/1‑149 (MQ=255)
                                                                                                     tcCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGc                                  >  5:103600/1‑149 (MQ=255)
                                                                                                     tcCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGc                                  >  7:20805/1‑149 (MQ=255)
                                                                                                      cgaaTGCGATCCTGGGCGCAGTGCAGGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAg                                          <  8:91434/138‑1 (MQ=255)
                                                                                                      ccAATGCGATCCAGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAg                                          >  7:91434/1‑140 (MQ=255)
                                                                                                       cAATGCGATCCTGGGCTTTTTTAAAGCCATGTCGCACAGCCTGGGAGGCTTCCTAGACTTGCCCCATTGCTTGTGAAAAGGGGTGCTGGTGGAGCACGTTTTTGCGTTGAACTACAGCTGGGCGCCGGAGCGTTTCATGGTGATCggg                                 <  4:182557/148‑3 (MQ=255)
                                                                                                        aaTGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGa                               >  8:130780/1‑149 (MQ=255)
                                                                                                              aTCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGtt                         <  4:223771/149‑1 (MQ=255)
                                                                                                                              aCGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCcgc         >  5:62387/1‑149 (MQ=255)
                                                                                                                                      tCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCCGCGGTCTCa  <  8:254422/148‑1 (MQ=255)
                                                                                                                                        gCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCCGCGGt      >  2:239293/1‑142 (MQ=255)
                                                                                                                                        gCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCCGCGGt      <  1:239293/142‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GCCTGATCAATGGCAACCTGGTGGAGATGATCGCCAACCCTACTGATATTGCGCTGCGCGAGAAGATCATGCTCGGCAGCATGCAGGCGGGGCTGGCGTTCTCCAATGCGATCCTGGGCGCAGTGCACGCCATGTCGCACAGCCTGGGCGGCTTCCTCGACTTGCCCCATGGCTTGTGCAACGCGGTGCTGGTGGAGCACGTGGTGGCGTTCAACTACAGCTCGGCGCCGGAGCGTTTCAAGGTGATCGCCGAGGTGTTTGGTATCGACTGCCGCGGTCTCA  >  NC_002947/3071821‑3072102

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