Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,495,872 | 0 | G | A | 30.4% | 81.6 / 20.3 | 46 | S74L (TCG→TTG) | ansP | L‑asparagine transporter |
Reads supporting (aligned to +/- strand): ref base G (23/9); new base A (0/14); total (23/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.81e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGAAGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGCCATT > REL606/1495822‑1495922 | tGGAAGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACg > 1:7097295/1‑51 (MQ=255) tGGAAGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACg > 1:2266780/1‑51 (MQ=255) ggAAGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGa > 1:4163942/1‑51 (MQ=255) aaGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGaaa < 1:390609/51‑1 (MQ=255) aaGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGaaa < 1:5702024/51‑1 (MQ=255) aGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGaaaa < 1:4319075/51‑1 (MQ=255) gggCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAAc < 1:2354603/51‑1 (MQ=255) gggCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAAc < 1:3069778/51‑1 (MQ=255) gggCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAAc < 1:4480975/51‑1 (MQ=255) gCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACaa < 1:1968215/51‑1 (MQ=255) gCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACaa < 1:1179828/51‑1 (MQ=255) ggTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGc < 1:5386025/51‑1 (MQ=255) ggTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGc > 1:6899710/1‑51 (MQ=255) gCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACaaa > 1:5977794/1‑51 (MQ=255) gCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACaaa > 1:3332681/1‑51 (MQ=255) aCCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAAtt > 1:2223557/1‑51 (MQ=255) cTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATtaaata > 1:2924/1‑51 (MQ=255) tCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATtaaataa > 1:1223046/1‑51 (MQ=255) ccAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAcca > 1:5144852/1‑51 (MQ=255) ccAATGCACGCAGAAAAAAAAACAAAAAAAAGCCACAAATTAAATAAAcca < 1:3595678/51‑1 (MQ=255) cAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccag > 1:2285895/1‑51 (MQ=255) cAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccag > 1:6078278/1‑51 (MQ=255) aaTGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccagt > 1:1696579/1‑51 (MQ=255) tGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccagtgc > 1:1093428/1‑51 (MQ=255) gCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccagtgcc > 1:7110540/1‑51 (MQ=255) gCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAAccagtgccag > 1:387727/1‑49 (MQ=255) cAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAACcagtgccagtgc > 1:3488583/1‑51 (MQ=255) cAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAACcagtgccagtgc > 1:6523747/1‑51 (MQ=255) aaaaTAAAAAACAAAAAAAAGCCACAAATTAAATAAAccagtgccagtg < 1:6752050/47‑1 (MQ=38) aaaaaaaaaCAAAAACAAGCCACAAATTAAATAAAcc < 1:1818277/34‑1 (MQ=25) aaaaaaaaaCAAAAACAAGCCACAAATTAAATAAAcc < 1:1671829/34‑1 (MQ=25) aaaaaaaaaCAAAAACAAGCCACAAATTAAATAAAc < 1:4134751/33‑1 (MQ=25) aaaaaaaaaCAAAAAAAAGCCACAAATTAAATAAAccagtgccagt < 1:636335/43‑1 (MQ=38) aaaaaaaaaCAAAAAAAAGCCACAAATTAAATAAAccagtgccag < 1:5860771/42‑1 (MQ=38) aaaaaaaaaCAAAAAAAAGCCACAAATTAAATAAAccagt < 1:1491008/37‑1 (MQ=37) aaaaaaaaaCAAAAAAAAGCCACAAATTAAATAAACcagtgccagtgc < 1:6652855/45‑1 (MQ=38) aaTAAAAAACGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTgg > 1:1455241/1‑51 (MQ=255) aaaaaaaaCAAAAAAAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGa < 1:735267/49‑1 (MQ=255) tAAAAAACGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGAc > 1:6392889/1‑51 (MQ=255) aaaaaaaCAAAAAAAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGAc < 1:3409409/50‑1 (MQ=255) aaaaaaCAAAAAAAAGCCACAAATTAAATAAAccagtgccagt < 1:1816292/43‑1 (MQ=38) aaaaaaCAAAAAAAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGa < 1:100474/49‑1 (MQ=39) aaaaaaCAAAAAAAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGa < 1:135030/49‑1 (MQ=39) aaaaacaaaaacAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGAccc < 1:5743900/51‑1 (MQ=255) aaaaCGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCg > 1:4300034/1‑51 (MQ=255) aaaaCGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCg > 1:404053/1‑51 (MQ=255) aaaCGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGc > 1:443150/1‑51 (MQ=255) aaaCGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGc > 1:3009501/1‑51 (MQ=255) aaaCGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGc > 1:3354373/1‑51 (MQ=255) gAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGCCAtt > 1:3649567/1‑51 (MQ=255) gAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGCCAtt > 1:5605621/1‑51 (MQ=255) gAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGCCAtt > 1:798730/1‑51 (MQ=255) | TGGAAGGGCGGTGTAGCACCAGCTCACCCAATGCACGCAGAATAAAAAACGAAAACAAGCCACAAATTAAATAAACCAGTGCCAGTGCTGGACCCGCCATT > REL606/1495822‑1495922 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |