Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 40.0% | 66.9 / 23.9 | 50 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (25/5); new base C (0/20); total (25/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-09 | |||||||||||
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:1076740/76‑1 (MQ=255) tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca > 1:556191/1‑76 (MQ=255) tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca < 1:1813861/76‑1 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac > 1:1360222/1‑76 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac < 2:229584/76‑1 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca > 1:626989/1‑76 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca < 2:1844581/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 2:1443262/76‑1 (MQ=255) gtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGACTTTTTATTTTACCTCCTAATCCTCCAccagcca < 2:1732882/76‑1 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag < 1:1995976/76‑1 (MQ=255) aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCCCCACCTAATCCCCCCCCAGCCAGTaa < 2:440123/76‑1 (MQ=255) tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAACCCCCCACCAGCCAGTAACtt < 2:1733841/76‑1 (MQ=255) cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct > 1:725384/1‑76 (MQ=255) cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct > 1:2707/1‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 1:526896/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 1:2045005/1‑76 (MQ=255) cccTGAAGAGTAACCCACCGCCCTTTTTATTTCCCCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTcgc < 1:1767967/75‑1 (MQ=255) gAAGAGTACCCCACCGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc < 2:626340/76‑1 (MQ=255) gAAGAGTACCCCACCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc < 2:552431/76‑1 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc > 2:1345153/1‑76 (MQ=255) agagTACCCCAACGGCCTTTTTATTTCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:13954/76‑1 (MQ=255) agagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt > 1:1549517/1‑76 (MQ=255) agagTAACCCAACGGCCTTTTTATTTCA‑CACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCtc > 1:1256824/1‑75 (MQ=255) gagTACCCCAACGGCCTTTTTATTTCCCCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 2:720195/76‑1 (MQ=255) gTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg > 2:1879259/1‑76 (MQ=255) taccccAACGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt < 2:1005240/73‑1 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 2:1791192/1‑76 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 2:1596604/1‑76 (MQ=255) cccAACGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:890256/76‑1 (MQ=255) cccAACGGCCTTTTTATTTCCCCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:147950/76‑1 (MQ=255) ccAACGGCCTTTTTATTTCCCCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc < 1:254424/76‑1 (MQ=255) ccAACGGCCTTTTTATTTCCCCACCTAACGCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc < 2:1991803/76‑1 (MQ=255) aCGGCCTTTTTATTCCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtg < 1:1590522/76‑1 (MQ=255) gcttttttATTTCACCTCCAATCCCTCCGCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:1913635/73‑1 (MQ=255) gcttttttATTCCCCCCCCAACCCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:326149/73‑1 (MQ=255) gcttttttATTCCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:1374910/73‑1 (MQ=255) gcttttttATTCCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:616813/73‑1 (MQ=255) gCCTTTTTATTCCCCCCCCTAGTCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 2:1233182/76‑1 (MQ=255) cttttttATTCCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:1831988/74‑1 (MQ=255) cttttttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 1:695037/74‑1 (MQ=255) cTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg > 1:905985/1‑76 (MQ=255) cTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg > 1:905989/1‑76 (MQ=255) ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc > 2:1784309/1‑76 (MQ=255) ttttATTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc < 1:1026159/76‑1 (MQ=255) ttttATTCCCCCACCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc < 2:797178/76‑1 (MQ=255) tATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCCCt > 1:723226/1‑76 (MQ=255) tCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa > 1:1200186/1‑76 (MQ=255) accccaccaAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa < 1:357220/73‑1 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 2:29371/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 1:881467/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1513780/1‑76 (MQ=255) aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa > 2:1045504/1‑76 (MQ=255) aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa > 2:1991128/1‑76 (MQ=255) aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa > 2:1279331/1‑76 (MQ=255) tcctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaa > 2:1222454/1‑76 (MQ=255) cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa > 1:161429/1‑76 (MQ=255) ctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAAt > 1:1443711/1‑76 (MQ=255) tccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCGCTGTGACCAATTTTAAAATa > 1:1239958/1‑76 (MQ=255) caccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGc > 2:1568387/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGC > REL606/2631697‑2631846 |
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: ‑ |