Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,004,302 G→A 100% S483F (TCT→TTT)  hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,004,3020GA95.9% 72.7 / ‑3.8 24S483F (TCT→TTT) hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (0/23);  total (1/23)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCG  >  minE/2004238‑2004353
                                                                |                                                   
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1050206/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:452573/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:501040/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:54612/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:36546/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2277618/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2231393/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2118754/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2038143/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2014924/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1891603/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1861324/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1788786/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:716486/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:739087/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1708773/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1359767/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1288814/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1286299/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:829500/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1275880/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1182671/67‑1 (MQ=255)
ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:943265/67‑1 (MQ=255)
 cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:767102/66‑1 (MQ=255)
 cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:961812/66‑1 (MQ=255)
 cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:606708/66‑1 (MQ=255)
     tttaTGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1212929/59‑1 (MQ=255)
                     cTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:2285193/46‑1 (MQ=255)
                        gTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaag                                                   <  1:1501174/43‑1 (MQ=255)
                                                 ataaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCg  >  1:481311/1‑67 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:913081/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:563701/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:573687/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:664355/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:104841/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:966843/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:71328/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:464150/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:391018/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:369247/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:241622/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1979603/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1864782/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:180036/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1775085/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1753465/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1725295/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1281199/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:124463/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1218999/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1164724/52‑1 (MQ=255)
                                                    aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:1143473/52‑1 (MQ=255)
                                                     ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCtt              <  1:2060463/51‑1 (MQ=255)
                                                                |                                                   
CCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCG  >  minE/2004238‑2004353

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