Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,575,074 | 0 | A | G | 68.2% | 21.6 / 18.8 | 22 | intergenic (‑27/+25) | ydeM/ydeN | anaerobic sulfatase maturase/sulfatase |
Reads supporting (aligned to +/- strand): ref base A (3/4); new base G (7/8); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTC > NZ_CP009273/1574938‑1575209 | cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGaa > 1:318921/1‑138 (MQ=255) cGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATttcgct > 2:195091/1‑139 (MQ=255) aCTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATttcgcttc > 1:182542/1‑139 (MQ=255) cAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGATGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTg < 2:225413/138‑1 (MQ=255) aaGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGc < 1:95400/138‑1 (MQ=255) tcAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGata < 2:563221/138‑1 (MQ=255) tAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGt > 1:463226/1‑139 (MQ=255) gagaTTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTtctc < 1:264137/138‑1 (MQ=255) tGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGa < 1:503712/68‑1 (MQ=255) tGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGa > 2:503712/1‑68 (MQ=255) gAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTtcgcttcgct > 1:403739/1‑89 (MQ=255) gAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTtcgcttcgct < 2:403739/89‑1 (MQ=255) gAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAg < 2:463226/137‑1 (MQ=255) tAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGgctgc < 2:301740/138‑1 (MQ=255) tCGCTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATaaa < 1:1200/139‑1 (MQ=255) cTATCTCGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAActc > 1:57411/1‑139 (MQ=255) tctcGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAActctct > 2:541419/1‑139 (MQ=255) cGCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTAc > 2:184017/1‑139 (MQ=255) gCTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTAcc < 2:79067/139‑1 (MQ=255) cTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCgg < 2:244408/107‑1 (MQ=255) cTCAATAAGGCGGCGGGAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCgg > 1:244408/1‑107 (MQ=255) aaTAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACg < 1:541419/138‑1 (MQ=255) cggGAAAATCCGCCGCATGAAGGTTTAGTTATATCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCAtctc > 1:495697/1‑139 (MQ=255) | CTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTC > NZ_CP009273/1574938‑1575209 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |