Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,072,253 C→T 13.8% T100T (ACC→ACT)  qedH‑II → quinoprotein ethanol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,072,2530CT13.8% 74.5 / 4.2 29T100T (ACC→ACT) qedH‑IIquinoprotein ethanol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (15/10);  new base T (2/2);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GCCCGCTGGACCTGCTCAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCA  >  NC_002947/3072123‑3072383
                                                                                                                                  |                                                                                                                                  
gCCCGCTGGACCTGCTCAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCAGATCAAGGACGGGGTGATGTACCTAAccccctc                                                                                                                               >  2:268264/1‑131 (MQ=255)
         aCCTGCTCAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTcgcg                                                                                                                      >  1:440873/1‑136 (MQ=255)
                cAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTc                                                                                                               >  1:179098/1‑136 (MQ=255)
                cAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTc                                                                                                               >  1:181756/1‑136 (MQ=255)
                      cAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTg                                                                                                         >  2:110091/1‑136 (MQ=255)
                           aCGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGc                                                                                                    >  2:515455/1‑136 (MQ=255)
                                aaGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATGAAGGACGGGGTGATGTACCTGACTGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCa                                                                                               >  2:312607/1‑136 (MQ=255)
                                aaGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACTGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCa                                                                                               >  2:636532/1‑136 (MQ=255)
                                    cgCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCgg                                                                                           <  1:398625/135‑1 (MQ=255)
                                    aGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCgg                                                                                           >  2:406782/1‑136 (MQ=255)
                                    aGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCgg                                                                                           >  2:12697/1‑136 (MQ=255)
                                          gcCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGaa                                                                                     >  1:630830/1‑136 (MQ=255)
                                            ccGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGcc                                                                                                            <  1:65354/111‑1 (MQ=255)
                                            ccGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGcc                                                                                                            >  2:65354/1‑111 (MQ=255)
                                                       gTTCTCCTTCGGCGGGGAGAAGCAGCGAGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGCGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGTCAATAc                                                                        >  1:5214/1‑136 (MQ=255)
                                                            ccTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGc                                                                   <  1:268264/136‑1 (MQ=255)
                                                                            gcagcGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGAc                                                   <  1:515455/136‑1 (MQ=255)
                                                                                   ggCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGcc                                            >  1:412616/1‑136 (MQ=255)
                                                                                         cagGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGTAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCctgct                                      <  2:131257/136‑1 (MQ=255)
                                                                                         cagGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCctgct                                      <  1:110091/136‑1 (MQ=255)
                                                                                               cAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGAcc                                >  2:78780/1‑135 (MQ=255)
                                                                                               cAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACg                                >  2:25037/1‑136 (MQ=255)
                                                                                                              gACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCgcggcg                 >  2:531815/1‑136 (MQ=255)
                                                                                                                cGGGGTGATGTACCTGACTGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTc               <  1:636532/136‑1 (MQ=255)
                                                                                                                cGGGGTGATGTACCTGACTGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTc               <  1:312607/136‑1 (MQ=255)
                                                                                                                 ggggTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTcg              <  2:630830/136‑1 (MQ=255)
                                                                                                                        tGTAGCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTa       <  1:271749/136‑1 (MQ=255)
                                                                                                                          tACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTACg     >  2:419478/1‑136 (MQ=255)
                                                                                                                             cTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCATGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCa  <  1:406782/136‑1 (MQ=255)
                                                                                                                             cTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCa  <  1:531815/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                  
GCCCGCTGGACCTGCTCAATGTCAATAACGTCAAGGAGCTGCGCCCGGTCTGGGCGTTCTCCTTCGGCGGGGAGAAGCAGCGCGGCCAGCAGGCCCAGCCGCTGATCAAGGACGGGGTGATGTACCTGACCGGCTCCTACTCGCGGGTGTTCGCCGTGGATGCCCGCACCGGCAAGAAACTGTGGCAATACGATGCACGTCTGCCGGATGACATCCGCCCCTGCTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCA  >  NC_002947/3072123‑3072383

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