| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,631,772 | 0 | A | C | 27.1% | 102.4 / 12.2 | 48 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
| Reads supporting (aligned to +/- strand): ref base A (28/7); new base C (0/13); total (28/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.02e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631697‑2631845 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca < 2:972722/76‑1 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac > 2:788828/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca > 1:122718/1‑76 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 1:1209804/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 1:35094/76‑1 (MQ=255) gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTa > 2:93834/1‑76 (MQ=255) gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTa > 2:1334426/1‑76 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt > 1:1680377/1‑76 (MQ=255) gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtc > 1:666254/1‑76 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct < 2:748637/76‑1 (MQ=255) cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct > 1:1606389/1‑76 (MQ=255) aTTTAACTGCCTGAAGAATACCCCAACGGCCTTTTTATTTTACCACCTAATCCTCCCCCAGCCAGTAACTTCTCtt < 2:1231312/76‑1 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 1:135720/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCCCCCCCAGCCAGTAACTTCTCttttt < 1:735567/76‑1 (MQ=255) aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc > 2:204516/1‑76 (MQ=255) aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCCCCACCAGCCAGTAACTTCTCTTTTtct < 2:122718/76‑1 (MQ=255) cTGCCTGAAGAGTACCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtctc > 2:862524/1‑76 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 2:687855/1‑76 (MQ=255) tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc > 1:1337552/1‑76 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCCCCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc < 1:1685344/76‑1 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc > 1:1583940/1‑76 (MQ=255) agagTACCCCACCGGCCTTTTTATTTCCCCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 2:284247/76‑1 (MQ=255) gagTACCCCAACGGCTTTTTTATTCCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 2:345115/76‑1 (MQ=255) gagTACCCCAACGGCCTTTTTATTTCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 1:403325/76‑1 (MQ=255) gTACCCCAACGGCCTTTTTATTTCACCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg < 2:1550169/76‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCCCCCCCTAGCCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 2:1583940/76‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:374162/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:587368/1‑76 (MQ=255) cccACCGCCGTTTTTATTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 1:819206/76‑1 (MQ=255) ccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc > 2:872675/1‑76 (MQ=255) aaCGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt > 2:1304347/1‑76 (MQ=255) aaCGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt > 2:1490421/1‑76 (MQ=255) cGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCGGCgtgt > 2:448091/1‑76 (MQ=255) ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt > 1:593230/1‑76 (MQ=255) gcttttttATTCCCCCTCCTATCCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:676738/73‑1 (MQ=255) gCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:892391/76‑1 (MQ=255) gCCTTTTTATTCCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:633552/76‑1 (MQ=255) ccTTTTTATTTCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:94483/76‑1 (MQ=255) ccTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa > 2:1731274/1‑76 (MQ=255) ccTGTTTATTTCCCCACCAAATCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:94467/76‑1 (MQ=255) tttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGc > 1:1663124/1‑76 (MQ=255) tttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa > 2:1758490/1‑76 (MQ=255) tATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt > 2:590838/1‑76 (MQ=255) tATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 2:1366712/76‑1 (MQ=255) tttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt > 2:255918/1‑76 (MQ=255) ttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa > 1:1841675/1‑76 (MQ=255) caccacCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAAc > 1:752420/1‑76 (MQ=255) ccaccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaa < 2:1514/76‑1 (MQ=255) cccTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAat < 1:972530/75‑1 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 2:897661/1‑76 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 2:1111000/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:754434/1‑76 (MQ=255) ccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAg > 2:1562348/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631697‑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: ‑ |