| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 37.0% | 32.7 / 15.0 | 27 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (2/15); new base G (10/0); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.82e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTA > REL606/1484030‑1484178 | tttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGa < 2:1310243/76‑1 (MQ=255) ctgctgCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAtttt < 2:26644/76‑1 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATtttt > 1:657892/1‑76 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtttt < 2:1544063/76‑1 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTACGGCAGGGGTTTTTTTCtttt > 2:886453/1‑76 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCtttttt > 1:1371342/1‑74 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACtttttt > 1:141005/1‑74 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGa < 1:530426/76‑1 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTTCGATTTc > 1:1326935/1‑76 (MQ=255) gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCAcc < 1:787787/76‑1 (MQ=255) gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAATTTTTACGATTTTACCATg > 1:1323790/1‑76 (MQ=255) gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACTTTTTACGATTTCACCATGTCGGt > 1:173390/1‑76 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa < 1:1396176/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa < 2:920009/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa < 2:1197755/76‑1 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCGGGAACGcc > 1:134760/1‑76 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGGAACGcccg > 2:1455688/1‑74 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTATACTTTTTACGATTTCACCATGTCGGGAAAGCCgg > 1:275774/1‑76 (MQ=255) ttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGTAAAGCCGGGAc > 2:369162/1‑76 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTTTTCTTTTTACGGTTTTACCATGTGGGGAAGGCCGGGAcgg > 2:1535552/1‑74 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTg < 2:800922/76‑1 (MQ=255) ggCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTTACCATGTCGGTAACGCCGGCACTGa > 1:1306168/1‑76 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGGAACGCCGGAACTGACGCt > 1:614852/1‑76 (MQ=255) ctgctgACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGGACTGAAGCTgg > 2:435417/1‑76 (MQ=255) acacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAata > 2:1241918/1‑76 (MQ=255) cacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAatat < 2:815382/76‑1 (MQ=255) cTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATcc < 2:779005/76‑1 (MQ=255) ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTg < 2:665994/76‑1 (MQ=255) tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc < 1:84035/76‑1 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 1:313625/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt < 2:1090405/76‑1 (MQ=255) gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt < 2:1268508/76‑1 (MQ=255) ggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa < 2:346020/76‑1 (MQ=255) | TTTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTA > REL606/1484030‑1484178 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |