| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 41.4% | 32.1 / 13.0 | 29 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (5/12); new base G (12/0); total (17/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.23e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA > REL606/1484031‑1484179 | ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAt < 2:1186160/76‑1 (MQ=255) tACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtt > 1:912701/1‑76 (MQ=255) ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACttt < 2:259243/76‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 1:1056506/76‑1 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTACGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTAAttttta > 1:1394781/1‑74 (MQ=255) tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGAttt > 2:754532/1‑76 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGAtttt > 2:1048530/1‑75 (MQ=255) tgtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACTTTTTTCGATTTCAc > 2:550001/1‑76 (MQ=255) tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa < 1:643901/76‑1 (MQ=255) tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGt < 1:379716/76‑1 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTc < 1:965567/76‑1 (MQ=255) cATTGTCTNTTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa < 1:1349465/76‑1 (MQ=255) tGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCGGGAACg > 1:735344/1‑76 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTTTCGGGAACGCCg > 2:281619/1‑76 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATATTTTTTACGATTTTACCATGTCGGGAAAGCCg > 1:1287699/1‑76 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGCTTTTTTTCTTTTTACGATTTTACCCTTTCGNGAAAGCCg > 1:1352401/1‑76 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGATTTTACCATGTCGGTAACGCCgg > 1:1212512/1‑76 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGATTTCACCATGTCGGGAAAGCCggca > 1:1489862/1‑74 (MQ=255) ttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATTTCGGGAAAGCCGGGAc > 1:307927/1‑76 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTGACGATTTTTCCACGTCGGGAAAGCCCGAACTg > 2:1068423/1‑76 (MQ=255) cTCGCTGCTGACACTGTTATGGCAGGGGGTTTAAACTTTTGACGATTTCACCCTGTCGGTAACGCCGGCACTTACg > 1:7797/1‑76 (MQ=255) cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg < 2:1249385/76‑1 (MQ=255) cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg < 2:671438/76‑1 (MQ=255) tCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGc > 2:1472789/1‑76 (MQ=255) tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc > 2:890061/1‑76 (MQ=255) gctgGCACTGTTATGGCAGGGGTTTTATATTTTTGACGATTTCACCATGTCGGTAACGCCGGCACTGACGCTGGCg > 2:1521192/1‑76 (MQ=255) gctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg > 1:1302351/1‑76 (MQ=255) ctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGa < 1:1139011/76‑1 (MQ=255) acTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc > 2:664070/1‑76 (MQ=255) ggCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCt > 1:1164401/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 2:270965/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 2:706309/1‑76 (MQ=255) ggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa < 1:943471/76‑1 (MQ=255) gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa < 1:1059884/76‑1 (MQ=255) | TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA > REL606/1484031‑1484179 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |