Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,621,712 | T→C | 10.5% | intergenic (+305/+4259) | PP_5654 → / ← PP_4091 | phosphotriesterase/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,621,712 | 0 | T | C | 10.5% | 99.3 / 3.2 | 38 | intergenic (+305/+4259) | PP_5654/PP_4091 | phosphotriesterase/transposase |
Reads supporting (aligned to +/- strand): ref base T (14/20); new base C (2/2); total (16/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01 |
CTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATA > NC_002947/4621564‑4621859 | cTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGTTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAAt < 6:132650/149‑1 (MQ=255) gTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAAtt > 3:299198/1‑149 (MQ=255) aaTTACGAGCGCACCACACAAACGCACCTAGCTAACCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa < 8:5151/149‑1 (MQ=255) cGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACAt < 7:55343/149‑1 (MQ=255) gAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTTCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACTACTGGCCTCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATc < 3:237598/149‑1 (MQ=255) cgcACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAg < 4:299198/149‑1 (MQ=255) accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTc < 3:31030/126‑1 (MQ=255) accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTc > 4:31030/1‑126 (MQ=255) accacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgtg > 2:94458/1‑148 (MQ=255) ccacACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTa < 4:53250/149‑1 (MQ=255) acaAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgt < 2:184719/142‑1 (MQ=255) acaAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAgt > 1:184719/1‑142 (MQ=255) cAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAg < 3:115601/137‑1 (MQ=255) aaGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGc > 2:247997/1‑129 (MQ=255) aaGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGc < 1:247997/129‑1 (MQ=255) aaaaaaaCGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGc > 1:64311/1‑147 (MQ=255) cAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATTTCGCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGgc < 1:94458/149‑1 (MQ=255) aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc > 2:275556/1‑122 (MQ=255) aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc > 2:280063/1‑122 (MQ=255) aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc < 1:280063/122‑1 (MQ=255) aaCAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAACGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAgcgc < 1:275556/122‑1 (MQ=255) aTAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGc > 3:260020/1‑149 (MQ=255) tGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa < 4:263372/73‑1 (MQ=255) tGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTa > 3:263372/1‑73 (MQ=255) aCTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCgg < 1:125240/148‑1 (MQ=255) ccTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGAt > 6:238612/1‑149 (MQ=255) cTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATc > 5:174368/1‑149 (MQ=255) cTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATc < 6:174368/149‑1 (MQ=255) tCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATgg < 2:109115/115‑1 (MQ=255) acaacTGGCCGCGAACTGTGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCcgc < 4:309998/149‑1 (MQ=255) aacTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATAGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGa > 1:139037/1‑149 (MQ=255) ccGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCt < 2:64311/149‑1 (MQ=255) cgcgAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTc < 4:260020/149‑1 (MQ=255) cgcgAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTc < 1:9885/149‑1 (MQ=255) ggAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGc > 8:93482/1‑149 (MQ=255) gCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGaa < 5:238612/149‑1 (MQ=255) atatGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAacac < 2:139037/148‑1 (MQ=255) atatGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTCATCGAACCCCGCGATCCCCTCGGCCGCAGCGAaccc > 2:83218/1‑146 (MQ=255) aTGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATa > 4:224415/1‑149 (MQ=255) | CTAGCGCAATGTTAAATTACGAGCGCACCACACAAACGCACCTAGCATCCCAAGAAACATAGAAAAAAACGGTGCGCAACAGATAGCCCCTGAGAAAACTTGTTGCCTCCAACTAACAACTGGCCGCGAACTGGGAGCTTATATGCAATGTCCTTAATT‑‑GCTACACATCTAGTGTATCGCAGGACAGCTGATCAAGCGCATCGATCCGCAAGGGGGGAATGGCGCTACAGCTACGACGATATCGGCCGCCTGATCGAAACCCGCGATCCCCTCGGCCGCAGCGAACACATCAAATA > NC_002947/4621564‑4621859 |
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: ‑ |