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 | 95.7% | 46.8 / ‑5.3 | 23 | intergenic (+152/+185) | codA/cynR | cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (13/9); minor base A (0/1); total (13/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.35e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAC > NC_000913/357472‑357731 | aGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGc > 1:723187/1‑139 (MQ=255) tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 1:306860/139‑1 (MQ=12) tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcgcg > 2:762755/1‑139 (MQ=12) cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAt < 1:1262453/139‑1 (MQ=18) agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAtt < 2:1286775/104‑1 (MQ=255) agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAtt > 1:1286775/1‑104 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg > 1:1233836/1‑128 (MQ=11) cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg < 2:1233836/128‑1 (MQ=1) cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc > 2:1416697/1‑139 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc > 1:65965/1‑139 (MQ=255) cgtaTTCCGGACGATGATCTGTATCAGCGCTTCCTGTGATTCATTACCGGCCACATCCGTGCATTTTTACAACGCCATGTGCAAAAAGGCCTTCGCAGCGGCTTTCTGCCAGCTGTAGTTTATGCCGTATGCGGCGTgc < 1:578074/136‑2 (MQ=12) aaTTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg > 1:251986/1‑57 (MQ=11) aaTTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 2:251986/57‑1 (MQ=11) ccGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg > 2:1219944/1‑49 (MQ=11) ccGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg < 1:1219944/49‑1 (MQ=14) gCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTATTTTATCCGCGACGCGGCGTGAACGCCTTATcgt > 1:1461056/1‑137 (MQ=14) tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCATCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa < 2:765947/139‑1 (MQ=18) tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa > 1:1067287/1‑139 (MQ=18) tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa > 1:1000195/1‑139 (MQ=18) gtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag < 1:589129/139‑1 (MQ=255) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg < 1:798362/107‑1 (MQ=255) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg > 2:798362/1‑107 (MQ=14) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCATTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc > 2:1407775/1‑139 (MQ=255) | AGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAC > NC_000913/357472‑357731 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |