Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,528,397 | T→A | 8.4% | 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,397 | 0 | T | A | 8.4% | 173.0 / 5.7 | 59 | intergenic (+28/‑48) | purB/PP_4017 | adenylosuccinate lyase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (28/25); new base A (2/3); total (31/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01 |
GAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGT > NC_002947/4528250‑4528539 | gAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTgg > 2:1335467/1‑150 (MQ=255) cgcTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCgg < 2:370616/150‑1 (MQ=255) cgcTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCgg < 2:762536/150‑1 (MQ=255) tgACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGttg < 1:621486/150‑2 (MQ=255) gACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGGCGAGCTCAAGCAGCTGACCCCCGCTACCTACTTCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCCGGCCCGGGGgttt > 2:628726/1‑150 (MQ=255) aCCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGttttt > 1:1038635/1‑150 (MQ=255) tCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCg > 1:1005810/1‑150 (MQ=255) aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc < 1:972427/150‑1 (MQ=255) aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAACGCCCGGCCAGGCCGGGCGtt > 1:720646/1‑137 (MQ=255) aCGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAACGCCCGGCCAGGCCGGGCGtt < 2:720646/137‑1 (MQ=255) gCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaa > 2:1449372/1‑150 (MQ=255) gCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaa < 2:1148468/150‑1 (MQ=255) ccTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaaa < 1:1338056/150‑1 (MQ=255) cATGCCAGCCGAAGCCAAGGCCGAGCCCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACg < 1:532150/150‑1 (MQ=255) aTGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGt > 2:810705/1‑150 (MQ=255) tGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGGCTCTCGCTACCTACATCGGTAACGCGGCAGCCCAAGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtt < 1:84780/150‑1 (MQ=255) tGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtt < 1:376006/150‑1 (MQ=255) gAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAAt < 2:1038635/150‑1 (MQ=255) ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc > 1:1583782/1‑120 (MQ=255) ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc < 2:1583782/120‑1 (MQ=255) ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc < 2:1583800/120‑1 (MQ=255) ccAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTc > 1:1583800/1‑120 (MQ=255) aaGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTg > 1:890417/1‑150 (MQ=255) gCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGtt < 1:485728/103‑1 (MQ=255) gCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCAGGCCGGGCGtt > 2:485728/1‑103 (MQ=255) cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAAt < 1:1193699/150‑1 (MQ=255) tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATc > 2:996467/1‑150 (MQ=255) gcagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa < 1:996467/150‑1 (MQ=255) gcagcTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGAGTATCATACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGa < 1:1335467/150‑1 (MQ=255) tGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTc < 1:810705/150‑1 (MQ=255) cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAgctgc > 1:512400/1‑150 (MQ=255) acATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc > 1:1281087/1‑127 (MQ=255) acATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc < 2:1281087/127‑1 (MQ=255) ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc > 1:1193889/1‑85 (MQ=255) ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTAc < 2:1193889/85‑1 (MQ=255) gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt < 2:1148052/100‑1 (MQ=255) gTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGtttttt > 1:1148052/1‑100 (MQ=255) gcgGCAGGCCAGGCCAAACGCATCTAAGCTCAGCGGCGCGTACAACGCCCGGCCTAGCCGGGGGTTTTTATACCCGTACGAAAACTACGGTGTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCg > 2:440490/1‑150 (MQ=255) gcgGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACTTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTTCAGCTGCCAGGCGGCAtcttt < 2:730880/150‑3 (MQ=255) cccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCcgcg > 1:1414104/1‑150 (MQ=255) ccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGa > 1:91911/1‑150 (MQ=255) aaaCGCATCTAAGCTGAGCGTCGCGTANAACGCCCGGCCTNGCCGGGNGTTTTTATTCCCGTACGAAAATTACTTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCTGCGGCATCTCGGCCCGCGAATTCATg > 1:1299078/1‑150 (MQ=255) aCGCATCTAAGCTGATCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATgcg > 2:843393/1‑150 (MQ=255) aTCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGac > 1:683140/1‑150 (MQ=255) tAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc > 1:506241/1‑117 (MQ=255) tAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTcggc < 2:506241/117‑1 (MQ=255) gcgcagcgTACAACGCGCGCCCTGGCCGGGCGTTTTTATTCCCGTACCAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACCCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaa < 1:628726/145‑1 (MQ=255) tCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGTTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAgaagaa > 2:1350316/1‑150 (MQ=255) tCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAgaagaa > 1:685146/1‑150 (MQ=255) cgcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGaagaag < 2:512400/150‑1 (MQ=255) gcgTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGc < 2:262828/150‑1 (MQ=255) gTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaa < 1:160645/143‑1 (MQ=255) gTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCagaa > 2:160645/1‑143 (MQ=255) tACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgc > 1:1034012/1‑150 (MQ=255) cAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCGCGCGAATTCATGCGCGACTAATGGAAGAAGAAGCCgctg > 1:407454/1‑150 (MQ=255) aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGACAATTACGTTTTTTCAATAGGTTCACCATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAAATCATGCGCGACTACTAGCGGAAGAAGCCgctgc > 2:1333489/1‑150 (MQ=255) aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAAt < 1:726514/117‑1 (MQ=255) aaCGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAAt > 2:726514/1‑117 (MQ=255) cccGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTTTTTTTTTTTAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGt > 2:156348/1‑150 (MQ=255) | GAAGCGCTGCTGACCTTCATCGACGGCCTGGACATGCCAGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGT > NC_002947/4528250‑4528539 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |