| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 157,994 | 0 | G | T | 30.6% | 174.9 / 34.8 | 72 | Q93K (CAA→AAA) | htrE | predicted outer membrane usher protein |
| Reads supporting (aligned to +/- strand): ref base G (8/42); new base T (22/0); total (30/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.27e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGAGAATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTTAATAA > REL606/157945‑158044 | ggAGAATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTtgt < 1:1818824/51‑1 (MQ=255) agaATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTtgt > 1:5158698/1‑49 (MQ=255) agaATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGc < 1:4806023/51‑1 (MQ=255) agaATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGc < 1:5221719/51‑1 (MQ=255) gaATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCg < 1:183206/51‑1 (MQ=255) aTTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGTGtt < 1:3896631/51‑1 (MQ=255) tttATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGttt > 1:2821011/1‑51 (MQ=255) tttATATGAAACTACAATAAATTCTTAAGTGTGATACAGGCTTGTGCGttt < 1:6421211/51‑1 (MQ=255) ttATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGtttt < 1:5052952/51‑1 (MQ=255) ttATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGtttt < 1:2571205/51‑1 (MQ=255) tataTGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGttttt > 1:5502066/1‑51 (MQ=255) atatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGtttttt > 1:5529813/1‑51 (MQ=255) atatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGgttttt > 1:6274705/1‑51 (MQ=255) atatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGtttttt < 1:4983219/51‑1 (MQ=255) atatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGtttttt < 1:6851475/51‑1 (MQ=255) tatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGttttttt < 1:1044674/51‑1 (MQ=255) tatGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGttttttt < 1:3663070/51‑1 (MQ=255) atGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGgttttttt > 1:2559598/1‑51 (MQ=255) atGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGtttttttt < 1:6940882/51‑1 (MQ=255) tGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGgtttttttt > 1:1409861/1‑51 (MQ=255) aaaCTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGttttttttt > 1:5569359/1‑49 (MQ=255) aaaCTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGttttttttt > 1:3960478/1‑49 (MQ=255) aaaCTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTcc < 1:2446052/51‑1 (MQ=255) aaCTGCAATAAATTCTTAAGTGTGATGCAGGCTTGTGCGTTTTTTTTTCCt < 1:6957001/51‑1 (MQ=255) aCTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGgttttttttt > 1:108136/1‑47 (MQ=38) aCTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCtt < 1:95690/51‑1 (MQ=255) aCTGCAATAAATTCTTAAGTGTGATACAGGCTCGTGCGTTTTTTTTTCCtt < 1:1289689/51‑1 (MQ=255) cTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGttttttttt > 1:74990/1‑46 (MQ=255) cTGCAATAAATTCTTAAGTGTGATACAGGCTTTTGCGgttttttttt > 1:6002044/1‑46 (MQ=38) cTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGttttttttttcctt < 1:1215461/51‑6 (MQ=255) cTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttc < 1:2582482/51‑1 (MQ=255) tGCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGtttttttttt > 1:3000290/1‑45 (MQ=255) gCAATAAATTCTTAAGTGTGATACAGGCTTTTGGGgtttttttt > 1:4757569/1‑44 (MQ=38) gCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaa < 1:3679548/51‑1 (MQ=255) gCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaa < 1:4323012/51‑1 (MQ=255) cAATAAATTCTTAAGTGTGATACAGGCTTGTTCGgtttttttt > 1:2042113/1‑43 (MQ=38) cAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaat < 1:5681291/51‑1 (MQ=255) aaTAAATTCTTAAGTGTGATACAGGCTTTTGGGtttttttttt > 1:1542322/1‑42 (MQ=255) aaTAAATTCTTAAGTGTGATACAGGCTTTTGGGttttttttt > 1:6753210/1‑42 (MQ=255) aaTAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaatt < 1:6883973/51‑1 (MQ=255) aTAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaattt < 1:1715139/51‑1 (MQ=255) aTAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaattt < 1:3104994/51‑1 (MQ=255) aTAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaattt < 1:4442821/51‑1 (MQ=255) tAAATTCTTAAGTGTGATACAGGCTTTTGGGgttttttttt > 1:704304/1‑40 (MQ=37) aaaTTCTTAAGTGTGATACAGGCTTTTGCGTTTTTTTTTTCtt > 1:6767495/1‑43 (MQ=38) aaaTTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaatttca < 1:2308643/51‑1 (MQ=255) aaaTTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaatttca < 1:7249853/51‑1 (MQ=255) aaTTCTTAAGTGTGATACAGGCTTTTGTGttttttttt > 1:5622674/1‑38 (MQ=38) aaTTCTTAAGTGTGATACAGGCTTTTGCGgttttttttt > 1:3722704/1‑38 (MQ=37) aTTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttcaatttcaat < 1:6068853/51‑1 (MQ=255) ttCTTAAGTGTGATACAGGCTTTTGGGttttttttt > 1:6316412/1‑36 (MQ=38) ttCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATa < 1:7162936/51‑1 (MQ=255) ttCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATa < 1:1548079/51‑1 (MQ=255) ttCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATa < 1:5721608/51‑1 (MQ=255) tCTTAAGTGTGATACAGGCTTTTGCGgtttttttt > 1:2360372/1‑35 (MQ=37) tCTTAAGTGCGATACAGGCTTTTGGGgttttttttt > 1:1069440/1‑35 (MQ=16) cTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATaaa < 1:5470811/51‑1 (MQ=255) cTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATaaa < 1:5196253/51‑1 (MQ=255) tAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATg < 1:5311496/51‑1 (MQ=255) tAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATg < 1:220238/51‑1 (MQ=255) aGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTa < 1:3489060/51‑1 (MQ=255) tgtgATACGGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAAt < 1:6865528/51‑1 (MQ=255) tgtgATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAAt < 1:6410509/51‑1 (MQ=255) gtgATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATg < 1:5340834/51‑1 (MQ=255) gtgATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATg < 1:6843766/51‑1 (MQ=255) tgATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGc < 1:192013/51‑1 (MQ=255) cAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTgg > 1:2669500/1‑51 (MQ=255) cAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTgg < 1:5241810/51‑1 (MQ=255) aGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGt > 1:6607077/1‑51 (MQ=255) ggCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGtt > 1:7197014/1‑51 (MQ=255) cTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaa > 1:687149/1‑51 (MQ=255) ttGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaat > 1:6741449/1‑51 (MQ=255) ttGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaat > 1:3892502/1‑51 (MQ=255) ttGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaat > 1:3990601/1‑51 (MQ=255) tgtgCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaata > 1:358515/1‑51 (MQ=255) tgtgCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaata > 1:3226484/1‑51 (MQ=255) gtgCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTtaataa > 1:4051156/1‑51 (MQ=255) | GGAGAATTTATATGAAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTAATGCTTTGGTTAATAA > REL606/157945‑158044 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |