Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,638,457 | C→T | 12.1% | E109K (GAA→AAA) | PP_1437 ← | heavy metal sensor histidine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,638,457 | 0 | C | T | 12.1% | 77.7 / 3.6 | 33 | E109K (GAA→AAA) | PP_1437 | heavy metal sensor histidine kinase |
Reads supporting (aligned to +/- strand): ref base C (13/16); new base T (2/2); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01 |
GTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGGCGCGGCTGCCCAGGTCCTCGCTGCGCAG > NC_002947/1638320‑1638604 | gtcgGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGa > 1:42292/1‑150 (MQ=255) gAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAgggccgg < 1:394643/150‑1 (MQ=255) gACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGAGAAAGGCCAGGCGCGCATCCGTGTGCAGTTCACGGACCCGCGATTCCAtgggcgcgcaa > 1:368637/1‑148 (MQ=255) aCATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTGACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGa < 2:553275/150‑1 (MQ=255) ccTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCa < 2:352329/150‑1 (MQ=255) aCAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTgg < 2:335111/150‑1 (MQ=255) aGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTg > 1:218275/1‑139 (MQ=255) aGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTg < 2:218275/139‑1 (MQ=255) gcgcATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGc > 1:440782/1‑150 (MQ=255) cATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAg > 1:240603/1‑150 (MQ=255) gTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAgg < 1:203576/150‑1 (MQ=255) gTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAgg < 2:440782/150‑1 (MQ=255) ggTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAg < 1:87929/150‑1 (MQ=255) ttACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGcccc < 2:489335/150‑1 (MQ=255) ccgccgCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATc > 1:231639/1‑150 (MQ=255) gccgcTGTCACGCGAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCAc > 2:38641/1‑150 (MQ=255) tttCACGCCTGGCGTGACAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGt > 1:51808/3‑150 (MQ=255) tGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGcaca < 2:57781/116‑1 (MQ=255) tGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGcaca > 1:57781/1‑116 (MQ=255) tCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTTACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTcc > 2:335126/1‑150 (MQ=255) tCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTTACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTcc > 2:335278/1‑150 (MQ=255) cgcGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTgggg < 1:420781/86‑1 (MQ=255) cgcGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTgggg > 2:420781/3‑88 (MQ=255) aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCACGGGGGGGTGGCGGCCGTCGAGATCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCcag > 1:168210/1‑150 (MQ=255) gCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGa < 2:42292/150‑1 (MQ=255) gTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATg < 2:231639/150‑1 (MQ=255) tGAAAGGCCAGGCGCGCATCGGTGTGCAGTTTACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATgg < 1:335278/150‑1 (MQ=255) tGAAAGGCCAGGCGCGCATCGGTGTGCAGTTTACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATgg < 1:335126/150‑1 (MQ=255) tgtgCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGGCGCGGCTGCCCAGGTCCTCGc > 2:45614/1‑150 (MQ=255) gtgCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGGCGCGGCTGCCCAGGTCCTCGCt < 2:240603/150‑1 (MQ=255) gTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAg < 2:549904/79‑1 (MQ=255) gTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAg > 1:549904/1‑79 (MQ=255) ttCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGGCGCGGCTGCCCAGGTCCTCGCTGCGCAg < 1:45614/150‑1 (MQ=255) | GTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGGCGCGGCTGCCCAGGTCCTCGCTGCGCAG > NC_002947/1638320‑1638604 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |