Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AM260480 | 580,260 | 0 | C | G | 22.7% | 46.2 / 11.6 | 22 | G69A (GGC→GCC) | h16_B0514 | transcriptional regulator, GntR‑family |
Reads supporting (aligned to +/- strand): ref base C (9/8); new base G (2/3); total (11/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCACGACAGGATTTCCGAGCGCGGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATATTTGCCTGCAGTTC > AM260480/580127‑580393 | gCACGACAGGATTTCCGAGCGCGGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGc < 2:925879/134‑1 (MQ=255) gCACGACAGGATTTCCGAGCGCGGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGc < 2:714956/134‑1 (MQ=255) aTTTCCGAGCGCGGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGAc > 2:368951/1‑134 (MQ=255) cgcgGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGGCCTGCTGGACCAGCACGGc < 1:623071/134‑1 (MQ=255) aaCCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGcc > 2:431608/1‑99 (MQ=255) aaCCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGcc < 1:431608/99‑1 (MQ=255) aaCCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCt > 2:239351/1‑134 (MQ=255) aTCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAg < 2:145504/134‑1 (MQ=255) ggCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAgcgcg < 1:368951/134‑1 (MQ=255) cGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCa < 2:627152/134‑1 (MQ=255) cgcACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGcc < 1:239351/134‑1 (MQ=255) aaGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGt < 2:878083/134‑1 (MQ=255) aGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGGCCTGCTGGACCAGCACGGCGATGGCCTGCCg < 1:609647/64‑1 (MQ=255) aGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGGCCTGCTGGACCAGCACGGCGATGGCCTGCCg > 2:609647/1‑64 (MQ=255) gTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCa > 2:392344/1‑134 (MQ=255) tgcccttgccctGCACGCGCTCCAGCAGGGCCTGCTGGAc < 1:674962/40‑1 (MQ=255) tgcccttgccctGCACGCGCTCCAGCAGGGCCTGCTGGAc > 2:674962/1‑40 (MQ=37) tGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGaaa > 2:50556/1‑133 (MQ=255) cACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACa < 1:392344/133‑1 (MQ=255) cgcTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATAtt > 1:850697/1‑134 (MQ=255) cgcTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATAtt > 2:540892/1‑134 (MQ=255) cagcagGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAg > 1:410885/1‑51 (MQ=255) cagcagGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAg < 2:410885/51‑1 (MQ=255) cagcagGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATATTTGCCt > 1:763602/1‑134 (MQ=255) cccTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATg > 1:978502/1‑51 (MQ=255) cccTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATg < 2:978502/51‑1 (MQ=255) cccTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATATTTGCCTGCAGTTc < 1:540892/134‑1 (MQ=255) | GCACGACAGGATTTCCGAGCGCGGCACCGCCGGCGCGCCATCGACGCCGCCGAAGCGGAAAAACCGCATCAGGCTGGCCCCGCTCAGTCCCTTGCGCACAAAGGTGCCCTTGCCCTGCACGCGCTCCAGCAGGCCCTGCTGGACCAGCACGGCGATGGCCTGCCGGATGGTGCCCAGCGCAATGCCGAACTCCCTGGCCAGCGCGCCTTCCGCAGGGATGGCGCTGCCGGGCGTCCACGCGCCCGAAACAATATTTGCCTGCAGTTC > AM260480/580127‑580393 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |