Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 357,606 | C→T | intergenic (+152/+185) | codA → / ← cynR | cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,606 | 0 | C | T | 97.0% | 77.3 / ‑6.4 | 33 | intergenic (+152/+185) | codA/cynR | cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (15/17); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCG > NC_000913/357481‑357734 | tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 1:691822/139‑1 (MQ=16) cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 2:840807/125‑1 (MQ=16) cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg > 1:840807/1‑125 (MQ=17) cTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcgcg > 1:618104/1‑118 (MQ=9) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 1:784826/81‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTa > 1:333154/1‑139 (MQ=11) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTCCGGCCAACAACcg > 2:784826/1‑81 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGtt > 2:253462/1‑139 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGttt > 2:907658/1‑139 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGGTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg > 2:748216/1‑125 (MQ=11) aGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGGTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg < 1:748216/125‑1 (MQ=11) cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg > 1:1288960/1‑139 (MQ=14) cATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg < 1:253462/139‑1 (MQ=255) cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggcg > 2:1031042/1‑139 (MQ=255) cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggcg > 2:1248855/1‑139 (MQ=14) gCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTAt < 2:327248/139‑1 (MQ=18) gCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTAt < 2:1311864/139‑1 (MQ=255) gCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg > 1:1291396/1‑139 (MQ=18) cTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATACCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg > 2:370012/1‑139 (MQ=18) tGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc < 1:907658/139‑1 (MQ=18) aTCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt < 2:333154/139‑1 (MQ=255) gtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag < 2:1288960/139‑1 (MQ=255) cAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag < 1:1248855/139‑1 (MQ=18) aGCGGTTACCGTCATTCATTTCCGCCAGCAACCGCGAATTTATCAAACGCCATGTGCGAAAATGCCTTCGCAGCGACTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTCATCCGGCCTACGTAGAGc < 2:1053102/139‑1 (MQ=25) gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg < 2:622641/107‑1 (MQ=14) gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg > 1:622641/1‑107 (MQ=14) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa < 1:636333/134‑1 (MQ=255) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa > 2:636333/1‑134 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCTGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa > 1:839662/1‑133 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg < 1:630488/101‑1 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg > 2:630488/1‑101 (MQ=14) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGcg < 2:839662/133‑2 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg < 1:1031042/139‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg < 1:1135981/139‑1 (MQ=255) | TATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCG > NC_000913/357481‑357734 |
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: ‑ |