Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,115,762 | C→A | 13.4% | R24R (CGG→CGT) | PP_4500 ← | TrpH family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,115,762 | 0 | C | A | 13.4% | 69.9 / 4.0 | 30 | R24R (CGG→CGT) | PP_4500 | TrpH family protein |
Reads supporting (aligned to +/- strand): ref base C (11/15); new base A (2/2); total (13/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 = 2.28e-01 |
ATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGC > NC_002947/5115616‑5115896 | atggTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCgg < 2:718174/150‑1 (MQ=255) cgcgCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc < 1:168510/150‑1 (MQ=255) cgcgCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc < 2:154459/150‑1 (MQ=255) cgccgcCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGATGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCa > 1:774006/1‑150 (MQ=255) aGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggggt < 1:198636/131‑3 (MQ=255) aGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggggt > 2:198636/1‑129 (MQ=255) cTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGcc > 2:394259/1‑150 (MQ=255) ccccGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCgg > 2:773549/1‑150 (MQ=255) gCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGc < 1:129858/150‑1 (MQ=255) cgcATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGa < 1:394259/150‑1 (MQ=255) aTACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGAGCGTGGACCCCGTGCTCATGAGCCCGGGCAACCATGACCGAAAGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAg > 1:52834/1‑150 (MQ=255) aCCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGATGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTg < 2:774006/150‑1 (MQ=255) ccAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGa > 1:719014/1‑130 (MQ=255) ccAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGa < 2:719014/130‑1 (MQ=255) ggTCAACGCAGGCCTGGCGCGCTTCGCGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGTGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGAt < 2:314024/150‑1 (MQ=255) gTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCGGCCAGCCCCGGTCAACCAGTACCGAGGGCGCCAGGGCGCCATCGTAGTCGGTGTTGTGACAGTGCAGATc > 1:14186/1‑150 (MQ=255) cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAgagga < 2:30878/133‑1 (MQ=255) cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAgagga > 1:30878/1‑133 (MQ=255) cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTc > 2:559093/1‑150 (MQ=255) ttCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCg > 2:537796/1‑150 (MQ=255) aGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGAtgt < 1:445533/150‑1 (MQ=255) cggtcAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGcc < 2:701823/150‑1 (MQ=255) cTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTc < 1:773549/150‑1 (MQ=255) gCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTACAGACCAACATTCCTAAAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTgg < 2:690027/150‑1 (MQ=255) aCCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCt < 2:52834/150‑1 (MQ=255) cccGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGc > 1:289759/1‑150 (MQ=255) ccgggCTCATGAGCACGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGAt < 2:129939/103‑1 (MQ=255) ccGGGCTCATGAGCACGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGAt > 1:129939/1‑107 (MQ=255) aTGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCt > 2:464045/1‑54 (MQ=255) aTGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCt < 1:464045/54‑1 (MQ=255) | ATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGGCTGC > NC_002947/5115616‑5115896 |
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: ‑ |