Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,008,326 G→A 100% C232C (TGC→TGT hybO ← hydrogenase 2, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,008,3260GA81.7% 54.8 / 5.4 38C232C (TGC→TGThybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base A (15/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGACGGCCATAGGCGAAGGTCGGACGGTTTTTGTC  >  minE/2008273‑2008376
                                                     |                                                  
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2200853/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:3418798/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:356216/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:447611/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2955458/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2618627/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2584995/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2287363/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2057804/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2053902/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:601788/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:2005836/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:1783154/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:669552/69‑1 (MQ=255)
tcatcaCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:1779154/69‑1 (MQ=255)
    caCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGAc                                 >  1:214639/1‑69 (MQ=255)
    caCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGAc                                 >  1:2220673/1‑69 (MQ=255)
    caCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGAc                                 >  1:611397/1‑69 (MQ=255)
    caCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGAc                                 >  1:1635985/1‑69 (MQ=255)
     aCCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGACg                                >  1:2707066/1‑69 (MQ=255)
         aaCTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGACGGCCa                            <  1:454048/69‑1 (MQ=255)
           ctctTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGACGGCCAt                           >  1:530328/1‑68 (MQ=255)
                     aaaCGACCAGCATCGAAGTGCGGGCGACGTTCACAGTGTTCGTGAATc                                     <  1:1170212/48‑1 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGt   >  1:2509324/1‑68 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:474994/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:514692/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:1022966/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:844139/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:3248583/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:3182209/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:2956789/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:2525967/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:249643/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:2113292/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:1918635/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTCGGACGGTTTTTGTc  >  1:1262780/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGTGTTCGTGAATCAGACGACCATAGGCGAAAGTAGGACGGTTTTTGTc  >  1:14325/1‑69 (MQ=255)
                                   gAAGTGCGGGCGACGTTCACAGGGTTCGTGAATCAGACGACCATAGGCGAAAGTCgg              >  1:1378699/1‑57 (MQ=255)
                                                     |                                                  
TCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTTCGCAGTGTTCGTGAATCAGACGGCCATAGGCGAAGGTCGGACGGTTTTTGTC  >  minE/2008273‑2008376

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