Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,072,445 G→A 19.0% K164K (AAG→AAA qedH‑II → quinoprotein ethanol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,072,4450GA19.0% 42.0 / 5.0 21K164K (AAG→AAAqedH‑IIquinoprotein ethanol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base A (2/2);  total (12/9)
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.71e-01

CTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATCCAGGC  >  NC_002947/3072346‑3072576
                                                                                                   |                                                                                                                                   
ctgcGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAaagaa                                                                                                 <  2:440873/136‑1 (MQ=255)
 tgcGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAaagaag                                                                                                >  2:410241/1‑136 (MQ=255)
 tgcGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAaagaag                                                                                                >  2:177023/1‑136 (MQ=255)
  gcGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAgg                                                                                               >  1:574169/1‑136 (MQ=255)
                                  ggCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAAGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCa                                                               >  2:531339/1‑136 (MQ=255)
                                  ggCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAAGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCa                                                               >  2:528019/1‑136 (MQ=255)
                                          ggTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAg                                                       >  2:379432/1‑136 (MQ=255)
                                            tgtTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACcac                                                                                                       >  1:347101/1‑86 (MQ=255)
                                            tgtTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACcac                                                                                                       <  2:347101/86‑1 (MQ=255)
                                               tcttcGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGAt                                                  >  2:65844/1‑136 (MQ=255)
                                                                            gCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCg                     >  2:166938/1‑136 (MQ=255)
                                                                                tGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCgtggt                 <  1:166938/136‑1 (MQ=255)
                                                                                 gAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCgtggtg                <  2:574169/136‑1 (MQ=255)
                                                                                  aacaacGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGtggtgg               <  1:25037/136‑1 (MQ=255)
                                                                                  aaCAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGtggtgg               <  1:78780/136‑1 (MQ=255)
                                                                                  aaCAAGGACACCGGCAAAGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGtggtgg               <  1:531339/136‑1 (MQ=255)
                                                                                  aaCAAGGACACCGGCAAAGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGtggtgg               <  1:528019/136‑1 (MQ=255)
                                                                                          cacCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATc       >  1:500794/1‑136 (MQ=255)
                                                                                               gCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATCCAGGc  <  1:626774/136‑1 (MQ=255)
                                                                                               gCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATCCAGGc  >  2:57704/1‑136 (MQ=255)
                                                                                               gCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATCCAGGc  >  2:626774/1‑136 (MQ=255)
                                                                                                   |                                                                                                                                   
CTGCGACGTAATCAACCGCGGCGTCGCGCTGTACGGCAACCTGGTGTTCTTCGGCACGCTGGACGCCAAGCTGGTGGCCCTGAACAAGGACACCGGCAAGGTGGTCTGGAGCAAGAAGGTCGCCGACCACAAAGAAGGCTACTCCATCAGCGCCGCGCCGATGATCGTCAATGGCAAGCTGATCACAGGCGTTGCCGGCGGCGAGTTCGGCGTGGTGGGCAAGATCCAGGC  >  NC_002947/3072346‑3072576

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: