| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 164,569 | 0 | C | G | 41.9% | 25.2 / 19.1 | 43 | intergenic (+35/‑161) | hrpB/mrcB | RNA‑dependent NTPase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB |
| Reads supporting (aligned to +/- strand): ref base C (8/17); new base G (0/18); total (8/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.32e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGC > NC_000913/164429‑164705 | gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc < 1:1260726/141‑1 (MQ=255)gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGcc < 2:602190/141‑1 (MQ=255) aCCGTGAGGTGCAAACAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGACC‑gg < 2:272750/141‑1 (MQ=255) aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑gg < 2:2035141/141‑1 (MQ=255) cGTGAGGTGCAAAAAGTGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAATTTATTCGTAT‑GGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCC‑gggg < 2:230571/141‑1 (MQ=255) cGTGAGGTGAAAAAAGAGATGAATGGGCTTTATCAAAAACATGTCTGGCCGGACGACCCGGAAAAT‑ACTGCACCGACGCGACGGACTAAAAAGTATTCGTAT‑GTTTGGTTTTATCCCTCTCCCTGTGGGATTGGGCC‑gggg < 2:1311231/141‑1 (MQ=255) ttaacggggGTTAACCAAAACTTGTCTGGCCGTACGACCCGGAAAAT‑ACTGCACCGACGCGACGGACTAAAAATTATTCTTAT‑GTTTGTTTTTCTCCCTCTCCGTGGGGTTGTGGGCG‑GGGGTGCGGGCATCAGCGCGCAc < 2:125912/136‑1 (MQ=255) gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACg > 2:1212206/1‑141 (MQ=255) cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTggaccc > 1:1436443/1‑141 (MQ=255) tttaacacAAACTTGTCGGGCGGGACGCCCCGCAAAAT‑ACGCCACCGACGCGAGGGACTAAAAATTATTGTTAT‑TGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCG‑GGGGTGTGGGCATCAGCGCGCACGTtcaccct < 2:1791016/134‑1 (MQ=16) tttAACCCAAACTTGTCGGGCCGG‑CCCCCCGGGAAATAAAGGCCCCGCGGCTACGCACTAAAATTTTTTGGTAT‑GGTTGTTTTTCACCCTCACCCTGTGGTTGTGGGCG‑GGGGTGTGGCCCTCAGCGCGCACGTtcaccct < 2:2441601/140‑1 (MQ=255) gTTATCCAAAACATGTGTGGCCGGCCGACCCGGAAAAT‑ACTGCACCGACGCGAGGGACCAAAATTTATTCTTAT‑GGTTGGTTTACTCCCTCGCCCTGTGGGATTGGCCG‑GGGGTGCGGCCCTCAGCGCGCACGTtcaccct < 2:470628/141‑1 (MQ=255) ttATCCCACACAGGTCGGGCCGGCCCCCCCGCCAAAA‑CGTCCCCCCACGCGACGGACGAAAAAGTATACGTAC‑GTTTGGTTTTCGCCCCCCCCCTGGGGGAGCGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTtcaccctc < 2:885826/141‑1 (MQ=255) ttATCCCAAACTTGTCTGGCCGGACGACCCGGCAAAT‑ACTCCCCCGACGCGCCGGACTAAAAAGTAGTCTTAT‑GGTTGGTTTTCTCCCTCTCCGTGGGGTTGCGGGCG‑GGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:345024/141‑1 (MQ=255) ttATCCCAAACATGTCTGGCCGGACGACCCGGAAAAT‑ACTGCACCGCCGCGAGGGACGAAAAAGTATTCTTAT‑GTTTTGTTTTCTCCATCACCCTGTGGGTTAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:1017715/141‑1 (MQ=37) ttATCCCAAACAGGAGTGGCCGGACGACCCGCAAAAC‑ACTGCCCCGACGCGACGGACGAAATATTTTTGTTAG‑GGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCC‑GGGTTGCGGGCCTCAGCGCGCACGTtcaccctc < 2:1013706/141‑1 (MQ=255) ttATCCAAACAATGTGTGGCCGTACCACCCGGAAAAT‑CCGGCACGGACGCGCGGGAAGAAAATGTTTTCGTTA‑GTTTGTTTTTCTCCCTCTCCCTTTGGGAGTGGGGG‑GGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:1495/141‑1 (MQ=16) ttATCCAAAACTGGTCTGGCCGGACCCCCCGGAAAAT‑ACTGCACCGACGCGCCGGACTAAAATGTATTGGTAT‑GGTTGTTTTTCTCCCTCTCCGTGTGGGTGAGGGCG‑GGGGTGAGGGCCTCAGCGCGCACGTtcaccctc < 2:2058067/141‑1 (MQ=255) ttAACCCAAACTTTTGTGGCCGCACTACCGGGAAAAA‑AATGCACCGACGCGCGGGCCGAAAATGTATTGGTAT‑GGTTGTTTTTCTCCCTCGCCCTGTGGGAGCGGGCG‑GGGGTGTGGGCATCAGCGCGCACGTtcaccctc < 2:2481997/141‑1 (MQ=255) ttAACCCAAACTGGTCTGGCCGGACGCCCCGGCAAAA‑ACTGCACCGACGCGACGCACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGTACGTtcaccctc < 2:1213853/141‑1 (MQ=38) ttAACCCAAACATGTCTGGCCGGACGACCCGGAAAAT‑ACTGCACCGACGCGAGGGACCAAAAATTTTTCGAAT‑GGTTGGTTTACACCCTCACCCTGTGGGTGTGGGCG‑GGGGTTAGGGCATCAGCGCGCACGTtcaccctc < 2:2486300/141‑1 (MQ=255) ttAACCCAAACATGTCGGCCCGGACGACCCGGAAATA‑ACGCCACCGCGGCCACGCACTAAAATTTATTCTTAT‑GGTTGTTTTTCACCCACGCCCTGTGGGTGCGGGCG‑GGGGTGCGGGCATCAGCGCGCACGTtcaccctc < 2:2516183/141‑1 (MQ=255) ttAACCCAAAAATGTCTGGCCGGCCCACCCGGCAAAT‑ACTCCCCCGACGCTACGCACTAAAATTTATTCGTAA‑GGTTGGTTTACACCCTCCCCCGGTGGGAGAGGGCG‑GGGGTGAGGGCATCAGCGCGCACGTtcaccctc < 2:886503/141‑1 (MQ=25) tGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAcccaa > 1:657636/1‑139 (MQ=37) tGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCGAAcccgg > 2:176163/1‑139 (MQ=37) ttCTGGCCGGACGCCCCGCCAAAT‑ACGCCACCGACGCGACGGACGAAAAAGTTTTCGTAT‑GGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct < 2:96546/140‑1 (MQ=37) gTCGGGCCGGACGCCCCGGAAAAT‑ACGGCCCCGACGCGACGGACTAAAATGTTTTGGTAC‑GGTTGTTTTACGCCCTCTCCCTGGGGGAGCGGGCG‑GGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct < 2:2042934/141‑1 (MQ=255) cGGACTACCCGCCAACA‑CCTGCACCGACCCGACGCACTAAAATTTTTTCGTAT‑GGTTGGTTTTCCCCCACCCCCGGTGGGAGTGGGCG‑GGGTTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa < 1:1555009/141‑1 (MQ=255) cTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCGAACCCACAACCCCAAGGGGGGGGGGTCCggtgg > 2:1706738/1‑137 (MQ=21) ccgccgCGACGTCCGAAAACTTATCCGTAA‑GGTGGGTTTCCTCCCACTCCGTGGGGTTGTGGGCC‑GGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac < 2:326758/141‑1 (MQ=25) gacgcgacgGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGACCCACCCCCCAACCCACTCCCACAAAGGGGGGGGGGCCGGTagaggcaca > 2:107178/1‑136 (MQ=21) gacgcgacgGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGAACCAAACCCAAAACCACACCCGCCAGGGAGTGGGGCCCCAAAGAGGACaa > 2:107182/1‑141 (MQ=21) gacggacgAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGGTGCCCCGAACCCATAACCACACACACAAGGGGGGGGGGGCCGGACGGGCACAAAattt > 2:2670130/1‑141 (MQ=16) ggacgAAAAAGTATTCATAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag < 2:462012/141‑1 (MQ=255) ctaaaaaGTATACGTAA‑GGTTGGTTTTCTCCCACGCCCGGTGGGAGTGGGCC‑GGGGTGAGGGCATCACCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga < 2:1537513/139‑1 (MQ=39) ctaaaaGTTTTCCTAACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCG‑GGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga < 2:2511994/139‑1 (MQ=38) acgTATTCGTAA‑GGTTGGTTTACTCCCTCTCCCTGGGGGAGCGGGCG‑GGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct < 2:2138995/139‑1 (MQ=39) tCGTAT‑GGTTGTTTTTCGCCCTCTCCCGGTGGGAGTGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct < 2:708164/141‑1 (MQ=39) tcttatggttgtttTACTCCCTCACCCTGTGGGTGTGGGCG‑GGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct < 1:1317300/135‑1 (MQ=38) tcggatggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct < 2:73595/135‑1 (MQ=255) tcgaacggttggttTTCGCCCTCTCCCGGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct < 2:708106/135‑1 (MQ=255) ttggttTCCGCCCACCCCGTGTGGGTGGGGGCG‑GGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa < 1:1212206/141‑1 (MQ=38) gtggccc‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt < 2:2296166/136‑1 (MQ=255) gagGCCG‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt < 2:2296112/141‑1 (MQ=255) ggggCCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt < 1:294861/140‑1 (MQ=39) gggCC‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc < 2:2642658/141‑1 (MQ=255) | GTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAA‑GGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGC > NC_000913/164429‑164705 |
| 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: ‑ |