Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,414,785 | A→G | 14.8% | intergenic (+327/‑18) | PP_2114 → / → PP_2115 | transposase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,414,785 | 0 | A | G | 14.8% | 51.7 / 3.7 | 28 | intergenic (+327/‑18) | PP_2114/PP_2115 | transposase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (8/15); new base G (2/2); total (11/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.13e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01 |
ACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAA > NC_002947/2414674‑2414933 | aCCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTc < 6:275080/148‑1 (MQ=31) ccAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCttgcg < 3:200060/149‑4 (MQ=11) cAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCAGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCCTCACCTGCCATCGc > 8:243638/1‑148 (MQ=9) aaaaaCGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTTCGGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCaa < 4:39071/149‑1 (MQ=12) aCGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCt < 6:170241/149‑1 (MQ=32) cGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTg < 3:162239/149‑1 (MQ=32) cgcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGtt < 8:1820/148‑1 (MQ=32) gcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCg < 1:90746/149‑1 (MQ=32) cgccGGCCCCAGCGCGGCGCTCCCGGGGCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGa < 3:245248/147‑1 (MQ=255) cGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCAAACGGGGGCC‑TCCGCGCTTGACGGGCATCCATGCCCCTCACCTGCCTTCGCAAGGCCTGCGTTAGGCCGCCGGaaat > 6:64750/1‑147 (MQ=11) gcgctgcgctCCGGGGTCCCCTCACTCCGGCCATGGTCCCGGCCGGACCGCGCCCACAGCCCCAACCTGGGGCC‑GAAGCGCATGAAGGGCATCCATGACCCTAACCAGACATACCAAAGCCTGCGTTGGGCCGCCAGAAGTAGCGaag > 4:297914/1‑148 (MQ=2) tgcgctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCCTCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTGCTGAAGTCGCGaagatc > 1:22963/1‑147 (MQ=17) cgctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaaga < 6:238881/149‑1 (MQ=32) cgctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa < 4:312387/147‑1 (MQ=32) cgctCC‑GGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaaga < 5:32712/148‑1 (MQ=17) ctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCGGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGc > 3:258469/1‑112 (MQ=17) ctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCGGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGc < 4:258469/112‑1 (MQ=17) gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatc < 5:64750/147‑1 (MQ=32) cAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaaga < 1:307465/149‑1 (MQ=32) gcCGAAGGCCCCATCCTTGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCAGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatca < 7:243638/149‑1 (MQ=17) ctgaGGCCCCAGGATGGGGCCTTCGGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaa > 1:214660/4‑149 (MQ=11) ccccATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatc > 1:136974/1‑146 (MQ=35) cccATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGcagatcaagatcaagatcaagatcaatctcaagatcaagatcacgatgaggatcaagatcaagc > 8:112351/1‑148 (MQ=11) ggCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTTCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAATTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaa < 8:226257/149‑1 (MQ=255) gCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatc > 3:146893/1‑146 (MQ=255) gCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatc > 3:141999/1‑146 (MQ=255) cc‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaaga < 2:22963/149‑1 (MQ=255) cAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatc > 8:42215/1‑148 (MQ=255) aGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAATGCCTGCGTTCGGCCTACTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaa < 3:276205/149‑1 (MQ=255) | ACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCC‑TCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAA > NC_002947/2414674‑2414933 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |