Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 35.2% | 108.7 / 17.6 | 54 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (29/6); new base C (0/19); total (29/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.64e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631698‑2631845 | aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac > 1:167154/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca > 1:985726/1‑76 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag > 2:1379001/1‑76 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 2:1526112/76‑1 (MQ=255) gtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcca > 1:1204368/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 1:765315/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 1:1273603/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 2:712086/76‑1 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACAAGCCAGaaa > 1:393497/1‑76 (MQ=255) aaaGGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa < 2:1672415/76‑1 (MQ=255) tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACtt > 1:1476830/1‑76 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 1:1377531/1‑76 (MQ=255) ttAATTGCCTGAAGGGTACCCAAGCGCCTTTTTTATTTCCCCACCTAAACCCCCCCCAGCCAGTAACTTCTCtttt < 2:1265851/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTCATCCTCCACCAGCCAGTAACTTCTCttttt > 1:1146206/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 1:1200019/1‑76 (MQ=255) tGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcg > 1:1425406/1‑76 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 1:1462407/1‑76 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 2:47460/1‑76 (MQ=255) cTGAAGAGTAACCCAACGGCCTTTTTATTTCCCCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTcgccg < 1:1482656/76‑1 (MQ=255) gAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCgccgcc > 1:513854/1‑76 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc > 2:1231680/1‑76 (MQ=255) gagTACCCCAACGGCCTTTTTATTCCACCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCta < 2:1702497/76‑2 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 2:330627/1‑76 (MQ=255) aaCCCAGCGGCGTTTTTTTTTCCCCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 1:357887/76‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 2:1792231/76‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:746463/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:161872/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 1:987989/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 1:987979/1‑76 (MQ=255) aCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCa > 1:1363359/1‑76 (MQ=255) ccAACGGCTTTTTTATTTCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc < 1:1427332/76‑1 (MQ=255) gCCTTTTTTTTTCCCCCCCAACTCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:976651/76‑1 (MQ=255) gCCTTTTTATTTCCCCACCTAACCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:942736/76‑1 (MQ=255) gCCTTTTTATTTCACCCCCTAATCCTCCCCCAGCCGGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:1548391/76‑1 (MQ=255) gCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:877662/76‑1 (MQ=255) gCCTTTTTACTTCACCTCCTAACCCGCCCCCAGCCAGTATCTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:1640653/76‑1 (MQ=255) ctttttttttTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:973251/74‑1 (MQ=255) ccTTTTTATTTCCCCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:973237/76‑1 (MQ=255) ccTTTTTATTTCACCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:1711552/76‑1 (MQ=255) ccTTTTTATTCCCCCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:60556/76‑1 (MQ=255) ccTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:1198472/76‑1 (MQ=255) ccTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:17141/76‑1 (MQ=255) cTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg < 2:1725321/76‑1 (MQ=255) ttttATTCCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc < 1:951426/76‑1 (MQ=255) tATTTCCCCCCCTACTCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 2:1153942/76‑1 (MQ=255) aTTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg > 2:1680818/1‑76 (MQ=255) aTTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg > 1:150530/1‑76 (MQ=255) tttCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt < 1:1303010/76‑1 (MQ=255) tttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt > 1:289151/1‑76 (MQ=255) ttCACCACCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa < 1:1087356/76‑1 (MQ=255) ttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa > 1:1706471/1‑76 (MQ=255) ttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa > 1:1706475/1‑76 (MQ=255) tCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa > 1:775282/1‑76 (MQ=255) cccccaccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa < 1:1743242/73‑1 (MQ=255) caccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa < 1:373139/76‑1 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 1:257945/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 1:200825/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:838891/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:1398537/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:544602/1‑76 (MQ=255) ccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAg > 1:1431639/1‑76 (MQ=255) | ATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631698‑2631845 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |