Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,568,967 | 0 | A | C | 47.2% | 28.8 / 20.6 | 36 | G1380G (GGT→GGG) | yfhM | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (14/5); new base C (0/17); total (14/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.24e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGT > REL606/2568892‑2569038 | cTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCa > 1:332537/1‑76 (MQ=255) cATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 2:437316/1‑76 (MQ=255) tCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCa < 2:624532/76‑1 (MQ=255) cGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAg < 1:99879/76‑1 (MQ=255) ggTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGcc < 1:84491/76‑1 (MQ=255) gTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCg < 2:576390/76‑1 (MQ=255) tttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGc > 2:26886/1‑76 (MQ=255) tttATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGc > 2:101634/1‑76 (MQ=255) aTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCtca < 2:565733/76‑1 (MQ=255) aTTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtca > 2:65841/1‑76 (MQ=255) ttttCTTCCAGCAAGGAGAGCATTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCtcat < 2:95794/76‑1 (MQ=255) ttttCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcat < 2:622390/76‑1 (MQ=255) ttcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCtcatca > 1:911344/1‑76 (MQ=255) tcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCAc > 2:941678/1‑76 (MQ=255) cttcCAGCAGGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:1062663/76‑1 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACt < 2:1090400/76‑1 (MQ=255) cttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:1054760/1‑76 (MQ=255) ccAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTAtt > 1:1039852/1‑76 (MQ=255) cAGCAAGGAGAGCATTACCGCTTTGCCGCGCAGTAAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTa < 2:1020174/76‑1 (MQ=255) cAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTa < 2:217773/76‑1 (MQ=255) cAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTa < 1:734678/76‑1 (MQ=255) aGCAAGGAGAGCATTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAc < 2:725648/76‑1 (MQ=255) aGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAc < 2:892204/76‑1 (MQ=255) aGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAc < 1:577925/76‑1 (MQ=255) gCAAGGAGAGCATTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAcg < 1:24637/76‑1 (MQ=255) cAAGGAGAGCATTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTAcgc < 2:783845/76‑1 (MQ=255) cAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcgc > 2:589819/1‑76 (MQ=255) aaGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcgcg > 2:408542/1‑76 (MQ=255) aGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACgcgg > 2:535118/1‑76 (MQ=255) agCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGttt < 1:693478/76‑1 (MQ=255) agCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGttt > 1:300519/1‑76 (MQ=255) cATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTc < 1:417147/76‑1 (MQ=255) aTTACCGCGTTGTCGCGCAGTGAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTca < 2:38707/76‑1 (MQ=255) tAACGCTTTGGCGCGCAGTAAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagc < 2:933313/76‑1 (MQ=255) aaCGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcg > 1:459605/1‑76 (MQ=255) gcgTTGCCGCGCAGTAAACTACCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcca < 2:808810/76‑1 (MQ=255) ccgTTGCCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcca < 2:949337/75‑1 (MQ=255) tgccgcgCAGTGACCTCCCTTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagcgc < 2:634502/73‑1 (MQ=255) tgccgccCAGTAACCTCCCGAAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccagcgc < 2:996519/73‑1 (MQ=255) gcAGTAAACTACCGTAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGc < 2:678501/76‑1 (MQ=255) tGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGcttct > 1:747715/1‑76 (MQ=255) aaTCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGt < 2:341437/76‑1 (MQ=255) | CTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGT > REL606/2568892‑2569038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |