Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,312,689 | T→A | 17.4% | intergenic (+30/‑132) | PP_2912 → / → hemB | hypothetical protein/5‑aminolevulinate dehydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,312,689 | 0 | T | A | 17.4% | 111.8 / 15.0 | 56 | intergenic (+30/‑132) | PP_2912/hemB | hypothetical protein/5‑aminolevulinate dehydratase |
Reads supporting (aligned to +/- strand): ref base T (22/24); new base A (3/7); total (25/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01 |
GTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT > NC_002947/3312550‑3312834 | gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa < 1:133629/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa < 8:247839/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa < 3:364859/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCcgcg > 2:303536/1‑147 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa < 4:28116/149‑1 (MQ=255) gagaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAaga < 7:282010/149‑1 (MQ=255) agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAagag < 1:199908/149‑1 (MQ=255) gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg > 2:58871/1‑149 (MQ=255) cGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa > 8:248989/1‑142 (MQ=255) cGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa < 7:248989/142‑1 (MQ=255) aCGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcg < 4:158971/146‑2 (MQ=255) aCGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTACCCCGCGAAGAGGcggc > 3:158971/1‑145 (MQ=255) cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg > 2:320221/1‑149 (MQ=255) ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca < 3:307722/145‑1 (MQ=255) ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCAGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca > 4:307722/1‑145 (MQ=255) gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa > 5:313517/1‑126 (MQ=255) gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa < 6:313517/126‑1 (MQ=255) gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg < 4:319359/149‑1 (MQ=255) gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg > 3:22030/1‑149 (MQ=255) gacgCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑cc > 5:77153/1‑149 (MQ=255) tcgtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGc > 5:150513/1‑149 (MQ=255) ctcggcgcGTGCTGCACGGCGACGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTTCCCGCGCAGAGGCCAACACCGGCAAACCCCCG‑‑‑CCGTTg > 2:167603/8‑149 (MQ=255) gtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCa < 1:320221/149‑1 (MQ=255) cgccgcGCAGCAC‑GCGCCGAGGCCGAGCAAGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCAGGCTGCAGTAGTTTATTCCTGTACCCCCTCTTCGCGGTGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCAtt < 2:121483/149‑1 (MQ=255) cGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGAc < 1:58871/149‑1 (MQ=255) gcaC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGa > 3:224476/1‑149 (MQ=255) gcaC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGa < 5:225418/149‑1 (MQ=255) aC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCATCGCGGGGTTTCCCGCGAAGAGGCCAACACAGGCAAACCCCGG‑‑‑CCGTTGCCATTGACGCAAcaca > 3:355766/1‑149 (MQ=255) c‑gcgCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCGAAGAGGGCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGaca > 3:47698/1‑148 (MQ=255) cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTG‑CATTGACGCAAGACACATtctc < 4:129861/149‑1 (MQ=255) gaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATtctc < 7:297287/149‑1 (MQ=255) gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCAtt < 4:224476/149‑1 (MQ=255) ccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt < 4:22030/149‑1 (MQ=255) gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTACCCGCGAAGCGGCCAACACACGCCAACCCCGG‑‑‑CCGTTGGCATTGACGCAAGACACATTCTAGTCATTTAGCCttttgcc > 6:215773/1‑146 (MQ=255) gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCCAAGAGGCCCACACAGGCAAACCACGG‑‑‑CCGTTGCCAATGCCGCAAGACACCTTGTCGTCATTTAGCCTTTTCACTGgggc > 4:249580/1‑149 (MQ=255) gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTTACCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCCTTGCCATTGACGCAAGACCCATTATAGGGATTTATCCTTTTAACTGgcgg > 2:152995/1‑148 (MQ=255) gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGa > 8:53253/1‑147 (MQ=255) gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa < 1:167603/149‑1 (MQ=255) gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg < 6:150513/149‑1 (MQ=255) gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg < 4:198717/149‑1 (MQ=255) gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg < 3:249580/149‑1 (MQ=255) gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACATGCAAACCACGGGACCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg < 1:303536/149‑1 (MQ=255) tATTCCTGTACCGCCGATTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTcc < 2:86224/148‑1 (MQ=255) ttCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCt < 4:47698/148‑1 (MQ=255) ggatctccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa < 5:129266/143‑1 (MQ=255) cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAg > 2:261972/1‑149 (MQ=255) cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGGTGCCATTGACGCAAGACACATTCTCCTCATTTAGCCTTTTCACTGGCGCCCGACAGGCATAAGGTATAAGCCGTTCAAATTCCCTg > 1:66974/1‑146 (MQ=255) tGTACCGCCTCTTCCCGGGGTTACCCGCGAAGAGCCCAACACAGGCAAACCACGG‑‑‑CCGTTTCCAATGACCCAAGACACAATCTCGTCATTTAGACTTTTCACTGGCGAGCCACAAGCATAAGGTATAAGCCGATCAAATTCCATgtac > 8:327621/1‑145 (MQ=255) gTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt > 5:313066/1‑149 (MQ=255) ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc < 1:261972/149‑1 (MQ=255) cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCAGTTCAAATTCCCTGGAGCTGTCa > 4:238526/1‑149 (MQ=255) gCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt < 2:123486/146‑1 (MQ=255) gCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGAGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt > 1:123486/1‑146 (MQ=255) ctctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt < 1:152995/149‑1 (MQ=255) gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac < 5:215773/148‑1 (MQ=255) ggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGGCAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAAtt > 1:330448/1‑149 (MQ=255) | GTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT > NC_002947/3312550‑3312834 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |