| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 30.6% | 59.9 / 16.4 | 36 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (5/20); new base G (11/0); total (16/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.27e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTT > REL606/1484031‑1484177 | ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGa < 1:1431249/75‑1 (MQ=255) aCTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAtt < 1:990488/75‑1 (MQ=255) tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAtttt < 1:1390170/75‑1 (MQ=255) ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACtt > 2:947618/1‑75 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACttt > 2:275989/1‑75 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACttt < 2:902422/75‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtttt < 2:1035729/75‑1 (MQ=255) ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGAc < 2:954716/75‑1 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGa > 2:1472621/1‑75 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGa < 1:1441061/75‑1 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTc > 1:770061/1‑75 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTc < 2:29613/75‑1 (MQ=255) tgtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAATTTTTACGAttttt > 2:612340/1‑73 (MQ=255) tgtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCa < 2:307692/75‑1 (MQ=255) gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCAt < 1:644543/75‑1 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGt < 1:529575/75‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTCCCCTGTCGGTa > 2:1384916/1‑75 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTTTTTTTACGATTTTACCATGTCGGTAACGcc > 2:610996/1‑75 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTATATTTTTTTCGATTTCACCAAGTTGGTAAccac > 1:589211/1‑72 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTGTCGGGAACGCCg > 2:646980/1‑75 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTGTCGGGAAAGCCggg > 2:1287232/1‑74 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTTACGATTTTACCCTGGCGGTAACGCCGGa > 2:1538472/1‑75 (MQ=255) tCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTAAGATTTCACCCTGTCGGGAACGCCGGGAc > 1:865827/1‑75 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt < 1:1184645/75‑1 (MQ=255) gCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa < 2:280785/75‑1 (MQ=255) tCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg > 1:93298/1‑75 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGGAACGCCGGAACTGACGc > 1:989660/1‑75 (MQ=255) ctgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg < 2:35453/75‑1 (MQ=255) tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg < 1:1605367/75‑1 (MQ=255) tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg < 1:443439/75‑1 (MQ=255) gctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc > 2:189543/1‑75 (MQ=255) acacTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCCGGCGAat > 2:1343977/1‑75 (MQ=255) acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAatat > 2:955426/1‑75 (MQ=255) tGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATcc < 2:1129876/75‑1 (MQ=255) tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGGATATCCGTg > 2:1322827/1‑75 (MQ=255) aTGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc < 1:79020/75‑1 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt < 2:792861/75‑1 (MQ=255) gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt < 1:955426/75‑1 (MQ=255) gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 2:1116827/1‑75 (MQ=255) ggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt < 1:761946/75‑1 (MQ=255) | TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTT > REL606/1484031‑1484177 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |