Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,528,353 | T→A | 9.6% | intergenic (+28/‑48) | purB → / → PP_4017 | adenylosuccinate lyase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,528,353 | 0 | T | A | 9.6% | 131.7 / 5.8 | 52 | intergenic (+28/‑48) | purB/PP_4017 | adenylosuccinate lyase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (21/26); new base A (2/3); total (23/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01 |
CTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCT > NC_002947/4528215‑4528492 | ctgACCTTCATCGACGGCCTGGACATGCCAGCCGGAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCg < 8:25693/148‑1 (MQ=255) tgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑tt < 6:224232/149‑1 (MQ=255) ggACCTTCATATAAGGCCTGGACATGCCATCCGAAGCCAAGGCCGAGCTCAAGCAGCTTACACCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAAAGCATCTAAGCTGATCGTCTCGTACAACGCCCGGCCTTTCCGGGAG‑tt < 4:10012/148‑1 (MQ=255) cTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTa > 5:190164/1‑149 (MQ=255) aTCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTc < 6:81224/148‑1 (MQ=255) tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTg > 7:162833/1‑130 (MQ=255) tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTg < 8:162833/130‑1 (MQ=255) tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTccc > 1:248672/1‑149 (MQ=255) cTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGaaa < 1:76879/149‑1 (MQ=255) tGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGaaaa > 1:233910/1‑149 (MQ=255) aTGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACg < 7:83489/149‑1 (MQ=255) ccGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTAAGtttttt > 1:304068/1‑148 (MQ=255) gAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGATGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCaa > 3:127565/1‑149 (MQ=255) aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCaa > 2:271501/1‑148 (MQ=255) cAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGt > 6:48540/1‑149 (MQ=255) aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGtt < 5:11619/149‑1 (MQ=255) aGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCG‑tt > 5:357917/1‑110 (MQ=255) aGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCG‑tt < 6:357917/110‑1 (MQ=255) cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTTGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACa < 1:78862/149‑1 (MQ=255) cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGaa > 5:160690/1‑149 (MQ=255) cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGaa < 3:113626/149‑1 (MQ=255) aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt < 7:283736/143‑1 (MQ=255) aaGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACAt > 8:283736/1‑143 (MQ=255) agaagcTTACCCCCGCTACCTACATGGGTAAGGCGGCTGCCCTGGCAAAAGGAATCTATGCTTAGCGTGGCGTAAAACGCCCGGCTTGGCCGGGCGCCTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc < 1:271501/146‑1 (MQ=255) cagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa < 2:233910/149‑1 (MQ=255) cTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAGCGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATAc > 6:159845/1‑149 (MQ=255) cccccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCa < 6:190164/149‑1 (MQ=255) ccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTg > 2:277761/1‑149 (MQ=255) ccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTg > 1:252603/1‑149 (MQ=255) ctacctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAg < 6:244568/149‑1 (MQ=255) cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgc < 5:48540/147‑1 (MQ=255) tacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTTCAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgctgct > 2:75234/1‑149 (MQ=255) tCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg < 8:294734/110‑1 (MQ=255) ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcg < 2:205701/148‑1 (MQ=255) ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggcg < 2:252603/148‑1 (MQ=255) cgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa < 6:54568/118‑1 (MQ=255) cgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa > 5:54568/1‑118 (MQ=255) aGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTAAATAGGTTGAACATGAATCCTGAAAATCCACTGCAGCTGCTAGGCGGCATCTCGGccc > 4:22621/1‑149 (MQ=255) cAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgcg > 1:147673/1‑148 (MQ=255) cAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCa < 3:22621/149‑1 (MQ=255) aaCGCATCTAAGATGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATg < 4:127565/149‑1 (MQ=255) gCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgcgc < 5:159845/149‑1 (MQ=255) gCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgcgc > 4:130884/1‑149 (MQ=255) cATCTTAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGcgcg < 5:92310/149‑1 (MQ=255) aTCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGa < 6:160690/149‑1 (MQ=255) tGAGCGTCGCGTACAAAGCCAGGCATGGCCGGGCG‑TATTTATTCCCGTAAGAAAATTACCTTTTTTCAATAGGTTGAAAATGAATCCTGATAATCCACTGCAGCTGCTAGGCGGCATCTCGGCCAGAGAATTCATGCGCGACTACTgg > 8:95436/1‑148 (MQ=255) gTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCag > 6:332048/1‑146 (MQ=255) gaataaaaaCGCCCGGCCAGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaag > 3:289165/8‑147 (MQ=255) gaataaaaaCGCCCGGCCAGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAgaagaa < 4:289165/142‑1 (MQ=255) ttcAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCTTCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCg < 1:75234/147‑1 (MQ=255) aCAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgc > 4:235877/1‑149 (MQ=255) cGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCTCTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctgct < 7:211262/149‑1 (MQ=255) | CTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCG‑TTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCT > NC_002947/4528215‑4528492 |
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: ‑ |