Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,082,692 | A→C | 11.1% | L133R (CTG→CGG) | purK ← | 5‑(carboxyamino)imidazole ribonucleotide synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,082,692 | 0 | A | C | 11.1% | 89.3 / 4.2 | 36 | L133R (CTG→CGG) | purK | 5‑(carboxyamino)imidazole ribonucleotide synthase |
Reads supporting (aligned to +/- strand): ref base A (20/12); new base C (2/2); total (22/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.34e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 |
GAACGGCACGAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGGCGCTGGGGTACACG > NC_002947/6082556‑6082838 | gAACGGCACGAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCGTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGc < 1:12400/150‑1 (MQ=255) cGGCACGAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCATCGGCCTGCCGGCCG‑TGCTgaaaac < 2:359451/150‑6 (MQ=255) aCGAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTgcg > 2:743322/1‑150 (MQ=255) gAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTgcggc > 1:254264/1‑150 (MQ=255) ccTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCa < 2:1045811/150‑1 (MQ=255) cTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCCCTCTGATGCTGGCCACTGCGGCGTCCAg > 1:840845/1‑150 (MQ=255) cTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTgcggcg > 2:389337/1‑145 (MQ=255) cTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTgcggcg < 1:389337/145‑1 (MQ=255) cTCCAGCAGGCATGGCACGCTGCCCAGCTAGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAg > 1:987669/1‑150 (MQ=255) tCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTg > 2:737014/1‑132 (MQ=255) tCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTg < 1:737014/132‑1 (MQ=255) cagcagGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGGCGAGGGTGGCCACTGCGGCGTCCAGATc > 2:140273/1‑150 (MQ=255) ggCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCt < 1:1042643/150‑1 (MQ=255) cAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGc > 2:665807/1‑150 (MQ=255) cGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCATCGGCCTGCCGGCCG‑TGCTgaaaac > 1:123431/1‑101 (MQ=37) cGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCATCGGCCTGCCGGCCG‑TGCTgaaaac < 2:123431/106‑6 (MQ=37) tACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTgg < 1:743322/150‑1 (MQ=255) aCCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTggg < 1:665807/150‑1 (MQ=255) cccACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTgggg > 2:467450/1‑150 (MQ=255) ccaccaCGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCAGTGCGGCGTCCAGATCGGCCTTCGAGAGGATGTAGGCGAAGGCCGGGGTggggg > 1:402836/1‑149 (MQ=255) caccacGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGAt > 1:390211/1‑150 (MQ=255) gCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGa > 1:838691/1‑150 (MQ=255) gcagcaCCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGc > 2:1206/1‑150 (MQ=255) ccTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCtttt > 1:187726/1‑150 (MQ=255) ttCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCg > 1:36200/1‑150 (MQ=255) cccTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAgg < 1:275767/150‑1 (MQ=255) tAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCg > 1:263345/1‑150 (MQ=255) aCCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGt > 2:324822/1‑140 (MQ=255) aCCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGt < 1:324822/140‑1 (MQ=255) cAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATgcg > 2:129466/1‑150 (MQ=255) ttttCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCgg > 2:27308/1‑150 (MQ=255) ttCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGgcg < 2:648152/150‑1 (MQ=255) gCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGGCGCTgg < 2:263345/150‑1 (MQ=255) a‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGGCGCTgggg < 2:74514/150‑1 (MQ=255) ccggcAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGGCGCTGGGGTaca > 1:863644/1‑150 (MQ=255) gtcAGGCCGATGCTGGCAACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGCGGTGGGGATGCGCAGGTCGCCGAACATGCTCTGTTCGAACGGGCGGTCAAGAGCGATGCGCTCCAATTTGGCGCTGGGGTACACg > 2:193203/3‑150 (MQ=255) | GAACGGCACGAAGCCCTCCAGCAGGCATGGCACGCTGCCCAGCTCGGCGAAGGTACCCACCACGTCTTCTGGCGTGCGCAGCACCTTCTGGCCCTTGCCGTCGTAACCCAGGGTGCGGGTTTTCAGCA‑CGGCCGGCAGGCCGATGCTGGCCACTGCGGCGTCCAGATCGGCCTGCGAGAGGATGTCGGCGAAGGCCGGGGTGGGGATGCCCAGGTCGCGGAACATGCTCTTTTCGAACAGGCGGTCACGAGCGATGCGCAGCGATTCGGCGCTGGGGTACACG > NC_002947/6082556‑6082838 |
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: ‑ |