| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,691,844 | 0 | T | G | 27.8% | 67.9 / 10.6 | 36 | N226T (AAC→ACC) | yiaY | predicted Fe‑containing alcohol dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (8/18); new base G (10/0); total (18/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.44e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT > REL606/3691780‑3691918 | cTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATc < 2:778168/76‑1 (MQ=255) tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCa < 2:665713/76‑1 (MQ=255) tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCa > 2:17430/1‑76 (MQ=255) gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCAt < 2:571036/76‑1 (MQ=255) gCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATgg < 1:869582/76‑1 (MQ=255) cATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGt < 1:1018473/76‑1 (MQ=255) aaGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTTGGGAATCATGGTCAc > 2:1143571/1‑76 (MQ=255) aGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACg > 2:900272/1‑76 (MQ=255) cATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCt < 1:383914/76‑1 (MQ=255) ttGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCctt > 2:136890/1‑76 (MQ=255) gCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCa < 2:638772/76‑1 (MQ=255) cTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAg < 1:772837/76‑1 (MQ=255) ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCGATCATGGTCACGGCTTTTAgg > 1:1044120/1‑75 (MQ=255) ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCCTGGTCCCGGCTTTCAGt > 2:1143592/1‑76 (MQ=255) tCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTTGGGAAACATGGTCAAGGGTTTTAGTg > 2:863883/1‑76 (MQ=255) cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGc < 2:608298/76‑1 (MQ=255) gcgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGcaca > 1:435108/1‑76 (MQ=255) gcgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCAAGGCTTTCCGGGcaac > 1:435111/1‑74 (MQ=255) gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGCGcaccca > 1:535463/1‑73 (MQ=255) gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAg > 1:535481/1‑76 (MQ=255) tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAAACAAGGTCACGGCCTTCAGTGCACAAACg > 2:42139/1‑76 (MQ=255) gCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAAACAAGGTCAAGGCTTTCAGTGCAAAAACGTCaa > 1:517605/1‑75 (MQ=255) aTTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCCGTGCACCAGCCTCAGTg > 2:581039/1‑76 (MQ=255) aTTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTg > 2:581047/1‑76 (MQ=255) tACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGGGAt > 1:785311/1‑76 (MQ=255) cTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCg < 1:35294/76‑1 (MQ=255) tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCgg < 1:884156/76‑1 (MQ=255) aTCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTg < 2:183879/76‑1 (MQ=255) tCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTgg < 2:326053/76‑1 (MQ=255) cTTCAACGGCTAACGGCAGGGTTTCGGCAAACATGGTCACGGCTTTTCGTGCACAAGCGTCAGAGATCGGCGTggc > 2:1054375/1‑76 (MQ=255) ttCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggcg < 1:835516/76‑1 (MQ=255) ggCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATa < 2:532543/76‑1 (MQ=255) ggtAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg < 1:123375/74‑1 (MQ=255) gCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg < 1:123363/76‑1 (MQ=255) gCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg < 1:1076612/76‑1 (MQ=255) tAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGaa < 2:489936/76‑1 (MQ=255) aCGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAAc > 2:152042/1‑76 (MQ=255) cGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACa > 2:229147/1‑76 (MQ=255) aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACg > 2:687523/1‑76 (MQ=255) gTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCt > 2:287972/1‑76 (MQ=255) | CTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT > REL606/3691780‑3691918 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |