Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | U00096 | 1,708,751 | 0 | A | G | 32.3% | 18.9 / 20.7 | 32 | E611E (GAA→GAG) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base A (9/12); new base G (6/4); total (16/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTG > U00096/1708611‑1708896 | cgcgCCTAACCGCGGAAGCTGGAACAGCTACAGGCTAATCCGGACCCAGAACAACAGGCCGATCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga < 2:1197234/149‑1 (MQ=14) ctaaGCGCGCAATCTGGAAGAGCATCAGGCTAATTCGGAACCAGCACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg < 2:1197238/147‑1 (MQ=14) aGCTGGAACACCATCAGCCTACTTCGGACCCAGAAACACATGCCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCgcgc < 2:1197249/149‑1 (MQ=14) cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc < 2:1197242/149‑1 (MQ=255) gAAGCAGAACCACAGGTCGATCCTCGAAAATCAGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc < 1:1197244/149‑1 (MQ=14) ggccGACCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAc < 2:1197245/146‑1 (MQ=255) gATCCGCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc < 1:1197255/96‑1 (MQ=255) gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa < 1:1197256/92‑1 (MQ=255) aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCaga > 1:1197332/1‑77 (MQ=255) aaaCCCGCCGTCGAAGCAGCTATCGCCCGCGCCATAGCTCGCAATCTGGAACAGCAACAGGCTAATGCGGAAGCAGAAGACGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa < 2:1197257/149‑1 (MQ=11) gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga < 1:1197261/71‑1 (MQ=255) cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt > 2:1197338/1‑75 (MQ=255) aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt > 1:1197341/1‑73 (MQ=255) aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGACGCTATTGCCAGTGCCAAAGCACGCAAGCAGGAACAGCAACAGGCTAAt > 2:2386926/1‑149 (MQ=14) gCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAgg > 1:1197345/1‑67 (MQ=255) aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc < 2:1197416/71‑1 (MQ=255) aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc > 1:1197416/1‑71 (MQ=255) ttGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGCTCGATCCGCGCAAGGACGCCGTCGCAGCCGCTATTGCCCGTGCCAAAGCCCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCa > 2:1197308/1‑149 (MQ=11) gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt > 1:1197352/1‑50 (MQ=255) gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa < 2:1197264/62‑1 (MQ=255) cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGtt < 2:1197352/50‑2 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAg > 1:1197310/1‑46 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAg < 2:1197310/46‑1 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCgcgc < 1:1197425/58‑1 (MQ=11) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCGGCTATTGCACGAGCCAAAGCACGCAAGCTg > 1:1197353/1‑106 (MQ=11) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCgcgc > 2:1197425/1‑58 (MQ=255) cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg < 1:1197265/46‑1 (MQ=38) aGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTAtt > 1:1197357/1‑77 (MQ=21) tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgctgc > 1:1197355/1‑58 (MQ=11) gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGca > 2:1197361/1‑69 (MQ=255) ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCaa > 1:1197364/1‑85 (MQ=255) ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGGGACAGATCACGCACGGGGTAGCAGAGACAGGATAATGCGGAGCGAGAAGAAGAGGTAGATCAGCGCAAGGCGGCCGTCGAAGCGGCTATTg > 2:1197365/1‑149 (MQ=2) | CGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTG > U00096/1708611‑1708896 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |