Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 48.9% | 49.4 / 24.1 | 47 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (3/21); new base G (23/0); total (26/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTG > REL606/3492011‑3492154 | ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg < 1:729689/76‑1 (MQ=255) ccTGTGACGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCATCAATCTCGTGTGAGTCAGTGGCTGGTGGAg < 1:1284526/76‑1 (MQ=255) tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGta < 2:453569/76‑1 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 2:895596/1‑76 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 1:736011/1‑76 (MQ=255) gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGNGATAAgca > 1:1232202/1‑76 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc < 1:1183466/76‑1 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc < 1:417878/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 1:190981/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 2:844268/76‑1 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATAAGCAGCAGaa > 2:10397/1‑76 (MQ=255) aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGc > 2:617434/1‑76 (MQ=255) aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGc > 2:810375/1‑76 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCt > 2:1082320/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTc > 2:151595/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGGAGCTGtg > 2:296173/1‑75 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTgg > 2:1368513/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTgg > 1:819507/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCCgg > 1:276499/1‑76 (MQ=255) gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcg < 1:1047788/76‑1 (MQ=255) gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCCgggg > 1:954358/1‑74 (MQ=255) gCTAAATTCGGATCAACAATGTNNTGTGAGTCAGNNTCTGGTGGTGTATAGGCAGCAGAAGCTGTCGCTGGcgtcg < 1:1262280/76‑1 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGgggcgt > 2:66946/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 2:747808/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg < 1:817215/76‑1 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAACTGTCGCTGGCGGCGTGa > 1:1047777/1‑76 (MQ=255) aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtt < 1:708189/76‑1 (MQ=255) tGCGGATCAACAATCTTGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 2:1249515/76‑1 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 1:68237/76‑1 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 2:1095624/76‑1 (MQ=255) gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa < 1:387583/76‑1 (MQ=255) gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa > 1:554236/1‑76 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGATTTTATGTaa > 1:561966/1‑76 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg < 1:1001150/76‑1 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAg > 2:1019439/1‑76 (MQ=255) caaTCTCTGGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGt < 2:554448/76‑1 (MQ=255) caaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAggg > 1:574132/1‑75 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgggg > 2:818017/1‑74 (MQ=255) aTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGgggg > 1:1167877/1‑76 (MQ=255) tctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGGTGGc > 1:481804/1‑76 (MQ=255) ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCCGAAGCTGTCGCTGGCGTCGGGATTTAATGTAAGGTGGcc > 2:161748/1‑76 (MQ=255) tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGt < 1:571179/76‑1 (MQ=255) gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 2:577466/76‑1 (MQ=255) aGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGGNTCNTGATTTAATGTAAGNTGGCCAGGTCAGa > 2:1264426/1‑76 (MQ=255) aGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGa < 1:1065643/76‑1 (MQ=255) aGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGGGGCCAGGGCAGa > 2:155919/1‑76 (MQ=255) cAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAttt < 2:1070985/76‑1 (MQ=255) cAGTGGCTGGGGGAGTACACGCCGCAGAAGCTGTCGCTGGCGTCGTGATTTTATGTAAGGGAGCCAGGGCAGAttt > 2:215343/1‑76 (MQ=255) ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGGTTTCCAt > 2:592911/1‑76 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg < 1:810375/76‑1 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg < 1:264350/76‑1 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg < 1:617434/76‑1 (MQ=255) | CCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTG > REL606/3492011‑3492154 |
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: ‑ |