Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,966,806 | C→A | 14.8% | L319M (CTG→ATG) ‡ | ydcP → | peptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,966,806 | 0 | C | A | 14.8% | 104.6 / 5.6 | 39 | L319M (CTG→ATG) ‡ | ydcP | peptidase |
Reads supporting (aligned to +/- strand): ref base C (16/17); new base A (2/4); total (18/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCG > NC_002947/3966664‑3966950 | gTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGt > 2:218562/1‑150 (MQ=255) gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCAGGCCGGTgaag < 1:510635/150‑4 (MQ=255) gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCAGGCCGGTgaag < 1:345803/150‑4 (MQ=255) ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCAGGCCGGTgca < 2:429181/147‑3 (MQ=255) ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCAGGCCGGTgaa > 1:429181/1‑145 (MQ=255) tCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGaa > 2:925/1‑150 (MQ=255) gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg < 1:925/150‑1 (MQ=255) gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg < 1:333447/150‑1 (MQ=255) aaCACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGCCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTg > 1:128731/1‑150 (MQ=255) acacTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGa > 1:126774/1‑150 (MQ=255) ttcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGt < 2:56212/144‑1 (MQ=255) ttcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGt > 1:56212/1‑144 (MQ=255) tcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg < 1:176632/150‑1 (MQ=255) ccTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGt > 1:269114/1‑150 (MQ=255) cccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTc < 1:486613/150‑1 (MQ=255) ccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAg > 1:95523/1‑122 (MQ=255) ccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAg < 2:95523/122‑1 (MQ=255) ggATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg > 2:94987/1‑137 (MQ=255) ggATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAg < 1:94987/137‑1 (MQ=255) aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCt > 1:78533/1‑150 (MQ=255) aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCt > 2:527282/1‑150 (MQ=255) aaCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACGCGCCAACCCACGCCACTCCAGCAAGGTTACGTTG‑GTAAAAGGTCGGCAAGCGTGACATGCAAATGGACACCGACGTACCGCt > 1:258407/1‑150 (MQ=255) ccTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGa < 1:218562/150‑1 (MQ=255) cttctcgGCAGCACCGTCTACTTTGTCACCGCACGCAAGTTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GTCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAAc < 2:60520/145‑1 (MQ=255) gACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAt < 2:269114/150‑1 (MQ=255) ctactTCGTGACCGACCGCAAGGCCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAtgt < 2:128731/150‑1 (MQ=255) ttCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCAGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg > 1:440622/1‑150 (MQ=255) gTGACCGACCGTAAGGTCGATATCGTCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTc < 1:561801/150‑1 (MQ=255) gaccgaccGCAAGGTCGATATCGGCGCCTTCGGCTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCaa < 1:527282/150‑1 (MQ=255) gaccgaccGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCaa > 2:121953/1‑150 (MQ=255) gTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc > 1:632692/1‑73 (MQ=255) gTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc < 2:632692/73‑1 (MQ=255) ccTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt < 2:25356/103‑1 (MQ=255) ccTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt > 1:25356/1‑103 (MQ=255) ccACC‑GGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCACAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACcgcgcg > 1:24700/2‑148 (MQ=255) aCC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt > 2:643751/1‑98 (MQ=255) aCC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCt < 1:643751/98‑1 (MQ=255) aCC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGc < 2:78533/150‑1 (MQ=255) cc‑GGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCACAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGCGGca < 2:24700/150‑2 (MQ=255) | GTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACC‑GGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCG > NC_002947/3966664‑3966950 |
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: ‑ |