| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,132,851 | 0 | A | G | 44.4% | 41.9 / 12.3 | 36 | N65N (AAT→AAC) | yceI | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (2/18); new base G (16/0); total (18/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTG > REL606/1132782‑1132925 | aTCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGtttttt > 2:1383475/1‑76 (MQ=255) ttCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCa < 1:655485/76‑1 (MQ=255) gTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTa < 1:686071/76‑1 (MQ=255) tGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaa < 2:355189/76‑1 (MQ=255) tGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaa < 2:1265251/76‑1 (MQ=255) gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTTTTAAAGGGaaa > 1:322397/1‑76 (MQ=255) gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGGTTTTTTTGTCAAAGGGaaa > 2:579597/1‑76 (MQ=255) gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaaa < 1:1160445/76‑1 (MQ=255) gATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTaaa < 2:1226075/76‑1 (MQ=255) tCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGt < 1:420859/76‑1 (MQ=255) gACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGGCGGGTTTTTTTTGTTAAAGGGAAAAATACCGTCg > 2:1358809/1‑76 (MQ=255) aCGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGa < 1:1165/76‑1 (MQ=255) cTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAGAAACGTCGAAAt > 2:377210/1‑76 (MQ=255) cTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGGAAAAAAACCGTCGAAAt > 2:452009/1‑76 (MQ=255) gtggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGGAAAAAAACCGTCGAAATC‑tt > 2:1422736/1‑76 (MQ=255) gtggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGGAAAAAAACCGTCGAAAAC‑tt > 2:416137/1‑76 (MQ=255) gtggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑tt < 1:561324/76‑1 (MQ=255) gtggtgTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑tt < 2:872992/76‑1 (MQ=255) gtgtTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTaa < 1:469952/76‑1 (MQ=255) tAATTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGGAAAAAAAACGTCGAAAACTTTTAAAgg > 2:813499/1‑76 (MQ=255) tAATTGTCACATTCACTTTATCGGCAGCCGGGGTTTTTTCGTCAAAGGGAAAAAAACCGCCGCAATC‑TTTAAAGGt > 2:467400/1‑76 (MQ=255) aaTTGTCACATTCACTTTATCGGCAGCCGGGCTTTTTTTGTCAAAGGTAAAAGTAACGTCGAAAAC‑TTTAAAAGTg > 2:635281/1‑76 (MQ=255) aaTTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTg < 2:78795/76‑1 (MQ=255) aaTTGTCACATTCACTTTATCGGCAGCCCGGTTTTTTTTTTCAAAAGGAAAAAAACCGTCGAAAACTTTTAAAggg > 1:1191378/1‑75 (MQ=255) aTTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTTTCAAAAGGAAAAAAACCGTCGAAAAC‑TTTTAAAGTGc > 1:256460/1‑76 (MQ=255) aTTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTTAAAGGGAAAAAAACCGTCGAAAAC‑TTTAAAGGAGc > 1:152574/1‑76 (MQ=255) aTTGTCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAAAACCGTCGAAATA‑TTTTAAGGGGc > 2:207730/1‑76 (MQ=255) tCACATTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAGGTAAAAATACCGTCGAAATC‑TTTTAAGGTGCCGta > 2:409516/1‑76 (MQ=255) tCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGGAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGta > 2:1475987/1‑76 (MQ=255) cacaTTCACTTTATCGGCAGCCGGGTTTTTTTTGTCAAAAGGAAAAAAACCGTCGAAATC‑TTTTAAGGTGGCGgat > 1:1449077/1‑76 (MQ=255) ttCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAAcc < 2:579160/76‑1 (MQ=255) cccTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAGCCAg < 2:1080063/74‑1 (MQ=255) tttATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCta < 2:1269184/76‑1 (MQ=255) tCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGc > 2:415902/1‑76 (MQ=255) gCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGt < 1:943054/76‑1 (MQ=255) cGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGc < 1:334745/76‑1 (MQ=255) gATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTg < 2:952091/76‑1 (MQ=255) | ATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGCCGGATTTTTTTCGTCAAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTGCTG > REL606/1132782‑1132925 |
| 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: ‑ |