Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 912,695 C→T P514P (CCG→CCA)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb912,6950CT80.0% 36.6 / 5.7 20P514P (CCG→CCA) hcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (16/0);  total (17/3)
Fisher's exact test for biased strand distribution p-value = 3.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGAATTTCTCGTTCAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAG  >  W3110S.gb/912637‑912742
                                                          |                                               
ggCCGAATTTCTCGTTCAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCggg                                              >  1:4163940/1‑62 (MQ=255)
    gAATTTCTCGTTCAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCgg                                          <  1:3663942/62‑1 (MQ=255)
    gAATTTCTCGTTCAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCgg                                          <  1:3643368/62‑1 (MQ=255)
                cAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATAtttt                              <  1:1297075/62‑1 (MQ=255)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:1665251/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:1755063/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:2311793/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:2639540/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:3453360/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:349185/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:1609845/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:1446765/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:36961/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:3789804/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:3801607/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:1170331/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:4410904/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:594528/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:722924/1‑51 (MQ=38)
                                                       aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag  >  1:923628/1‑51 (MQ=38)
                                                          |                                               
GGCCGAATTTCTCGTTCAGCACCGCCAGCAGGTCTGGTGTCAGGAAACCAGGTGCAGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAG  >  W3110S.gb/912637‑912742

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