Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,544,494 | 0 | C | G | 23.4% | 93.7 / 24.4 | 47 | intergenic (+53/+61) | mdfA/ampG | multidrug MFS transporter/muropeptide permease AmpG |
Reads supporting (aligned to +/- strand): ref base C (16/20); new base G (7/4); total (23/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.18e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGGGGGTGGCGGGGGCGTTGG‑CGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATC > NC_002947/1544357‑1544639 | ggcggGGGTGGGGGGGGCGTGGG‑CGCGGTGGCCGCGTTGTTGTGGGGTGGGGGGCCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑ttt < 6:137982/149‑1 (MQ=38) gggggTGGCGGGGGCGTTGG‑CGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTg > 7:35123/1‑147 (MQ=255) ggggggggggggggGTGGG‑CGGGGGGCCGGCGGTGTTTTGGGGTGGGGGGGCGTCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAg < 2:188534/149‑1 (MQ=12) gggggggCGTTGG‑CGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCgg < 3:181927/147‑1 (MQ=255) ggggggcCGTTGG‑CGCGGTGGCCGGTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGACCC‑TGTGTAGGAGCgg < 8:123901/147‑1 (MQ=25) gggggCGGTGGCCTCTGGGGCCGCTTTGGTTTGTGGTGGGTGGGCGGCCTGGGGATGTTGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCggg < 7:5544/149‑1 (MQ=37) cggtgGCCGCTGGGGTTTGGGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACACAGGGCCGCGCCGC‑TTTGTCGGAGCGGGTTTACCCGCGAaga < 8:35123/147‑1 (MQ=37) tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑ttt < 4:299821/123‑1 (MQ=255) tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑ttt > 3:299821/1‑123 (MQ=255) gTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGc > 6:394429/1‑118 (MQ=255) gTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGc < 5:394429/118‑1 (MQ=255) tgtgGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑ttt < 6:40892/107‑1 (MQ=255) tgtgGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑ttt > 5:40892/1‑107 (MQ=255) gggtgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGGGTAAAACCCCGCCGACACACGGCCGCGCCGC‑GTTGTAGGAGCGGGGTTTACCGCGGAGAGGCCAGCACACGGAATACAAAGc > 6:117956/1‑149 (MQ=255) gggtgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTATTCGCCGGTTAACCCCCTCCGACCCCGGGGCCCGCCGC‑GTTGTAGGGGCGGGGTTACCCGCCGAGGGGCCCGCACCCGCCATAaaaacc > 8:126256/1‑147 (MQ=255) gggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCCGGTAAACCCCCAGCGACACCCGGGCGCGCCGC‑GTTGTTGGGGCGGGGTTTACCGCGCAGGGGGCCGCCCAGGCAAAAAAAAACtt > 6:338010/1‑147 (MQ=255) cTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCgg < 2:188259/100‑1 (MQ=255) cTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCgg > 1:188259/1‑100 (MQ=255) gCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgt < 5:117956/148‑1 (MQ=255) cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg < 7:119033/109‑1 (MQ=255) tCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc < 3:304082/144‑1 (MQ=255) gCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTc < 3:52197/148‑1 (MQ=255) ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGG‑GTTGTAGGGGCGGGTTTACCCGGGAAGAGGCCCGCGCAGGCAATACAAAGCTTCAGGGCCGCTGAGCCACCGAGTGGTCGTCGGCCACCGGCGTCCCCCGCCCGGCTGCCt > 2:265956/1‑149 (MQ=37) ctgtaaaaCCCGCTCCAACACCGGGCGGCGCCGC‑TTTGTCGGCGGGGGTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATg < 5:22399/144‑1 (MQ=255) ggTAAACCCGCTCCTACACAGGGCCGCGCCGCATTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGc > 1:259124/1‑149 (MQ=255) aCCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 4:195442/1‑148 (MQ=255) aCCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg < 3:195442/148‑1 (MQ=255) aCCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 2:231843/1‑148 (MQ=255) cccGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcccg < 2:237975/130‑1 (MQ=25) cccGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcccg > 1:237975/1‑130 (MQ=25) ccGCTGCTTCACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACTAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTcc < 6:121889/135‑1 (MQ=255) ccGCTCCTACACAGGGCCGCGCCGC‑ATTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCGTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTcc > 5:121889/1‑135 (MQ=255) ccGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg > 5:120181/1‑146 (MQ=25) ccGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCa < 6:120181/148‑1 (MQ=37) cTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc < 1:265956/149‑1 (MQ=255) tCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc > 8:274339/1‑148 (MQ=21) ccTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg > 7:252826/1‑148 (MQ=37) cTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat > 5:78870/1‑149 (MQ=255) cTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg < 1:231843/147‑1 (MQ=255) cTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGCCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat < 5:338010/149‑1 (MQ=255) tACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat > 5:376774/1‑148 (MQ=25) caGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa < 7:312818/87‑1 (MQ=255) caGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa > 8:312818/1‑87 (MQ=255) aGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca > 6:32691/1‑145 (MQ=255) aGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca < 5:32691/145‑1 (MQ=255) ggCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggat > 3:350739/1‑147 (MQ=255) ggCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCAGATGCCACTGCAGGGAgctcaggatc > 8:394200/1‑148 (MQ=255) | GGCGGGGGTGGCGGGGGCGTTGG‑CGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATC > NC_002947/1544357‑1544639 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |