Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,552,005 | G→A | 14.3% | intergenic (+55/+126) | pydA → / ← rutR | NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,552,005 | 0 | G | A | 14.3% | 61.4 / 4.3 | 28 | intergenic (+55/+126) | pydA/rutR | NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family |
Reads supporting (aligned to +/- strand): ref base G (10/14); new base A (2/2); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01 |
TGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCC > NC_002947/4551864‑4552150 | tGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCggg > 1:389411/1‑150 (MQ=255) cATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGcc < 2:760569/150‑1 (MQ=255) aTCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCg > 2:451364/1‑150 (MQ=255) aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc < 1:405279/150‑1 (MQ=255) aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc < 1:47584/150‑1 (MQ=255) cacaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCca < 1:726917/150‑1 (MQ=255) aGGATACCGGCAAGCCGTACCTGAACGGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGTGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca < 2:360646/150‑1 (MQ=255) gATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGt < 2:659939/150‑1 (MQ=255) aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc < 2:389411/150‑1 (MQ=255) gCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGcc > 2:570091/1‑150 (MQ=255) tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa < 1:570091/150‑1 (MQ=255) gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCgg < 1:451364/150‑1 (MQ=255) caaCCCCTACCGTGAGGCCAGTTGAAGCTGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc < 2:694265/149‑1 (MQ=255) cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCg < 1:548985/135‑1 (MQ=255) cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCg > 2:548985/1‑135 (MQ=255) gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTcg < 1:538663/69‑1 (MQ=255) gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTcg > 2:538663/1‑69 (MQ=255) gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCTCTTCGCGGCCCGATCGCGac > 2:241024/1‑148 (MQ=255) gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCTCTTCGCGGCCCGATCGCGac < 1:241024/148‑1 (MQ=255) tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc < 2:388649/85‑1 (MQ=255) tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc < 2:388727/85‑1 (MQ=255) tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc > 1:388727/1‑85 (MQ=255) tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc > 1:388649/1‑85 (MQ=255) ttGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCt > 2:494376/1‑149 (MQ=255) ttGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCt < 1:494376/149‑1 (MQ=255) tGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg < 1:204169/94‑1 (MQ=255) tGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg > 2:204169/1‑94 (MQ=255) ccTTGTGTCGCGATCGGGCCGCACAGCGGCCCCACCGGTCTTTGCCCCGATACCGACATCGAGGGGCCAGGTCGCGGCCCGATCGCCACCCAACGCCGCTCCCACAAAAAGACCACCTCCACCTCGAATGCCGGATCACGGCTCCCGccc > 1:748043/1‑150 (MQ=255) ttGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCAGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGCCTCCAGCCCGa > 2:368035/1‑150 (MQ=255) tgtCGCGATCGGGCCGCACATCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACCCAAGGCCGCTCCCACAAAAAGACCTCGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcc > 2:315476/1‑150 (MQ=255) | TGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCC > NC_002947/4551864‑4552150 |
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: ‑ |