| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 22.0% | 116.3 / 11.4 | 50 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (6/33); new base G (11/0); total (17/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.31e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTT > REL606/1484037‑1484176 | ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtata < 1:1326573/76‑1 (MQ=255) ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGCCACTGTTTTGGTTGGGTTTTTTTAttt > 2:920666/1‑76 (MQ=255) ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACtt > 2:1819677/1‑76 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 1:894811/76‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 2:670818/76‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCGATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 1:208928/76‑1 (MQ=255) ggCCCACTAATGTGCTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCCGGGATTGTATACTTTTg < 1:167217/76‑1 (MQ=255) cccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTtttttactttttac > 2:389982/1‑76 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGAt > 1:1036021/1‑76 (MQ=255) tANTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAntt < 2:1672307/76‑1 (MQ=255) tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAttt < 2:1395662/76‑1 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGTCGATTTCa > 2:813917/1‑76 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCa < 2:458810/76‑1 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGAAGGGGTTTTTTACTTTTGACGCtttta > 2:1250664/1‑74 (MQ=255) ggtgGTTTGGCATTGTCTATTTCGGCTCTCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCAc < 1:327753/75‑1 (MQ=255) ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTNTTATGGCAGGGGTTTTATACTTTTTACGCTTTCACCAt > 1:1994716/1‑76 (MQ=255) ggTTTGCCCTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTTCGATTTTACCAt > 1:1893833/1‑76 (MQ=255) gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCAtt > 1:220567/1‑75 (MQ=255) gTTTGGCATTGTCTAGTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATg < 1:1693125/76‑1 (MQ=255) ttGGCATTGTCTATTTCGTCTCGCTGCNGANAGNGTGATGGNAGAGATTTTATACTTTTGACGATTTCACCATGTc < 1:1772531/76‑1 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTc < 2:609738/76‑1 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTc < 1:1848605/76‑1 (MQ=255) tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCg < 2:907035/76‑1 (MQ=255) ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCgg < 1:662754/76‑1 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa < 2:1360570/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa < 1:471708/76‑1 (MQ=255) tGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTAAGATTTTTCCATGTCGGGaaag > 2:1488868/1‑74 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGc < 1:644821/76‑1 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGc < 1:1296580/76‑1 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTCACCATGTCGGGAAAGcc > 1:1813228/1‑76 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATTTCGGGAACGcc > 1:1229140/1‑76 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGGAACGcc > 1:1322426/1‑76 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCgg < 1:1549699/76‑1 (MQ=255) aTTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTTCGATTTCACCATGTCGGGAACGcccgg > 1:596592/1‑73 (MQ=255) tCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt < 1:576162/76‑1 (MQ=255) ggCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa < 2:593556/76‑1 (MQ=255) gCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTAGGTAACGCCGGAACTGAc < 1:1368321/76‑1 (MQ=255) cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg > 2:761540/1‑76 (MQ=255) gctgctGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg < 1:888460/76‑1 (MQ=255) gctgctGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg < 1:192115/76‑1 (MQ=255) tgctgACACTGTTGTGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc < 1:1059797/76‑1 (MQ=255) gctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg > 2:110748/1‑76 (MQ=255) acacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAata < 1:877459/76‑1 (MQ=255) acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc < 1:1171587/76‑1 (MQ=255) acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc < 1:1171572/76‑1 (MQ=255) ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTg < 2:835563/76‑1 (MQ=255) ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCCGTAACGCCGGAANTGACGCTGGCNAATATCCGTg < 2:2002813/76‑1 (MQ=255) aTGGCAGGGATTTTATACCTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcg > 1:1925856/1‑76 (MQ=255) ggCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCt < 2:1351216/76‑1 (MQ=255) ggCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCt < 1:410341/76‑1 (MQ=255) gCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGGGCTg > 1:1316721/1‑76 (MQ=255) cAGGGATTTTATACTTTTGACGCTTTCACCATGCCGGTAACGCCGTAACTGACGCTGGCGAATAGCGGTGCGCTGt > 2:130081/1‑76 (MQ=255) cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt > 2:353145/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 2:1054560/1‑76 (MQ=255) | CTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTT > REL606/1484037‑1484176 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |