Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,622,914 | A→G | 9.1% | intergenic (+1507/+3057) | 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,622,914 | 0 | A | G | 9.1% | 112.1 / 2.3 | 46 | intergenic (+1507/+3057) | PP_5654/PP_4091 | phosphotriesterase/transposase |
Reads supporting (aligned to +/- strand): ref base A (24/16); new base G (2/2); total (26/20) | |||||||||||
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.82e-01 |
AGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCACTGCT > NC_002947/4622767‑4623062 | aGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGAt > 3:136596/1‑149 (MQ=255) cTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGt < 5:262234/149‑1 (MQ=255) gCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGAt < 2:299659/149‑1 (MQ=255) gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc < 6:175273/134‑1 (MQ=255) gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc < 6:179587/134‑1 (MQ=255) gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc > 5:179587/1‑134 (MQ=255) gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc > 5:175273/1‑134 (MQ=255) cGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCAcgacga < 3:74282/149‑1 (MQ=255) cccGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCATATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTgc < 1:51012/149‑1 (MQ=255) gCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCACGCATTGCAAACGTGATCCGGTGCACGCCCCGGTCGTGCGCCCCCag > 1:131864/1‑148 (MQ=14) gCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCAACTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGCTGCGGTGCACGACGAGGTCCTGCGCAacc > 8:71653/1‑147 (MQ=255) aaTCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAgg > 4:92825/1‑149 (MQ=255) gCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAAc > 5:130253/1‑149 (MQ=255) cAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACt > 7:7370/1‑149 (MQ=255) cAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACt > 7:2717/1‑149 (MQ=255) aGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTg < 4:136596/149‑1 (MQ=255) cttgTCCTGGCCGAACAAAGTGATCGCAACCCACTGTTC‑TTCGGCAGGGGGCACCTGCAACATTTCAAACTTAAAGGGCGGGGTTTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCAc < 2:131864/146‑1 (MQ=2) cTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCAc > 6:126720/1‑149 (MQ=39) cTGCCCGACCAACGGGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc < 8:209762/149‑1 (MQ=21) cccGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAcgcg > 4:269568/1‑149 (MQ=37) cGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACgcgc > 3:169396/1‑148 (MQ=37) gACCAACGTGAACTCAACCACCTGTACTTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc > 1:20376/1‑144 (MQ=255) gACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc < 2:20376/143‑1 (MQ=37) aaCGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGa > 3:98710/1‑148 (MQ=37) gTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGccac > 8:284852/1‑148 (MQ=25) tCAACCACCTGTAC‑TACGGGAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCATGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCgga > 7:4402/1‑148 (MQ=21) tCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGt > 7:184837/1‑149 (MQ=34) aaCCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCg > 4:90284/1‑149 (MQ=34) tGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGcc > 4:271767/1‑149 (MQ=37) gTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCg > 3:215433/1‑149 (MQ=37) tacGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAGGCGCTACGACACCAGCTGTCGCCTGGCCGGaaa < 3:185740/149‑1 (MQ=255) cGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGGTTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGaaa < 5:119730/147‑1 (MQ=255) cGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGaaa < 7:46643/147‑1 (MQ=39) cGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt > 5:328270/1‑148 (MQ=39) cTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACcg < 5:126720/149‑1 (MQ=39) cTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACcg > 4:145586/1‑149 (MQ=39) tGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTAACTAcc > 3:230864/1‑147 (MQ=255) ccccAGAGCAACCTCAACGGCCGGGTGTTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGa < 4:169396/147‑1 (MQ=255) cAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGc < 6:328270/149‑1 (MQ=39) gAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAg > 1:148359/1‑149 (MQ=39) ccGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGtgctgc < 8:184837/149‑1 (MQ=39) ccGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGATCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGtgctgc < 6:130253/149‑1 (MQ=255) gggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt > 7:137081/1‑116 (MQ=255) gggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt < 8:137081/116‑1 (MQ=255) ggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCAc > 8:183210/1‑149 (MQ=40) ggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCAc < 3:145586/149‑1 (MQ=40) aTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTTCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCActgct > 1:309472/1‑149 (MQ=255) | AGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTAC‑TACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCACTGCT > NC_002947/4622767‑4623062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |