Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,708,751 | 0 | A | G | 32.4% | 34.8 / 21.8 | 37 | E611E (GAA→GAG) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base A (5/18); new base G (9/2); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCC > NC_000913/1708609‑1708867 | cccTCGCCAAAGCGAGCAAGCTGGAACAGCATCAGGCTAATTCGGAAGCAGAACAACAGGTCGATCCTCGCAAGTCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaa < 2:1050323/149‑1 (MQ=14) atgctcgcAAGCTGGAAGAGCAACAGGATAATGAGGAACCAGAACACCTGGTCGATCCTCGTAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc < 2:1050328/144‑1 (MQ=14) aacccgcAAGCAGGAACATCAACAGGAAAATTAGGAACCAGCACAAGATGTCGATCCTCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCg < 2:1050331/145‑1 (MQ=14) cgcAATCTGGAAGATCATCAGGATAATTCGGAACCAGAAAAACATGCCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa < 2:1050325/146‑1 (MQ=14) cAAGCTGGAACAGCATCAGGCTAATTAGGAACCAGAACACCATGCCGATCCGCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc < 1:1050330/149‑1 (MQ=14) aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCGAGCAAGATGGAACCGCAGCAGGAGACTGAGGAGCCAGAAGCAAAGATAGATCCGCGCAAAgcggcc > 2:1050391/1‑149 (MQ=11) gccgAACCAGAACCCCTGGCTGATCCTCGAAAAACAGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt < 2:1050344/146‑1 (MQ=11) ggAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCaga > 2:1050393/1‑105 (MQ=255) aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGAGCGCAACCTGGAACAGAAACAGGCAAATGCCGATCCAGAAGATCAGGTAGATCAGCTCAAAGCGGCCGTCGAAGCGGCTATTGCCCGTg > 2:1050444/1‑149 (MQ=11) acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa < 2:1050332/149‑1 (MQ=255) aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGAAACAGGCTAATGCGGTACCAGCAGAACAGCTCGATCCGCGCAAAGCCGCCGTCGCAGCAGCGATTGCCCGTGCCaaa > 1:2600368/1‑149 (MQ=14) acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACa > 1:1050446/1‑99 (MQ=255) gATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTGATGCGGAACCAGAAGAACa < 1:1050338/92‑1 (MQ=255) gATCCGCGCAAAGCAGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc < 1:1050336/103‑1 (MQ=255) gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGcc > 1:1050410/1‑76 (MQ=255) gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga > 1:1050409/1‑79 (MQ=21) aGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTAGACAGAAGAACCGGGGGAGCCGCGCAAAGCCGCCGTCTCAGCAGCTCTTTCCAGTGAGATAGCACTCAAGATCGAAAAGc > 1:1286488/1‑149 (MQ=2) gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc < 1:2881422/67‑1 (MQ=255) gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt < 1:1050346/84‑1 (MQ=255) gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt < 1:1050347/84‑1 (MQ=255) gccgTCGAAGCCGCTATCGCCCGTGACAAAGCGCGCAAGCTGGAACAGCAACATTCTAATGCGGAACCCCATCAACAGTTCGTTCCTCGTAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACa < 2:1050348/149‑1 (MQ=11) gccgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa < 1:1050345/83‑1 (MQ=255) aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga > 1:1050418/1‑64 (MQ=255) gCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc < 1:1050351/84‑1 (MQ=18) gCTATAGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg < 1:1050352/88‑1 (MQ=255) cccGTGCCAAAGCGCTCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATc < 2:1050426/66‑1 (MQ=14) cccGTGCCAAAGCGCACAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATc > 1:1050426/1‑66 (MQ=14) ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAg > 1:1050429/1‑58 (MQ=255) cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa < 1:1050428/83‑1 (MQ=11) cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTGCCGTCGa > 2:1050428/1‑83 (MQ=11) cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc < 1:1050354/97‑1 (MQ=18) cgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAAc < 1:1050355/117‑1 (MQ=35) cgcAAGCTGGAACAGCAACATGATAATGATTATCCTGAACAAGAGTTCTATCCTCGTAAAGACGACGTCGATGCAGCTATTTCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc < 2:1050376/149‑1 (MQ=11) cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgctgc > 1:1050498/1‑65 (MQ=25) aGCAACAGGCTAATGCGGATCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTc < 1:1050436/56‑1 (MQ=11) tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGCCGCGGCTATTGCCCg > 1:1050514/1‑63 (MQ=11) gAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGaa > 1:1050442/1‑42 (MQ=255) | CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCC > NC_000913/1708609‑1708867 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |