Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | U00096 | 1,708,766 | 0 | C | T | 22.9% | 37.8 / 17.6 | 42 | V616V (GTC→GTT) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base C (13/14); new base T (4/4); total (23/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGC > U00096/1708623‑1708903 | cgcAAGCTGGAACACCAACAGGCTAATCCGGACCCAGCACCACAGGCCGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc < 2:1396545/149‑1 (MQ=14) tGGAACACCAACAGGAAAATTCGGAACCAGAACCAGTGGTCGATCCTCGAAAAGCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa < 2:1396548/149‑1 (MQ=14) gagcttcAGGCTAATCCGGACGCAGAACCCGTGGCGGATCCCCGTAAAGCGGCCGTGGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccgc < 2:1396564/143‑1 (MQ=11) caccaccAGGATAATGAGGAACCAGCACACGATGCTGATCCTCGTAAATCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccgc < 2:1396571/143‑1 (MQ=11) gCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGAGCGAACAGCGACAGCAACGGGCTGATGCGGAAACAGAAGAGACGGAGGATCCGAGAAAAGcagccg > 2:1396674/1‑149 (MQ=11) cccAGCACATCTGGCCGATCCTCGCAAAGCCGCCGTTGAAGCCGCTATTGCCCGTGCCCAAGCTCGCAAGCTGGAAGAGCAACAGGCTAATGCGGAAGCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGccc < 2:1396553/148‑1 (MQ=11) acccgaGGCCGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt < 2:1396554/121‑1 (MQ=14) aTCCGCGAAAAGCCGCCGTCGAATCCGCTATAGCCCGGGCCTAAGCGCTCAAGCTGGAACAGCTACAGGCTAATGCGGAACCAGAAGAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg < 1:1396561/149‑1 (MQ=14) tCCTCGAAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCTAATCTCGCAATCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGc < 2:1396562/149‑1 (MQ=11) aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTGGATCCGCGCAAAg > 1:1396645/1‑99 (MQ=255) gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAATCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaaa < 2:1396652/90‑1 (MQ=11) gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaaa > 1:1396652/1‑90 (MQ=255) gAAGCTGCTATAGCCAGAGCCATAGCGCGCAATCTTTAACAGCATCAGTCTAATTCGGTACTAGTAGAAGATGTCGATCCTCGCAAATCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTaa < 2:1396589/149‑1 (MQ=11) aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCa > 2:1396658/1‑84 (MQ=255) gcgagtGCCCGGGCCACCGCGGGCAAGCTGGAACAGCAACAGGCTACTCCGGAAGCAGCACCCCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa < 1:1396586/144‑1 (MQ=11) gCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAgcagcc > 1:1396660/1‑87 (MQ=11) cTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcg < 2:1396725/74‑1 (MQ=255) cTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcg > 1:1396725/1‑74 (MQ=255) ttGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATc < 1:1396730/69‑1 (MQ=255) ttGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATc > 2:1396730/1‑69 (MQ=255) tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCATCAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg < 2:1396580/149‑1 (MQ=255) ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAAGAACAGGTGGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAt < 2:1396665/88‑1 (MQ=11) ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAt > 1:1396665/1‑88 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAgctgcc < 2:1396755/67‑1 (MQ=11) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAgctgcc > 1:1396755/1‑67 (MQ=255) gcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt < 1:1396583/79‑1 (MQ=34) ggAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGATCCGCGCAAAgc > 1:1396765/1‑54 (MQ=255) aaCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc < 1:1396587/148‑1 (MQ=21) gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgga > 1:1396677/1‑104 (MQ=25) gCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAgg > 1:1396675/1‑104 (MQ=25) aaCAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa > 1:1396683/1‑44 (MQ=38) aCAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGc > 2:1396678/1‑86 (MQ=11) aGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGc > 1:1396682/1‑64 (MQ=255) ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTa > 1:1396681/1‑101 (MQ=35) gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCACAAAGCAGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGCt > 2:1396679/1‑83 (MQ=11) cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc > 1:1396726/1‑40 (MQ=40) tAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgc > 2:1396779/1‑39 (MQ=255) cAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTa > 1:1396687/1‑45 (MQ=255) gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt > 1:1396689/1‑107 (MQ=12) gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGATGGAACAGCAACAGGCTAATGCCGAGAAAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCTGatg > 1:3111119/1‑147 (MQ=255) agaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTaa > 1:1396690/1‑84 (MQ=37) gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc < 1:1396588/103‑1 (MQ=255) | CGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGC > U00096/1708623‑1708903 |
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: ‑ |