Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 46.3% | 56.6 / 11.9 | 41 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (2/20); new base G (19/0); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.40e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 | acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt < 1:250368/76‑1 (MQ=255) acgCCCTGTGGCGNCGCACAAAACAAGGNATNTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt < 2:1165950/76‑1 (MQ=255) cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtg < 1:799747/76‑1 (MQ=255) cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGa < 2:524793/76‑1 (MQ=255) ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg > 2:1292934/1‑74 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGgata > 2:218220/1‑76 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATAc > 1:877074/1‑76 (MQ=255) cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag < 2:111971/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 1:1013598/76‑1 (MQ=255) aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAg < 1:873380/76‑1 (MQ=255) aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGc < 1:921262/76‑1 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt < 1:947504/76‑1 (MQ=255) aCAAGGTATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTg < 1:1190119/76‑1 (MQ=255) aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTg > 2:1059266/1‑76 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGt < 1:812562/76‑1 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCCGAAGCGGt > 1:50267/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCTACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTc > 1:796714/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTc < 1:406701/76‑1 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTAGccg > 1:497217/1‑74 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTg > 1:497224/1‑76 (MQ=255) tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgg > 2:923128/1‑75 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGcggc > 1:263335/1‑74 (MQ=255) cTAAATGCGGATGAACAATCTCGTGTGAGTCANTGCCTGGGNGAGTATACGCGGCAGAAGCNGGCGCTGGcgtggg > 2:1080873/1‑73 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGgggcgt > 1:969100/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 2:153973/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 1:795812/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 1:234810/76‑1 (MQ=255) gaaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 2:1127901/75‑1 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 2:713922/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 1:178998/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGGGATTTAATGt > 2:34634/1‑76 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGt > 2:703842/1‑76 (MQ=255) tCAACAATCTCGTGTGACTCAGTGGCTGGGGGAGTATACGCCGCAGAAGCGGTCGCCGGCGGCGTGATTAAAtnga > 2:1223034/1‑73 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGTTTTATGTaa > 2:1161052/1‑76 (MQ=255) acaaTCTCGTGTGAGTCAGTGGCTGGGCGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAAGGCAAgg > 1:1276163/1‑76 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGGTTTAATTTTAgggg > 2:248789/1‑74 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgggg > 1:641215/1‑74 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGGGGTTTAATGTAAgcgg > 1:267382/1‑73 (MQ=255) tctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTAATGTAAGGGGGc > 1:609241/1‑76 (MQ=255) ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGcc < 1:555391/76‑1 (MQ=255) ttgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 2:1065943/75‑1 (MQ=255) gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 1:108199/76‑1 (MQ=255) gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAg < 1:294068/76‑1 (MQ=255) aGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGTTTTATTTAAGGGGGGCCGGGCAGAgttt > 2:917912/1‑75 (MQ=255) gTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCa < 1:840757/76‑1 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa > 2:105676/1‑76 (MQ=255) ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCNGGCGTCGTGATTTAATGTAAGGNGGCCAGGTGANATTTCAAt > 2:1086846/1‑76 (MQ=255) cgggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt < 1:804403/74‑1 (MQ=255) ggtggAGTATACGCAGCAGAAGCTGTCGCTGGCNTCGTGATTTAATGTAAGGTGGCNAGGTCAGACTNCAATCTGa > 2:1095111/1‑76 (MQ=255) gtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAc < 1:978813/76‑1 (MQ=255) tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc < 2:635228/76‑1 (MQ=255) | ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |