Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 60,574 | 0 | T | G | 28.8% | 127.4 / 10.0 | 52 | Y244S (TAC→TCC) | imp | organic solvent tolerance protein precursor |
Reads supporting (aligned to +/- strand): ref base T (2/35); new base G (15/0); total (17/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.03e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAA > REL606/60525‑60624 | cgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTa < 1:2394086/51‑1 (MQ=255) gTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAg < 1:2043964/51‑1 (MQ=255) tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGa < 1:3372717/51‑1 (MQ=255) gATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGaa > 1:4115778/1‑51 (MQ=255) tggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACt < 1:2862770/51‑1 (MQ=255) tggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACt > 1:5262151/1‑51 (MQ=255) ggtggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTc < 1:3705265/51‑1 (MQ=255) tggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaa < 1:4496855/51‑1 (MQ=255) tggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaa > 1:4619550/1‑51 (MQ=255) ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa < 1:3603023/51‑1 (MQ=255) ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa < 1:5417895/51‑1 (MQ=255) ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaaa < 1:2932757/51‑1 (MQ=255) gCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCaaa > 1:956808/1‑50 (MQ=255) cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt < 1:6332557/51‑1 (MQ=255) cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt < 1:5030349/51‑1 (MQ=255) cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt < 1:4923006/51‑1 (MQ=255) aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagta < 1:1509713/51‑1 (MQ=255) aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagta < 1:6776533/51‑1 (MQ=255) tCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaaa > 1:3396212/1‑47 (MQ=255) ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagtagt < 1:6505784/51‑1 (MQ=255) ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCaaa > 1:5661693/1‑46 (MQ=255) cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt < 1:3635103/51‑1 (MQ=255) cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt < 1:2872718/51‑1 (MQ=255) atTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggt < 1:5562149/51‑1 (MQ=255) tttGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCaaa > 1:5585170/1‑41 (MQ=38) gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggtggtgg < 1:2022789/51‑1 (MQ=255) cgcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTggtggtggt < 1:6357801/51‑1 (MQ=255) cgcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTggtggtggt < 1:3492088/51‑1 (MQ=255) cgcgATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggggg > 1:2041163/1‑47 (MQ=255) gcgATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtggggg > 1:5555698/1‑49 (MQ=39) gcgATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTAGTtggtggggg > 1:6029401/1‑49 (MQ=38) gcgATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTtggtggggg > 1:5610398/1‑49 (MQ=255) cgATGTTCCAGTAATATGGCAGGGAGAACTCaaa > 1:4491505/1‑34 (MQ=255) cgATGTTCCAGTAATATGGCAGGGAGAAATCaaa > 1:965425/1‑34 (MQ=37) gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa < 1:1378084/51‑1 (MQ=255) gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa < 1:1391493/51‑1 (MQ=255) gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTa < 1:837210/51‑1 (MQ=255) aTGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTAc < 1:6523665/51‑1 (MQ=255) tGTTCCAGTAATATGGCAGGGAGAACTCaaa > 1:6967430/1‑31 (MQ=255) gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACtt < 1:2219346/51‑1 (MQ=255) ttCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtggggg > 1:3841973/1‑43 (MQ=38) ccAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc < 1:4120052/51‑1 (MQ=255) ccAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc < 1:5679701/51‑1 (MQ=255) atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCgg < 1:308700/51‑1 (MQ=255) atGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa < 1:2599683/51‑1 (MQ=255) atGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa < 1:7282140/51‑1 (MQ=255) tGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGAt < 1:72198/51‑1 (MQ=255) ggCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATc < 1:2191725/51‑1 (MQ=255) cAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCaa < 1:34187/51‑1 (MQ=255) aGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAg < 1:4133036/51‑1 (MQ=255) gTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa > 1:3891494/1‑51 (MQ=255) gTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa > 1:6280803/1‑51 (MQ=255) tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaaa < 1:829727/51‑1 (MQ=255) tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaaa > 1:3197095/1‑51 (MQ=255) | CGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAA > REL606/60525‑60624 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |