Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | U00096 | 1,708,751 | 0 | A | G | 21.1% | 56.7 / 10.1 | 39 | E611E (GAA→GAG) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base A (10/20); new base G (5/3); total (15/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATT > U00096/1708609‑1708895 | acctcgcgAAAGCGCTCAAGCAGGAACATCATCAGGCTACTCAGGACGCAGCACCAGAGGCCGATCCTCGTAAATCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaa < 2:1427741/145‑1 (MQ=9) ctgcAAAGCGCGCAAGCTGGAACAGCTACAGTCTAATGCGGAAACAGAACAACATGTCTATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc < 1:2622052/146‑1 (MQ=14) aaGCTGGCAGACCATCAGGCAAATGACGACCCAGCACCCCTGGGCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTTCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg < 1:1427728/149‑1 (MQ=14) cTGGAACAGCATCAGGATAATCCGGACGCAGAACACCAGGGCGATCCGCGAAAAGCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaa < 1:1427715/149‑1 (MQ=14) aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg < 1:1427717/149‑1 (MQ=255) aCAGGCTAATTCGGTACCACCACATCTTGCCGATCCTCGCAAATCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCg < 2:1427729/149‑1 (MQ=14) cAGGCTAATCAGGACCCAGCACCACATGCCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGa < 2:1427731/149‑1 (MQ=14) atacttaggAAGCAGAAAAAGTTGTCGATCCTCGCAAACCCGCCGTTGAAGCAGCTTTAGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc < 2:1427725/142‑1 (MQ=11) aaCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc > 2:2010044/1‑149 (MQ=255) aCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGccc < 2:1427723/149‑1 (MQ=255) accaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa < 2:1427724/147‑1 (MQ=255) aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa > 2:1427874/1‑89 (MQ=255) gATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg < 1:1427735/93‑1 (MQ=255) gATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg < 1:1427734/94‑1 (MQ=255) gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa < 1:1427737/92‑1 (MQ=255) ccccaacccgccgccgAAGCCCCTATCGCCCGCGTCAAAGCTATCTAGATGGTACAGCAACCTGCTAAACAGTAAACAGAACCACCTCTCGATCCGCTCATATCTGCCGTTGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGa < 2:1427772/142‑1 (MQ=2) ccaccccgccgTCGAAGCCGCTATCGCCCGCGCCAAATTTCTAAAGCTTTAAGTGGAACAGGGTTATGTTGAACCATAACAACATGGCGTCGCTTGCTAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc < 2:1427768/144‑1 (MQ=2) accgccgccCTCGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCTTTAACTGCAACTGGCTAATGCGGAACCAGAAGAACAGGTCTATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa < 2:1427738/146‑1 (MQ=11) aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa < 1:2693607/149‑1 (MQ=255) cccccccgccgAAGCCGCTCTCGCCCGCGCCAACGCGCTCTAGCTGTAACTGCAACGGGCTAATGCGGAACCACAACACCATGCCGATCCGCGCATCTCCGCAGCCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGc < 2:1427744/144‑1 (MQ=9) cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCagc > 2:1427828/1‑147 (MQ=11) aaGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt < 1:1427746/73‑1 (MQ=255) gCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga > 1:1427821/1‑62 (MQ=255) cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc > 1:1427747/1‑74 (MQ=255) cTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAAc > 1:1427907/1‑63 (MQ=255) cTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAAc < 2:1427907/63‑1 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt < 1:1427842/48‑1 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt > 2:1427842/1‑48 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc > 1:1427933/1‑54 (MQ=11) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc < 2:1427933/54‑1 (MQ=11) gcAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTc < 2:1427849/68‑1 (MQ=11) gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTc > 1:1427849/1‑68 (MQ=25) gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTc > 1:1427850/1‑68 (MQ=25) aaGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgc > 2:1427843/1‑59 (MQ=255) aGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAgctgc > 2:1427846/1‑61 (MQ=11) aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgc > 1:1427851/1‑52 (MQ=255) gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc > 1:1427801/1‑69 (MQ=35) tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAACCGCGCAATGCCGCCTTccgtgacg > 2:1427752/1‑142 (MQ=14) cGGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAtt < 2:1427755/149‑1 (MQ=255) | CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATT > U00096/1708609‑1708895 |
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: ‑ |