Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,790,681 | A→T | 13.1% | intergenic (‑29/+28) | PP_5073 ← / ← hemE | hypothetical protein/uroporphyrinogen decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,790,681 | 0 | A | T | 13.1% | 125.8 / 11.0 | 47 | intergenic (‑29/+28) | PP_5073/hemE | hypothetical protein/uroporphyrinogen decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (25/15); new base T (3/3); total (29/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.66e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01 |
GCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG > NC_002947/5790542‑5790825 | gcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGGGACAACGGGTTGTCATGATAATCGGATACTGTTGCGGTTTCG‑CCAAACATGGTGAAGCCGCTATCATTGAATTCAAACTTAACGTACAAAAAACCCGGCCTCGAGGCCGGGTTTa > 1:292494/1‑150 (MQ=255) ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGaa < 2:418525/139‑1 (MQ=255) ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGaa > 1:418525/1‑139 (MQ=255) tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcag > 2:686010/1‑150 (MQ=255) gCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTc < 1:475834/150‑1 (MQ=255) aTAACGGATTTTCATGATAATCGGGTCCTGTTTTGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGcc < 2:796058/150‑1 (MQ=255) tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTg > 2:1024582/1‑150 (MQ=255) cATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc < 2:693933/135‑1 (MQ=255) cATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc > 1:693933/1‑135 (MQ=255) cATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGACATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGc > 1:316459/1‑150 (MQ=255) aTCGGGTCCTCTTGCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGt < 1:290225/150‑1 (MQ=255) tCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt > 1:649569/1‑150 (MQ=255) ggTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACTTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGGGCACGGCGTtgat < 1:347241/150‑1 (MQ=255) tGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATg < 1:1024582/150‑1 (MQ=255) gtttcgATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg > 1:49278/1‑105 (MQ=255) gtttcgATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg < 2:49278/105‑1 (MQ=255) tttcgATGGGACCAGAT‑TAGCGACGATCACTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCa < 2:292494/150‑1 (MQ=255) aCCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttt < 2:609601/68‑2 (MQ=37) aCCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttt > 1:609601/1‑67 (MQ=37) gAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAAc > 2:888421/1‑150 (MQ=255) gAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAAc < 2:1113478/150‑1 (MQ=255) gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGTGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGc > 1:890721/1‑150 (MQ=255) tCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTgct < 1:686010/150‑3 (MQ=255) ttCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGc > 1:41727/1‑150 (MQ=255) tCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGcc > 2:359695/1‑150 (MQ=255) aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa > 1:521138/1‑150 (MQ=255) aCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg < 2:41727/150‑1 (MQ=255) tACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAacac > 1:198810/1‑150 (MQ=255) cAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt > 1:124845/1‑150 (MQ=255) agtaaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt < 1:929721/81‑1 (MQ=38) agtaaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt > 2:929721/4‑84 (MQ=38) aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGTCGACAGTTCATGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCAGCGTGATGCCGTCGCCAAGGTTGAACACATg > 2:160789/1‑150 (MQ=255) aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa < 1:1206629/137‑1 (MQ=255) aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa > 2:1206629/1‑137 (MQ=255) aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc > 2:260416/1‑150 (MQ=255) aaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGcc > 2:1047637/1‑150 (MQ=255) aaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCa > 2:1161666/1‑150 (MQ=255) aaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCa < 1:1161666/150‑1 (MQ=255) aCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACt < 1:983247/45‑1 (MQ=25) aCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACt > 2:983247/1‑45 (MQ=25) aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAg > 1:108057/1‑150 (MQ=255) aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAg > 2:858773/1‑150 (MQ=255) aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAg > 1:192102/1‑150 (MQ=255) ccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAAc > 2:1260818/1‑107 (MQ=255) ccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAAc < 1:1260818/107‑1 (MQ=255) cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc > 1:64290/1‑150 (MQ=255) aTCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTg > 1:579873/1‑150 (MQ=255) aTCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTg > 1:157673/1‑150 (MQ=255) aTCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATAAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTg < 1:850853/150‑1 (MQ=255) | GCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGAT‑TAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTG > NC_002947/5790542‑5790825 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |