Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,526,631 | A→C | 10.3% | E259A (GAA→GCA) | ddlB → | D‑alanine‑‑D‑alanine ligase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,526,631 | 0 | A | C | 10.3% | 103.2 / 4.9 | 39 | E259A (GAA→GCA) | ddlB | D‑alanine‑‑D‑alanine ligase B |
Reads supporting (aligned to +/- strand): ref base A (20/15); new base C (2/2); total (22/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
CCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGC > NC_002947/1526485‑1526774 | ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCCCCCAGCCTTCGATGCCAATGGCATCGCAg < 2:308647/150‑1 (MQ=255) tGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCt > 2:247873/1‑150 (MQ=255) ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTgg > 1:158841/1‑150 (MQ=255) gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg > 1:203199/1‑150 (MQ=255) gtgtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGAc < 1:247873/150‑1 (MQ=255) ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCACGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg > 2:56005/1‑150 (MQ=255) aCGACTATGACGCCAAGTACATCTTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGc > 2:346259/1‑150 (MQ=255) cTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGTGGTTGGACGTGATGCAGGa < 2:472279/150‑1 (MQ=255) tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc < 2:292193/150‑1 (MQ=255) cAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg > 2:493016/1‑150 (MQ=255) cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGtt > 2:236764/1‑150 (MQ=255) gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg > 2:579152/1‑150 (MQ=255) aaCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGCTCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg > 2:426429/1‑150 (MQ=255) cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGAGGTGATGCAGGACGAGCAGGGCCGGTTCTGgctgct > 2:330449/1‑150 (MQ=255) aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGCCCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa < 1:485491/150‑1 (MQ=255) ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt > 2:410957/1‑150 (MQ=255) aCCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAaca < 1:579152/150‑1 (MQ=255) aCCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAaca < 1:236764/150‑1 (MQ=255) cGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACAcc > 2:64927/1‑150 (MQ=255) cATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGc > 2:555665/1‑150 (MQ=255) cTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACCCCGCACCAGGCATGACa < 2:269522/150‑1 (MQ=255) aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGTCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATc < 1:449250/150‑1 (MQ=255) aacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGc < 1:555665/150‑1 (MQ=255) aggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt < 2:307156/127‑1 (MQ=255) aggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt > 1:307156/1‑127 (MQ=255) gACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc < 2:432320/88‑1 (MQ=255) gACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc > 1:432320/1‑88 (MQ=255) gACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGt > 1:585120/1‑150 (MQ=255) aCCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGt > 1:564010/1‑51 (MQ=255) aCCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGt < 2:564010/51‑1 (MQ=255) ccagcCTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg > 2:599993/4‑150 (MQ=255) gCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt > 2:176910/1‑116 (MQ=255) gCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt < 1:176910/116‑1 (MQ=255) cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg > 1:16631/1‑61 (MQ=255) cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg < 2:16631/61‑1 (MQ=255) cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACt > 2:403033/1‑150 (MQ=255) ggCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTg < 1:64927/150‑1 (MQ=255) gCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGc < 1:403033/150‑1 (MQ=255) gCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGc > 1:467373/1‑150 (MQ=255) | CCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGC > NC_002947/1526485‑1526774 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |