Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 977,876 | T→A | 13.3% | intergenic (+39/‑141) | cysE → / → iscR | serine acetyltransferase/DNA‑binding transcriptional regulator IscR |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 977,876 | 0 | T | A | 13.3% | 112.9 / 9.1 | 45 | intergenic (+39/‑141) | cysE/iscR | serine acetyltransferase/DNA‑binding transcriptional regulator IscR |
Reads supporting (aligned to +/- strand): ref base T (20/19); new base A (2/4); total (22/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01 |
ACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTAC > NC_002947/977733‑978007 | aCCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCAc < 1:470196/150‑1 (MQ=255) cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc > 2:843981/1‑150 (MQ=255) gATGGGCAGTGATTATTGTGCCAAGGAACTGCCGGCCTGGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGcccc < 2:797075/150‑1 (MQ=255) gATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGcccc < 2:295487/150‑1 (MQ=255) gTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGc < 2:533328/150‑1 (MQ=255) ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTa < 1:454324/140‑1 (MQ=255) ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTa > 2:454324/1‑140 (MQ=255) cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgccc > 2:756377/1‑150 (MQ=255) cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGcac > 1:313057/1‑120 (MQ=255) cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGcac < 2:313057/122‑3 (MQ=255) ccTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCa > 1:721282/1‑150 (MQ=255) cTTGCCGGAGGACGACTTCAGCGAGGCTACCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCaa > 1:23199/1‑150 (MQ=255) cGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCg > 2:728835/1‑150 (MQ=255) gaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGAc < 2:788230/150‑1 (MQ=255) gaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGAc > 1:670646/1‑150 (MQ=255) gacgacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaa > 2:295610/1‑150 (MQ=255) cgacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAg < 2:721282/150‑1 (MQ=255) gacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCCTTGAACGAGCCGAACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGc > 2:210387/1‑150 (MQ=255) cTTCCGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAg < 1:847229/150‑1 (MQ=255) ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTc > 1:614479/1‑103 (MQ=255) ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTc < 2:614479/103‑1 (MQ=255) ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc > 1:770734/1‑113 (MQ=255) ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc < 2:770734/113‑1 (MQ=255) gCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGt > 1:130772/1‑150 (MQ=255) tGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAg > 2:1193/1‑150 (MQ=255) cccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt < 1:843981/150‑1 (MQ=255) ccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTg > 1:560793/1‑150 (MQ=255) gTGGCCCAACGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACtt > 1:502667/1‑150 (MQ=255) cccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAAt > 1:326333/1‑150 (MQ=255) ccAGCGGGACACTCAGCCGCATTGATCGCGCGGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTATGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATg < 2:125839/150‑1 (MQ=255) cgcgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGTTTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg < 2:502667/150‑1 (MQ=255) cgcgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg < 1:295610/150‑1 (MQ=255) gATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCTCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGcaca < 2:23199/150‑1 (MQ=255) gccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAg < 1:45947/94‑1 (MQ=255) gccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAg > 2:45947/1‑94 (MQ=255) gccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCTCACACATCCt < 2:130772/150‑1 (MQ=255) tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTg < 1:728835/150‑1 (MQ=255) aTGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGt > 2:226550/1‑150 (MQ=255) ggtggCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTAc < 1:1193/150‑1 (MQ=255) gCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTaa < 2:173684/86‑1 (MQ=255) gCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTaa > 1:173684/1‑86 (MQ=255) ggggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa < 2:112818/87‑1 (MQ=255) ggggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa > 1:112818/1‑87 (MQ=255) ggggCGTGAGCACACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa < 2:463442/62‑1 (MQ=38) ggggCGTGAGCACACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa > 1:463442/1‑62 (MQ=38) | ACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTAC > NC_002947/977733‑978007 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |