Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,665,195:1 | (T)7→8 | 100% | intergenic (+75/‑120) | dmsD → / → clcB | twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,665,188 | 1 | . | T | 89.4% | 139.1 / 11.8 | 47 | intergenic (+68/‑127) | dmsD/clcB | twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter |
Reads supporting (aligned to +/- strand): ref base . (3/2); new base T (19/23); total (22/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG > NC_000913/1665089‑1665281 | ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGGTTACGCGCAATACGGCAATAATAGGTTAAAGTGTCAC‑g > 1:2366344/1‑100 (MQ=255) ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑t > 2:264880/1‑100 (MQ=255) gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttt > 1:1330021/1‑100 (MQ=255) tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑ttttt < 2:1500892/100‑1 (MQ=255) cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt < 1:1036426/100‑1 (MQ=255) cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt > 1:462145/1‑100 (MQ=255) cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt < 2:338569/100‑1 (MQ=255) gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttt < 1:1974232/99‑1 (MQ=255) cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑ttttttta < 2:748621/100‑3 (MQ=255) ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttttat < 1:141041/100‑4 (MQ=255) ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑t‑tttttttatct < 2:1127938/100‑6 (MQ=255) ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTaaa < 2:2065755/100‑1 (MQ=255) tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca < 2:1168198/100‑1 (MQ=255) tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca < 2:33704/100‑1 (MQ=255) ntttCGATAGGATCACCGTAATATTGCCGGAAGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGAGCCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgca > 2:2208132/2‑99 (MQ=255) gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcac > 1:394868/1‑100 (MQ=255) tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacg < 1:1246962/100‑1 (MQ=255) tttCGATAGGATCACCGTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacg < 2:1333412/100‑1 (MQ=255) ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAgcacgc > 1:288755/1‑100 (MQ=255) gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GTTTT‑TTTTTTATCTCTTAAAgcacgc > 1:821895/1‑99 (MQ=255) ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt > 1:1697237/1‑100 (MQ=255) gcaTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt > 1:1624358/3‑100 (MQ=255) aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACTCACTGCtttt > 1:865884/1‑100 (MQ=255) tCACCGTAATTTTGCCGGATGGTTACGATTTCGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCttt < 1:163375/98‑1 (MQ=255) aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg < 1:2382925/100‑1 (MQ=255) ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg > 2:762420/1‑100 (MQ=255) cGTAATATTGCCGGAGGGTGTTGATTACGCGCACTCTTTCAATAATAGGTTACATTGTCTCTT‑‑‑T‑TTTTTTATCTCTTTAATCACGCAGTGCTTTTg < 2:1036092/96‑1 (MQ=255) tcttGCCTGCTGTTTACGATTACGCGCAATCCTTCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACTTACTTCTTTTGCGGCTGGcc < 2:884786/98‑1 (MQ=255) tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc < 2:621358/100‑1 (MQ=255) tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc < 1:400636/99‑1 (MQ=255) ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc < 2:311680/100‑1 (MQ=255) tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGAACGCACTGCTTTTGCGGCTGGCctct > 1:1838646/1‑100 (MQ=255) gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTTTTTCTTTTATCACGCACTGCTTTTGCGGCt > 2:1437973/1‑92 (MQ=255) cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt < 2:1236185/99‑1 (MQ=255) aTGGTGACGACTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAATCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc < 2:1573989/100‑1 (MQ=255) gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTACGGCTGGCCTCTTTTGACGCa > 1:370572/1‑100 (MQ=255) gTGACGATTACGCGCAATCCGGCAATAATAGGATACAGTGGCACGG‑‑‑T‑TTTTGTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa > 1:2281700/1‑100 (MQ=255) gTGACGATTACGCGCAATCAGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa < 2:650361/100‑1 (MQ=255) cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa > 1:1213753/1‑99 (MQ=255) aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt < 1:1143039/98‑1 (MQ=255) tACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc < 1:264883/100‑1 (MQ=255) aCGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg < 1:1555146/100‑1 (MQ=255) aaTCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATATCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTAGCCCgtg > 2:224061/1‑100 (MQ=255) aTCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt > 2:1724079/1‑100 (MQ=255) aTCCGCCAATAATAGGTTACAGTCCCAC‑G‑‑‑TCTATCTTAACTCTTAAAGCACCCACTCCTCATGCGGCTAGCCTCTTTTGCCGCCACATatccgc > 2:2223249/1‑89 (MQ=255) tCCGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt < 1:413462/100‑1 (MQ=255) cGGCAATAATAGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc > 1:1246367/1‑100 (MQ=255) ataataGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa < 1:221997/96‑1 (MQ=255) taataGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc < 2:288752/100‑1 (MQ=255) aataGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa > 2:401270/1‑100 (MQ=255) aGGTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGTCTCTTTTGCTGCAAAATAGTCGCCCGTGTTTCATTGCCCAtt > 2:1686521/1‑100 (MQ=255) ggTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt > 2:1890481/1‑100 (MQ=255) gTTACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc > 2:806650/1‑100 (MQ=255) ttACTGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt < 1:35568/100‑1 (MQ=255) ttACAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt > 1:1909645/1‑100 (MQ=255) aCAGTGTCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc < 1:289537/100‑1 (MQ=255) gtCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg < 1:1681953/100‑1 (MQ=255) tCAC‑GT‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg < 1:617683/99‑1 (MQ=255) | TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAC‑G‑‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG > NC_000913/1665089‑1665281 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |