| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,544,492 | 0 | G | C | 32.0% | 33.0 / 19.3 | 25 | intergenic (+51/+63) | mdfA/ampG | multidrug MFS transporter/muropeptide permease AmpG |
| Reads supporting (aligned to +/- strand): ref base G (9/8); new base C (3/5); total (12/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG > NC_002947/1544352‑1544635 | cTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCggcgcggc > 7:20156/1‑147 (MQ=255) ggggggggggggGCGTTGGGGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAg < 8:20156/149‑1 (MQ=25) ggggggggggggGCGTGGGCGCTGGGGCCGCTGTGGTTGGTGGGGGGTGGGCGGCCTGGGGAACGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTa < 1:38711/149‑1 (MQ=21) gggggggggggCGTGGGCGCGGGGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGa < 6:80046/143‑1 (MQ=38) gggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGc < 4:196883/145‑1 (MQ=25) gggggggCGTGGGCGCTGTGGCCCCTGGGGTGTGGGGTGGGGGGCCGGCCTGGGGAAGGGTCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCgg < 7:63100/147‑1 (MQ=18) gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTa > 5:215112/1‑142 (MQ=255) gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTa < 6:215112/142‑1 (MQ=255) cggggggCGGTGGCGCTGTGGCCGCTGGGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCgg < 6:110393/147‑1 (MQ=21) gtgggCGCTGTGGCTGCTGTGGTTGGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTAccc < 3:36215/146‑1 (MQ=255) gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGGGT‑TAAACCCCTCCGACACAGGGCCGCGCCGCGTTTTTTGAGGGGGGTTTACCCCCGAA‑‑AGGccc > 3:121109/1‑148 (MQ=255) gcggcCTGGGTAAGGGGGTGTAGCTCTATAATCGGTATTGCCGGGGCGGGCCTCTTCGCGGTT‑AAACCCGCTCCTACACAGGGCCGCGCGGC‑TTTGTAGGAGGGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAgg < 4:121109/149‑1 (MQ=37) aaGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACTTAGCGGGT‑GAACCCGCTACAACACAGGGCCGCGCCGC‑TTTGTAGGCGCGGGGTTACACGCGAAC‑AGGCCAGCACAGGCAAAACAAACCATAAAGGCCGCTCAga > 6:228845/1‑147 (MQ=255) cggaaGCTCGAAAACCGGTATTGCCGGGGCCGGCCTCTTCGCGGGA‑AACCCCGCTCCACCACCGGCCCGCGCCGC‑TTTGTAGCAGCGGGTTAACCCGCGAAGCAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa < 1:236583/145‑1 (MQ=255) gccTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑GTTGTAGGAGCGGGGTTACCCGCGAAG‑AGGCCAGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCCCCGAGTGGTCGTCCGCCACCGGCGGCCcccgccccgc > 4:283138/1‑149 (MQ=25) tctTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc > 5:334017/1‑149 (MQ=255) tCGCGGGT‑AAACCCGCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt > 3:76475/1‑148 (MQ=21) ccGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCa < 2:17575/80‑1 (MQ=255) ccGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCa > 1:17575/1‑80 (MQ=255) gCTCCTACAAAGCGGCGCGGCCC‑TGTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg > 8:300695/1‑149 (MQ=21) tCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc > 2:124110/1‑148 (MQ=255) tCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc > 1:36938/1‑148 (MQ=255) cTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat > 5:130653/1‑149 (MQ=255) cacaGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag < 5:228845/149‑1 (MQ=255) cacaGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca < 6:334017/148‑1 (MQ=255) | CTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGT‑AAACCCGCTCCTACACAGGGCCGCGCCGC‑TTTGTAGGAGCGGGTTTACCCGCGAAG‑AGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG > NC_002947/1544352‑1544635 |
| 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: ‑ |