| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 164,569 | 0 | C | G | 52.9% | 0.9 / 15.0 | 17 | intergenic (+35/‑161) | hrpB/mrcB | RNA‑dependent NTPase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB |
| Reads supporting (aligned to +/- strand): ref base C (4/4); new base G (0/9); total (4/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.94e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCT > NC_000913/164434‑164669 | gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt > 1:961504/1‑141 (MQ=255)gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt > 1:1316632/1‑141 (MQ=255) cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTTTTCGTGAGGTG‑GGGTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTggaccc > 2:1300388/1‑141 (MQ=37) cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc > 2:1300402/1‑141 (MQ=255) gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTTTTCGTAAGGTT‑GTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct < 2:588571/141‑1 (MQ=38) ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGTTT‑GGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:950349/141‑1 (MQ=255) ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTCCACCGACGCGACGGACGAAAAATTATTCGTAAGGTT‑GGTTTGCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtcaccctc < 2:887041/141‑1 (MQ=38) ttATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTCCACCGACGCGCCGGACTAAAAATTATTCTTAGGGTT‑GGTTTACTCCCACTCCCGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtcaccctc < 2:775599/141‑1 (MQ=255) ttATCCAAACCAGGGCGGGCCGGACGACCCGGAAAATCGGGCCCCGCCGCGCCGGACG‑AAAACTTTTTGTTAGGTTGGGTTTCCGCCCCCCCCCGGGGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtcaccctc < 2:1156203/141‑1 (MQ=255) ttATCCAAAACAGGTGGGGCCGGACGACCCGGCAAATACTGCCCCGCCGCGACGGACGAAAAATTATACGTACGGTT‑GGTTTTCTCCCTCACCCTGTGGGAGTGGGCGGGGTTGTGGCCCTCAGCGCGCACGTtcaccctc < 2:1904971/141‑1 (MQ=255) taaccccAAACTTGTGGGGCCGCACGACCGGGCAAATCCTCCACCGCCGCGACGTACCAAAATGTTTTCTTAAGGTG‑GGTTTACTCCAACGCCGGGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:775502/137‑1 (MQ=255) taacccaaaaCATGTCTGGCCGGCCGACCCGGAAAATCCTCAACCGACGCGAGGTACTAAAAATTATTGTTATGGTT‑TTTTTTCTCCCTCACCGTGGGGTTGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTtcaccctc < 2:245097/134‑1 (MQ=255) aacccAAACATTTCTGGCGGGACCACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTATGGTT‑GTTTTTCGCCCTCTCCGTGTGGGTGAGGGCCGGGTTGCGGGCATCAGCGCGCACGTtcaccctcac < 2:643296/139‑1 (MQ=255) ccGACGCGACGTACGAAAAAGTATTCTTAAGGTT‑GTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac < 2:563351/141‑1 (MQ=38) aaaacGTATTCGAATGGTG‑GTTTTTCTCCCACACCCGGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata < 2:961504/141‑1 (MQ=38) tCGTATGTTT‑GTTTTTCTCCCGCCCCCTGTGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATAtcttcttct < 2:1322635/141‑1 (MQ=38) tCGTAAGGTT‑GGTTTTCACCCTCTCCCGGGGGGAGAGGGCGGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct < 2:247202/141‑1 (MQ=39) | GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTT‑GGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCT > NC_000913/164434‑164669 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |