Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 31.4% | 67.8 / 10.2 | 35 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (18/6); new base C (0/11); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.97e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGC > REL606/2631697‑2631846 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca > 1:1089897/1‑76 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac < 2:1087052/76‑1 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca > 1:1422413/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca > 1:516088/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCATccacca > 1:1688981/1‑76 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 2:895815/76‑1 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 2:1004730/76‑1 (MQ=255) gtgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagc > 1:911237/1‑76 (MQ=255) gtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcca < 2:1059871/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 2:1658595/76‑1 (MQ=255) gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtc > 2:714111/1‑76 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCtt < 2:20821/76‑1 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 1:1613671/1‑76 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 1:33463/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 2:12672/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 2:1200646/1‑76 (MQ=255) cTGCCTGAAGAGTACCCCAACGCCCTTTTTATTCCCCCCCATAACCCTCAACCAGCCAGTAACTTCTCTTtctctc < 1:1295494/76‑1 (MQ=255) cTGCCTGAAGAGTAACCCAACGGCTTTTTTATTTCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTtctc < 2:1284006/76‑1 (MQ=255) cTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtctc > 1:459351/1‑76 (MQ=255) tGCCTGAAGAGTACCCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcg < 1:1055619/76‑1 (MQ=255) tGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcg > 2:874388/1‑76 (MQ=255) cgCCTGAAGAGTACCCCAACGGCCTTTTTATTTCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTcg < 2:1301495/75‑1 (MQ=255) cTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccg > 1:330160/1‑76 (MQ=255) gTACCCCACCGGCCTTTTTATTTCACCACCTAATCCACCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg < 2:1402679/76‑1 (MQ=255) tAACCAAACGGCCTTTTTATTTCACCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt < 2:500415/76‑1 (MQ=255) accccAACGCCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 1:767891/74‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 1:15487/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:1526972/1‑76 (MQ=255) ccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc > 2:133485/1‑76 (MQ=255) aCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtg < 2:1630175/76‑1 (MQ=255) gCCTTTTTATTTCCCCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTNGCCGCCCTGCGTCAGCGTGttt < 2:232300/76‑1 (MQ=255) ccttttttATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:1085055/73‑1 (MQ=255) tttATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa < 1:1670089/76‑1 (MQ=255) tATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt > 1:916681/1‑76 (MQ=255) tttCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt < 1:93806/76‑1 (MQ=255) ccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata > 2:1468768/1‑76 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 2:1052551/1‑76 (MQ=255) cTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat < 1:1548912/76‑1 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 2:611964/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1371505/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1244291/1‑76 (MQ=255) caccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGc > 1:831813/1‑76 (MQ=255) caccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGc > 1:348946/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGC > REL606/2631697‑2631846 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |