Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,708,844 | 0 | C | G | 20.4% | 85.7 / 14.6 | 56 | A642A (GCC→GCG) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base C (27/16); new base G (5/6); total (34/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.24e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATT > NC_000913/1708719‑1708991 | cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAAACGCAATCGAATCCTCTATTGCAAGAGACAAAGAAACCAAGCTCGCCCAGAAAAATACTAATGCAGAGCCAGAAGAACAGGTGGATcc > 2:1493609/1‑149 (MQ=2) ggCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTAGCAGCGGCTATTGCCCGAGAGCGGGAGAGCGACCTGGAGCGGAAACATGGTAATGACGGTCCAGAACATCAGGTTGATCAGCGCAAAGCGGCCGtggccg > 2:1493708/1‑144 (MQ=2) aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg > 1:1493504/1‑104 (MQ=18) aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg < 2:1493504/104‑1 (MQ=18) cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTTCTGAAGAGCCACGAGAACGGGTAGATCCTAGCAATGCGGCCGTCGCGGCTGCTAtagaa > 2:1493630/1‑142 (MQ=11) ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCACCAGGCAAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGGCCAAGCCGCTATAGCCCGAGCCAAAgcgc > 2:3723309/1‑149 (MQ=14) gATCCTCGCAAAGCTGCCGTCGAAGCAGCTATAGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt < 1:1493517/99‑1 (MQ=14) gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa < 2:1493522/149‑1 (MQ=255) gATCCCCGCAAACCAGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa < 1:1493518/84‑1 (MQ=14) acccccgcAACCCCGCCGCCGAAGCCGCTATTGATTGTTCTAAAGCTCGCAAGCCGTAACAGAAACAGTCTATTTCCGAACAAGAAGATCTTGCCTATTCTCGCAAAGCTGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCaaa < 2:1493577/144‑1 (MQ=9) cccccccgTCGACGCCGCTACCGCCCCCGCCAAAGCTCTCAAGCTTTAACAGCATCAGCCTACTGCGTACACATCTCAACCGGTCGACGCGCGCAAAGCGTCCTTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc < 2:1493589/144‑1 (MQ=255) cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGCACAGCacc > 1:1493659/1‑147 (MQ=255) gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc < 1:1493530/86‑1 (MQ=255) gccgTCGAAGCAGCTATTGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa < 1:1493525/83‑1 (MQ=11) tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgc > 1:1493646/1‑91 (MQ=255) tattcGCTATTGCCCGTGCCAAAGCACGCAAGCTTGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt < 2:1493534/145‑1 (MQ=255) aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAGCAGGTCGATCCGCGCAAAGCAGCAGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAGCCGGCTAAt > 2:1493594/1‑149 (MQ=14) aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc > 1:1493658/1‑109 (MQ=255) gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTGGATCCGCGCAAAgc > 1:1493662/1‑83 (MQ=255) tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgc > 1:1493700/1‑81 (MQ=255) aTTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAGAACGAGGTCGAGACGAGGAAAGCAGCCGTCGAGGACTCTATTGCCCGCGGCAAAGAAAGAAAACGGGAACAGCAACCGGCTAGTGCGGGGcc > 2:1493665/1‑149 (MQ=9) gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGCAGCGGCTATTGCCCGCGCAAAAGCCCGAAAACCCGCACAGCAACAGGCTGATGCCGAGcc > 2:1493702/1‑146 (MQ=11) cccccgCCAAAGCCCGCAAGCTTGAACACCAACAGCCTAATGCCTACCCAGAAGAACAGCGCGACCCGCGCTAAGACGCGTGCGAATCCGATATCGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga < 1:1493599/143‑1 (MQ=16) cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTCGCCGCCGCTATTGCCCGCGCACAGGCCCGCAAACCCGCCCAGCAGCAGGCTAATGCCGAGCCcgaa > 2:1493668/1‑149 (MQ=11) aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAgcgg > 2:1493679/1‑68 (MQ=255) aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAg > 1:1493680/1‑67 (MQ=255) aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa > 1:1493544/1‑44 (MQ=255) aaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa < 2:1493544/44‑1 (MQ=255) accACGCAAGCTGGAACAGCAACAGGCTAATGCCGATCCATAACACCAGGGCGTTCCTCGCATAGCCTCAGTCTAAGCCTCTATTGCCCGATTCATAGCTCGCAATCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg < 2:1493545/147‑1 (MQ=14) gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATcc > 1:1493686/1‑52 (MQ=255) aGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg < 1:1493547/66‑1 (MQ=25) gCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTCGAAGCGGCTATTGCCCTCGACAAAGCGCGCAATCGTGAACAGCAACCGGCTAGTGAGGATACAGAAGAACAGGTTGATCCgcgcg > 2:1493692/1‑148 (MQ=14) tGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc < 2:1493550/102‑1 (MQ=255) tGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc > 1:1493550/1‑102 (MQ=12) tGGAACAGCAACAGGCTAATGCCGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt < 1:1493548/71‑1 (MQ=21) ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg < 1:1493552/37‑1 (MQ=255) ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg > 2:1493552/1‑37 (MQ=255) ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc > 1:1493613/1‑44 (MQ=255) ggAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgcgc > 1:1493651/1‑81 (MQ=255) aaCAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc < 1:1493554/102‑1 (MQ=255) aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc > 2:1493554/1‑102 (MQ=12) cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc > 1:1493698/1‑50 (MQ=255) cAACAGGCTAATGCCGAGCCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg < 1:1493556/67‑1 (MQ=255) aaCAGGCTAATTCCGAGCCAGAACAACAGGTCGATCcgc < 1:3655954/39‑1 (MQ=255) cAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGAGCGCAAACTGTAACAGAAAACTGATAATCAGGAGGAAGAACAAAAGCTTCATACACACAAAGCTGCAGTGGa > 2:1493706/1‑149 (MQ=9) aGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAc > 1:1493558/1‑96 (MQ=16) aGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCTATCGGGAACAGCAACCGGCTAATGCGGGGCCAGATGAACAGGTTGATCCGCGCACAGCTGCCGTCgca > 2:1493557/1‑147 (MQ=14) ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGTGCCAGAAGAACAGGTTGATc > 1:1493622/1‑128 (MQ=255) gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg < 2:1493559/85‑1 (MQ=11) gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg > 1:1493559/1‑85 (MQ=11) cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg < 2:1493558/96‑1 (MQ=255) tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa < 1:1493560/85‑1 (MQ=255) tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa > 2:1493560/1‑85 (MQ=255) tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc > 1:1493676/1‑39 (MQ=255) aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGTGCAAGAAGAACAGGTTGATCAGCGAAAAGCTGCAGTAGAAGCgcc > 2:1493736/1‑147 (MQ=14) gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCg > 1:1493726/1‑59 (MQ=255) ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt > 1:1493694/1‑149 (MQ=255) | CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATT > NC_000913/1708719‑1708991 |
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: ‑ |