Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,474,268 | G→T | 11.7% | T423K (ACG→AAG) ‡ | PP_3964 ← | transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,474,268 | 0 | G | T | 11.7% | 138.7 / 7.3 | 49 | T423K (ACG→AAG) ‡ | PP_3964 | transposase |
Reads supporting (aligned to +/- strand): ref base G (18/25); new base T (3/3); total (21/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01 |
GATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC > NC_002947/4474123‑4474416 | gATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt > 3:69263/1‑147 (MQ=255) aaaCTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAgcgc < 4:250464/149‑1 (MQ=255) tCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCTGACGCGCTAcc > 3:71631/1‑149 (MQ=255) tCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTAcc > 8:34816/1‑149 (MQ=255) aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCTGACGCGCTAcctgga < 4:71631/148‑5 (MQ=18) aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCg > 2:163818/1‑149 (MQ=255) aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCg < 5:33590/149‑1 (MQ=255) cGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCa < 5:212193/149‑1 (MQ=255) gCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGTCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAg < 7:192412/149‑1 (MQ=255) ccGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga < 8:116761/137‑5 (MQ=18) ccGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga > 7:116761/1‑133 (MQ=18) gCTTACCGCTGCGCGGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTaa < 3:146850/149‑1 (MQ=255) aCCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGAGCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATc < 2:171489/147‑1 (MQ=255) ccGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGa < 1:163818/148‑1 (MQ=255) ccGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAg < 6:165154/149‑1 (MQ=255) tGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTAGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCt > 6:268906/1‑149 (MQ=255) gATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGt < 1:34876/149‑1 (MQ=255) ccGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga < 7:228526/114‑5 (MQ=14) ccGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga > 8:228526/1‑110 (MQ=14) gCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGa > 5:196930/1‑149 (MQ=255) gCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGa < 8:158664/149‑1 (MQ=255) acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg < 4:14783‑M2/85‑1 (MQ=255) acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg > 3:14783‑M2/59‑143 (MQ=255) aaaTAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTCGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATc < 5:268906/149‑1 (MQ=255) aTAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCa > 2:100014/1‑149 (MQ=255) caaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCAt > 4:17623‑M2/55‑149 (MQ=255) ACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCAtt > 4:69263‑M1/1‑148 (MQ=255) CCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt > 4:236404‑M1/1‑90 (MQ=255) CCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt < 3:236404‑M1/90‑1 (MQ=255) tGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAg < 7:34816/149‑1 (MQ=255) ccATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGcc > 7:71768/1‑148 (MQ=255) aTGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAAc < 1:100014/149‑1 (MQ=255) cGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATc > 3:20760/1‑147 (MQ=255) gTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATg < 4:20760/147‑1 (MQ=255) aCCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt < 8:71768/149‑1 (MQ=255) ctggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa < 6:196930‑M2/143‑1 (MQ=255) ctggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa < 1:249590‑M2/143‑1 (MQ=255) tGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATc > 5:6270‑M2/2‑149 (MQ=255) ttgtCGATGGGCACAGGCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCg < 7:229544/148‑1 (MQ=255) tCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt < 5:73290/137‑1 (MQ=255) tCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt > 6:73290/1‑137 (MQ=255) gATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCAcc < 8:82684/148‑1 (MQ=255) gATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCg > 3:209745/1‑149 (MQ=255) tGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGc < 4:209745/148‑1 (MQ=255) acaGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc > 7:161987/1‑149 (MQ=255) acaGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc > 3:49560/1‑149 (MQ=255) cccATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTc < 4:49560/149‑1 (MQ=255) cAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTc > 7:225345/1‑102 (MQ=255) cAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTc < 8:225345/102‑1 (MQ=255) cgcgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa < 4:63578/109‑1 (MQ=255) cgcgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa > 3:63578/1‑109 (MQ=255) gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCtt > 5:158004/1‑148 (MQ=255) gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTtc < 3:17623/149‑1 (MQ=255) | GATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC > NC_002947/4474123‑4474416 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |