Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,580,849 | T→A | 15.2% | intergenic (+27/‑56) | benR → / → PP_3160 | BenABC operon transcriptional activator/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,580,849 | 0 | T | A | 15.2% | 69.6 / 7.1 | 33 | intergenic (+27/‑56) | benR/PP_3160 | BenABC operon transcriptional activator/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (15/13); new base A (3/2); total (18/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01 |
CAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAG > NC_002947/3580708‑3580989 | cAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCa < 3:12399/149‑1 (MQ=255) tgtgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggt < 5:176196/149‑3 (MQ=255) tgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGc < 2:52595/149‑1 (MQ=255) tgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCAGTTATGCCAGTTGc > 2:61442/1‑149 (MQ=255) gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGtt < 4:208604/146‑1 (MQ=255) gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGtt > 3:208604/1‑146 (MQ=255) gTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAAccacca > 6:103145/1‑149 (MQ=255) cTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggtc < 4:214889/119‑4 (MQ=255) cTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggtc > 3:214889/1‑116 (MQ=255) tGGGGCGGTTTTCCGAGATTCATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt < 8:70119/127‑1 (MQ=255) tGGGGCGGTTTTCCGAGATTCATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt > 7:70119/1‑127 (MQ=255) gaTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg < 4:127225/144‑1 (MQ=255) gaTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg > 3:127225/1‑144 (MQ=255) aTCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCAt > 4:21722/1‑149 (MQ=255) gTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGtt < 7:19824/97‑1 (MQ=255) gTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGtt > 8:19824/1‑97 (MQ=255) tttGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCACTTGCGAGCGCTTAGCTGTTAAACAACCATTTGTTCTATGCCGGACGAATGGTCACAACCATGAAGCAGTCg > 7:234276/1‑148 (MQ=255) ttCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCCGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAat > 7:1657/1‑145 (MQ=255) ttCACAGACATTCAAGCCGAGCGCTTGACGGCGTTTGGCTTGTAGACCTGGAATTTATGCGAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCAGAGCCATGCAGCAGTCGCAGAAAGAAat < 8:1657/145‑1 (MQ=255) acaGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAatttattt < 3:21722/148‑1 (MQ=255) acaGACATTCAAGCGGCGAGAGTGACGGCGTTTGGCTTGTAGACCTAGCATTTATGACATTTGAGAGCGCTTTGCTGTTAAACCACCATTATTTCTATGCCGGACGAATGGTCACAGACATGCAGCAGTCGCAGAAACAAATTTATTTc > 1:195820/1‑149 (MQ=255) gACATTCAAGCGTCGCGGGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGGGCTTTGCTGTTAAACGACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTAc < 4:67958/149‑1 (MQ=255) gACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTAc < 5:103145/149‑1 (MQ=255) cgcgcgTGAAGGAGTTTGGCTTGTAGACCTGGCATTTATGCCTGTTGCGAGCGATTTGCTGTTAAACCACGATTTGTTCTATGCCGGACGAATGGTCACAACCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGAc < 8:234276/149‑1 (MQ=255) gCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCg < 3:17740/125‑1 (MQ=255) gCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCg > 4:17740/1‑125 (MQ=255) ttGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGtt > 1:25862/1‑149 (MQ=255) ctcgcaacTGGCATAAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCTTTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTg > 2:51349/8‑149 (MQ=255) tAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgc > 3:162149/1‑149 (MQ=255) aacTGGCATAAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgcg > 7:16496/3‑147 (MQ=255) aCCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGgcg > 6:263274/1‑147 (MQ=255) ccTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGGCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAg < 2:195820/149‑1 (MQ=255) tGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAg > 2:201081/1‑147 (MQ=255) | CAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGCGTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAG > NC_002947/3580708‑3580989 |
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: ‑ |