Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,965,101 | 0 | T | G | 30.0% | 65.6 / 10.5 | 40 | N19H (AAC→CAC) | thyA | thymidylate synthase |
Reads supporting (aligned to +/- strand): ref base T (9/19); new base G (12/0); total (21/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.16e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGTCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTG > NC_000913/2964956‑2965248 | gtgtCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGttttt < 1:1389179/149‑1 (MQ=255) gCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCtgtg < 2:1130341/149‑1 (MQ=255) gCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCtgtg < 2:702188/149‑1 (MQ=255) gAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGGCttag > 2:1730746/1‑147 (MQ=255) cagcagTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCa < 1:1085152/149‑1 (MQ=255) cagcagTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCa < 1:1086160/149‑1 (MQ=255) agTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCACAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCAATCGTATAGCaaatt > 2:1107966/1‑145 (MQ=255) tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTg < 2:293494/149‑1 (MQ=255) aTGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCa < 1:744715/149‑1 (MQ=255) atgatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAg > 2:692350/1‑149 (MQ=255) gatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGtt < 1:1032142/149‑1 (MQ=255) atgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGACAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTc > 2:1685952/1‑149 (MQ=255) gAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGGCGTCGTGGGGGACTTTTTTGATCAGTTGTaa > 1:892219/1‑149 (MQ=255) aCGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTTTTCCTTCGTCGTGCACTTTTTTCATAAGTTGTAAAt > 1:2100482/1‑149 (MQ=255) cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTTTGTCTTGGTCGAGCACATTTTTAATAAATTATAAATACTGTTTCaa > 1:1276057/1‑148 (MQ=255) cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGGCATAGTCGAGCACGTTTTTCATAAGGTGTAAATAATTTTTAAt > 2:700557/1‑149 (MQ=255) cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTATGTTTTCCTTCGTGGAGCACGTTTTGAATAAGTTATAAATAATTTTTAAt > 2:699991/1‑149 (MQ=255) aCGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATTTGTGGCATACTCGACCAAATTTTGGATCAGGTATAAATACTGGTTAAttg > 2:1017186/1‑147 (MQ=255) tGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAgg < 1:1426583/149‑1 (MQ=255) cACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATTGTTCCTCAGGaca > 2:2061962/1‑147 (MQ=255) cACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATGTGTGCCATCGGCGAGCACTTTTTGCATCAGTTCTAAATAAAGTTTAAAGGTTACTCAGGaaa > 2:1006192/1‑149 (MQ=255) aGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACgtg < 1:537667/149‑1 (MQ=255) gggAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTTGCTTAGTCGACCACGTTTTTGAAAAGTTATAAATAAATTTTAATGGTTGCTCAAGAAACATGTTg > 2:240335/1‑149 (MQ=255) tCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCtgtg < 1:1274478/149‑1 (MQ=255) aGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTTTGCCTTCGTGGAGCCCGTTTTTCATAAGTTCTAAATACTGTTTGAAGTTTACTCAGGTAACGTGTTTCGGTGGGGTggggcg > 1:1067811/1‑144 (MQ=255) ggTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTTTGTGTGGCGTCGTCGGGCACATTTTTGATAAGTTGTAAATAATTTTTAAATGTTGCTAAGGGAAACTGTTTGTGTGGGGGGcggcga > 1:763935/1‑149 (MQ=255) ggTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTGATGGTTGCTAAGGAAACGTGTTGCGGTGGGCTGcgacga > 2:523996/1‑149 (MQ=255) ttAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATTTGTGCCTTAGTCGAGCACTTTTTTGATCAGTTCTAAATAATGTTTGATTGTTACTAAGGAAACGTGTTGCTGTGGGGTGCGACGata > 2:96299/1‑149 (MQ=255) ttAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGCCGTAGTGGTGCACAGTTTTGAACAGGTCTTAAAAATGTTTGATTGTTGCTGAGGAAACGTGTTTGTGTGGGGGGCGGCGata > 2:785203/1‑149 (MQ=255) tAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATAAATTTTGAAGGTTCCTCAGGAAACGTGTTGCGGTGGGCGGCGACGatat > 2:779426/1‑149 (MQ=255) aCGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGccc < 2:993512/149‑1 (MQ=255) gatgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACcatc < 2:538300/149‑1 (MQ=255) tgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACcatcat < 1:1170552/149‑1 (MQ=255) ccAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt < 1:699991/149‑1 (MQ=255) ccAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt < 1:700557/149‑1 (MQ=255) cAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATc < 1:565139/149‑1 (MQ=255) aaaTGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATcaca < 1:303801/149‑1 (MQ=255) ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCg < 2:2100482/149‑1 (MQ=255) tcaggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGa < 1:268779/146‑1 (MQ=255) cGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAAtt > 1:1947129/1‑149 (MQ=255) cGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAAtt > 2:1852933/1‑149 (MQ=255) tttttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTg < 2:1837559/149‑1 (MQ=255) | GTGTCGCCCTGCAGAAACCACAGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTG > NC_000913/2964956‑2965248 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |