Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb912,7190GA61.5% 13.6 / 25.3 26G506G (GGC→GGT) hcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (9/1);  new base A (16/0);  total (25/1)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCAA  >  W3110S.gb/912692‑912765
                           |                                              
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:349185/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:923628/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:722924/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:594528/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:4410904/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:3801607/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:3789804/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:36961/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:1170331/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:3453360/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:2639540/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:2311793/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:1755063/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:1665251/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:1609845/1‑51 (MQ=38)
aGTTGGGCCGGTGACGATATTTTTCACACCCAGAGAAAGCAGCGTcagcag                         >  1:1446765/1‑51 (MQ=38)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:3080788/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:2522901/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:2446506/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:3773145/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:2293733/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:4455991/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:4494445/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:4617588/1‑62 (MQ=255)
           tGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCa   >  1:1532968/1‑62 (MQ=255)
            gACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCaa  <  1:4573093/62‑1 (MQ=255)
                           |                                              
AGTCGGGCCGGTGACGATATTTTTCACGCCCAGAGAAAGCAGCGTCAGCAGAATGACGATCGCTTTCTGTTCAA  >  W3110S.gb/912692‑912765

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: ‑