| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,921,703 | 0 | T | G | 29.0% | 149.7 / 24.5 | 69 | intergenic (‑88/+111) | ppiC/ECB_03654 | peptidyl‑prolyl cis‑trans isomerase C (rotamase C)/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (11/38); new base G (20/0); total (31/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.32e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGG > REL606/3921653‑3921753 | ttAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGt < 1:5313103/51‑1 (MQ=255) aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTAc < 1:6307821/51‑1 (MQ=255) gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACa < 1:5562144/51‑1 (MQ=255) cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACat < 1:618317/51‑1 (MQ=255) cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACat < 1:169498/51‑1 (MQ=255) cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACat < 1:5549869/51‑1 (MQ=255) tGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACatat < 1:6505620/51‑1 (MQ=255) gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACatata > 1:5726258/1‑51 (MQ=255) tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATg < 1:5437480/51‑1 (MQ=255) tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATg < 1:3399297/51‑1 (MQ=255) tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACatat > 1:2496384/1‑48 (MQ=255) ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATgg < 1:6168469/51‑1 (MQ=255) ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATgg < 1:175772/51‑1 (MQ=255) ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACata > 1:4918141/1‑46 (MQ=255) cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATgg > 1:5792594/1‑50 (MQ=255) gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCa < 1:6376339/51‑1 (MQ=255) cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaa < 1:115357/51‑1 (MQ=255) cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaa < 1:5970304/51‑1 (MQ=255) aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaa < 1:2929511/51‑1 (MQ=255) aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaa > 1:5964740/1‑51 (MQ=255) aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaa < 1:2962954/51‑1 (MQ=255) gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaaa < 1:5406214/51‑1 (MQ=255) gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaaa < 1:1224858/51‑1 (MQ=255) aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAg < 1:4021063/51‑1 (MQ=255) aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAg < 1:191841/51‑1 (MQ=255) aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAg < 1:1733260/51‑1 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGc < 1:5091459/51‑1 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGc < 1:1919126/51‑1 (MQ=255) tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaaa > 1:4655770/1‑48 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGcc < 1:2099435/51‑1 (MQ=255) gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaaaa > 1:2477630/1‑48 (MQ=255) cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGccc < 1:7046320/51‑1 (MQ=255) cTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGCAAAAGccc > 1:2280193/1‑51 (MQ=255) cTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaaaa > 1:3909484/1‑47 (MQ=255) cTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAAccc > 1:3650224/1‑51 (MQ=255) tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCg > 1:2386602/1‑51 (MQ=255) aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCgg < 1:972019/51‑1 (MQ=255) aCGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGccc > 1:654295/1‑48 (MQ=255) cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGaa < 1:118658/51‑1 (MQ=255) gTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGaaa > 1:161101/1‑51 (MQ=255) gTAATTGTCACTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaa > 1:3931552/1‑51 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGaaaa < 1:3065121/51‑1 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:1213324/1‑51 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:6730878/1‑51 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:2194677/1‑51 (MQ=255) tAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:2660287/1‑51 (MQ=255) ttGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa > 1:1273399/1‑48 (MQ=255) gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg < 1:5544287/51‑1 (MQ=255) gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg < 1:6674725/51‑1 (MQ=255) tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCggg < 1:5820122/51‑1 (MQ=255) gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCtt < 1:6885845/51‑1 (MQ=255) gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCtt < 1:6923392/51‑1 (MQ=255) gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCtt < 1:1056676/51‑1 (MQ=255) gtcagGGAGTGGCTGGGACATATATGGCAAAAACCCGGaaaa > 1:1021114/1‑42 (MQ=255) tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCggg > 1:7120038/1‑47 (MQ=255) tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTc < 1:5728018/51‑1 (MQ=255) tcagGGAGTGGCTGGGACATATATGGCAAAAACCCGGaaaa > 1:551939/1‑41 (MQ=255) cagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCtt > 1:6747703/1‑49 (MQ=255) cagGGAGTGGCTGGTACATATATGGCAAAAACCCGGAAAATCCGGGCtt > 1:1229076/1‑49 (MQ=255) agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGc < 1:963112/51‑1 (MQ=255) gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa < 1:4967507/51‑1 (MQ=255) gggAGTGGCTGGTACATATATGGCAAAAGCCCGCAAAATCCGGGCTTCGCa < 1:5784771/51‑1 (MQ=255) ggAGTGGCTGGTACATATATGGCAAAAACCCGGAAAATCCGGGCtt > 1:4696967/1‑46 (MQ=255) gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg < 1:2790584/51‑1 (MQ=255) gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg > 1:6636839/1‑51 (MQ=255) gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg > 1:6171490/1‑51 (MQ=255) gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCt > 1:1618354/1‑51 (MQ=255) gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCt > 1:1758112/1‑51 (MQ=255) tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCa > 1:1592776/1‑51 (MQ=255) ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCaa > 1:5171160/1‑51 (MQ=255) tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg > 1:3768041/1‑51 (MQ=255) tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg < 1:5895208/51‑1 (MQ=255) | TTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGG > REL606/3921653‑3921753 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |