Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 1,537,424 | 0 | C | G | 35.0% | 27.1 / 11.1 | 20 | A391G (GCG→GGG) | mdfA | Multidrug transporter MdfA |
Reads supporting (aligned to +/- strand): ref base C (13/0); new base G (0/7); total (13/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCC > ACB143/1537290‑1537551 | cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCg > 7:274959/1‑136 (MQ=255) aCAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCgggg > 3:169063/1‑136 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGc > 3:102443/1‑136 (MQ=255) gcgCACGGGGCTGGTTTCGGCCATGCTGGGTATGATCACGATTGCATTGCTGGGGTTTGGCGGAGCGCTGCGGGCGAGGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGggggcgggggcg < 8:301228/134‑1 (MQ=255) aaGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑gg > 5:34064/1‑136 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGcc > 7:31965/1‑136 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGCGGGTGGCGGGGGCGTT‑GGCGCTGTGccc > 3:332916/1‑136 (MQ=255) tGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTtgt > 1:118763/1‑136 (MQ=255) ggggATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTtgtg > 1:141236/1‑136 (MQ=255) atgatCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGtggtgg > 3:319500/1‑136 (MQ=255) gatCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTggtgggt > 5:256514/1‑136 (MQ=255) aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggc > 3:63667/1‑136 (MQ=255) gATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcgg > 6:54761/1‑136 (MQ=255) gctgctGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATaa > 1:200512/1‑136 (MQ=255) tggaaTTTGGCCGGGGCGGGGGGGGGGGGGGGGGGGTT‑GGCGTGGTGGCCGGGGTGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 7:112712/132‑1 (MQ=255) gtgcTTTTGGCCGGGGGGGGGGGGGTGGGGGGGGCGTG‑GGGGGTGGGGCGGCTGTGGTTTGTGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 4:257908/135‑1 (MQ=255) gggcttttgGCCCTGGGGGGGGGGGGGGGGGGGGCGGG‑GGGGCGGTGGCGGGTGGGGTTTGGGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 3:57505/131‑1 (MQ=255) gttttGCGGTGGGGGGGGGGGTGGGGGGGGGGTTGGGGGGTGGGGCCGGGGTGTTTTTGGGTGGGGGGGGGGCAGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAc < 8:208444/134‑1 (MQ=255) ttGCGCTGGGGGGGGGGTGGGGGGGGGCGTT‑GGGGCGGGGGCCGCGGGGGTTTGTGGTGGGTGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGc < 3:40519/135‑1 (MQ=255) ttGCCGGGGGGGGGGGGGGGGGGGGGGGGTTGGGCGGTGGGGCGGCTGTGGTTGGTGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg < 6:44899/135‑1 (MQ=255) gcgtGGGGGGGGGGGGGGGGGGGGGGTG‑GGGGGTGTGGCGGGGGTGGTTGGTGGTGGGTGGGGGGCGTGGGGAAGGGGCGGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt < 6:89861/133‑1 (MQ=25) gggggggggggggggggg‑ggCGGTGGGGCCGTGGGGGTGGGGGGGGGGGGGGGGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgcc < 7:237361/136‑1 (MQ=21) ggggggggggggggggTG‑GGCGCGGGGGCGGCTGTGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgcc < 1:2321/136‑1 (MQ=37) | CGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTT‑GGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCC > ACB143/1537290‑1537551 |
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: ‑ |