Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 29.1% | 117.6 / 14.5 | 55 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (35/4); new base C (0/16); total (35/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.63e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCA > REL606/2631697‑2631847 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca < 2:15721/76‑1 (MQ=255) tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca > 2:1121656/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca < 1:1301815/76‑1 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 2:1117755/76‑1 (MQ=255) gtgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagc > 1:360118/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccggccag > 2:1484696/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag > 1:769004/1‑76 (MQ=255) aGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCTTTTTTATTTCACCCCCTAACCCCCCCCCAGCCAGt < 2:1257915/76‑1 (MQ=255) aGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGt > 1:1399671/1‑76 (MQ=255) gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTa > 1:1297601/1‑76 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa > 1:750105/1‑76 (MQ=255) aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAAc > 1:1314944/1‑76 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCCCCACCTAATCCTCCCCCAGCCAGTAACt < 2:281024/76‑1 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt > 1:1229034/1‑76 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt > 1:157254/1‑76 (MQ=255) gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCTTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTtc < 2:1224164/76‑1 (MQ=255) cccATTTAACTGCCTAAAGAGTAACCCAACGGCCTTTTTATTTCCCCCCCTAATCCCCCCCCAGCCAGTAACTtct < 2:452606/76‑1 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCtt < 2:202515/76‑1 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt > 1:326871/1‑76 (MQ=255) aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt > 1:353456/1‑76 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 1:1280013/1‑76 (MQ=255) tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 2:789645/1‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 2:1452330/1‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 2:718519/1‑76 (MQ=255) aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc > 1:721184/1‑76 (MQ=255) cTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTtctc < 2:3992/76‑1 (MQ=255) tGCCTGAAGAGTACTCCAACGGCCTTTTTATTTCCCCCCCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTcg < 2:1410748/76‑1 (MQ=255) tGCCTGAAGAGTACCCCAACGGCCTTTTTATTTCACCACCTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTcg < 2:158004/76‑1 (MQ=255) agagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt > 1:936865/1‑76 (MQ=255) gagTAACCCACCGGCCTTTTTATTCCACCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 1:324979/76‑1 (MQ=255) agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGc > 1:98985/1‑76 (MQ=255) gTACCCCACCGCCCTTTTTATTCCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg < 1:1121656/76‑1 (MQ=255) gTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg > 1:670507/1‑76 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 2:797685/1‑76 (MQ=255) cccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg > 2:1427688/1‑76 (MQ=255) caaaCGGCCTTTTTATTGCACCACCTACCCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAgt < 1:779053/74‑2 (MQ=255) ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt > 2:1112692/1‑76 (MQ=255) cttttttATTCCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:189124/74‑1 (MQ=255) ccTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa > 1:118333/1‑76 (MQ=255) cTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg < 1:1409188/76‑1 (MQ=255) tATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 1:106863/76‑1 (MQ=255) tATTTCACCACCTAATCCTCCACCAGTCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 2:1135096/76‑1 (MQ=255) tATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt > 1:854186/1‑76 (MQ=255) ttttCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg < 1:1429023/75‑1 (MQ=255) ttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa > 2:377115/1‑76 (MQ=255) tccaccccCTACTCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa < 1:395699/74‑1 (MQ=255) tccccaccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa < 1:223265/73‑1 (MQ=255) accaccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACa > 1:1166943/1‑76 (MQ=255) caccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa > 1:695542/1‑76 (MQ=255) ccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata > 2:178021/1‑76 (MQ=255) ccTAATCCTCCACCAGCCAGAAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata > 2:694728/1‑76 (MQ=255) ccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata < 1:1136757/76‑1 (MQ=255) ccTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata < 1:1059347/76‑1 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 1:1067569/1‑76 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAatat > 1:755741/1‑76 (MQ=255) cTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat < 2:4677/76‑1 (MQ=255) tAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt < 2:185198/76‑1 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1147085/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:728784/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGCAACAAATATTa > 1:814198/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 2:413391/1‑76 (MQ=255) ctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGCCAGCGTGTTTAGCCACTGTAACAAATATTAAAAt > 2:973147/1‑76 (MQ=255) accaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCa > 1:1352340/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCA > REL606/2631697‑2631847 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |