Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 34.3% | 63.0 / 14.6 | 35 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (3/20); new base G (12/0); total (15/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.45e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492009‑3492157 | gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggg > 2:837956/1‑76 (MQ=255) cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg > 2:804461/1‑75 (MQ=255) ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg < 1:771130/76‑1 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt < 1:495056/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 1:92315/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 2:405032/76‑1 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGtataa > 2:685491/1‑75 (MQ=255) gCGTCGCACAAAACAAGGGATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACg > 2:724195/1‑76 (MQ=255) gCGTCGCACAAAACAAGGGATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTAGCTGGGGGAGTATACg > 2:724174/1‑76 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag < 1:1222364/76‑1 (MQ=255) cacaAAACAAGGCATGTGGCTGAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag < 1:1055476/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 1:863845/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 1:954400/76‑1 (MQ=255) aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAg < 1:1118631/76‑1 (MQ=255) aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAg > 1:830092/1‑76 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt < 1:651107/76‑1 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTc > 2:44547/1‑76 (MQ=255) gcgtgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt < 1:733012/73‑1 (MQ=255) gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt < 1:1043986/76‑1 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTg > 1:528491/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTgg < 2:534383/76‑1 (MQ=255) tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGc < 2:468273/76‑1 (MQ=255) gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcg < 1:1040994/76‑1 (MQ=255) gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTAAAAGCAACAAAAACTGTCGCTGGcgtcg > 1:84260/1‑76 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgg > 2:974115/1‑75 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg > 2:1050588/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 2:1101682/76‑1 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGa > 1:259529/1‑76 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 1:438244/76‑1 (MQ=255) gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa < 1:342379/76‑1 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa < 1:569020/76‑1 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa < 1:1121247/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 1:53565/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 2:22209/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 1:1007030/76‑1 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTTGCTGGGGGCGTGGTTTTATGTAAg > 2:464811/1‑76 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGGGGAATAAACGCAACAAAAACTGTCGCTGGCGGCGTGAATTAAAGTAAg > 2:194662/1‑76 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg < 2:338805/76‑1 (MQ=255) tctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGc < 1:478785/76‑1 (MQ=255) ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTTATGTAAGgggggc > 2:704001/1‑74 (MQ=255) ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTTATGTAAGGGGGcc > 1:147931/1‑76 (MQ=255) gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTAATGTAAGGTGGCCAgg > 2:764240/1‑76 (MQ=255) gtgtGAGTCAGTGGCTGGGGGGGGATACGCCGCAGAAGCTGTCGCTGGCGTCGGGGTTTAATGTCAGGGGGCCAgg > 2:764249/1‑76 (MQ=255) tgtgAGTCAGTGGCTGGTGGAGTATACGAAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGt < 2:144379/76‑1 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt > 2:11590/1‑76 (MQ=255) tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc < 1:31404/76‑1 (MQ=255) | GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492009‑3492157 |
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: ‑ |