Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,002,474 | A→T | 10.9% | intergenic (‑118/+28) | amaA ← / ← amaB | L‑pipecolate oxidase/L‑piperidine‑6‑carboxylate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,002,474 | 0 | A | T | 10.9% | 118.7 / 7.3 | 46 | intergenic (‑118/+28) | amaA/amaB | L‑pipecolate oxidase/L‑piperidine‑6‑carboxylate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (13/28); new base T (2/3); total (15/31) | |||||||||||
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.26e-01 |
GACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCC > NC_002947/6002329‑6002622 | gACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGcc > 2:575813/1‑150 (MQ=255) aGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCgg < 2:36509/150‑1 (MQ=255) tCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCTGGCCCCCggttc < 1:667733/150‑4 (MQ=255) cGGACATGGCCGGCTCCAGCAATGTGTTCTTGTTGGGGCCGCAGGGCGGCCGTTTGGCGGGTAAACCCGCTCCTACAGGGGTCGGTGGCGCTTGGACGGCCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCAt < 2:67276/150‑1 (MQ=255) ggACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATg < 2:742192/150‑1 (MQ=255) gACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCTGGCCATgc < 1:575813/150‑1 (MQ=255) ggCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgc < 2:739438/145‑1 (MQ=255) ggCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgc > 1:739438/1‑145 (MQ=255) gCCGGCTCCAGCAATTTGTTCTTGTGGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATc < 2:539119/150‑1 (MQ=255) ccGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCa > 2:252070/1‑150 (MQ=255) cGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGGCGGTGCCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCaa < 2:822436/150‑1 (MQ=255) ggcgagttGTTCTGGTTGGGGCCGCAGGGAGGCGCTTTCGCGGGTCAACCGGTGTCTACGGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGGGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCa < 2:794995/144‑1 (MQ=255) tttGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGa > 1:633079/1‑150 (MQ=255) cTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGa < 1:252070/150‑1 (MQ=255) cTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCCCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGa < 1:103445/150‑1 (MQ=255) cTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCCCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGa > 2:103445/1‑150 (MQ=255) ggCAGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGTGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTgc < 1:249166/150‑1 (MQ=255) cggcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCa < 2:633079/150‑1 (MQ=255) gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGc < 2:347641/150‑1 (MQ=255) ttCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCGAGCTGCCGCTggcg < 1:632936/150‑1 (MQ=255) cccGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAGCCGGGGGCCTGGCCCCCggtt > 2:368375/1‑81 (MQ=255) cccGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAGCCGGGGGCCTGGCCCCCggtt < 1:368375/83‑3 (MQ=38) ccGCTCCTACAGGGGTCGGTGTCGCTTGGACGAGCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTc > 1:828057/1‑104 (MQ=255) ccGCTCCTACAGGGGTCGGTGTCGCTTGGACGAGCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTc < 2:828057/104‑1 (MQ=255) ccGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTNGGAGCGNGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACg > 1:605978/1‑150 (MQ=255) tcctccAGGGGTCGGTGGCGCTTGTACGAACTCTGTCGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGgtgt < 2:28442/150‑1 (MQ=255) tCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGgtgt < 2:828792/150‑1 (MQ=255) cTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATAACCCCAACACGATTCCCTGCGCCAGCGGCAGCTCGCGCGAGTACTTGAGGGTCtta < 2:32899/150‑2 (MQ=255) tACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTgg > 2:216367/1‑150 (MQ=255) ggggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTg < 2:381784/150‑1 (MQ=255) gggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGagtagt < 1:172627/133‑1 (MQ=255) gggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGagtagt > 2:172627/1‑133 (MQ=255) gggTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTgg < 2:787091/150‑1 (MQ=255) gAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTcc < 2:749374/150‑1 (MQ=255) gcccTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGCCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg < 1:698259/149‑1 (MQ=255) cTGTAGGAGCGGGTCTACCCGAGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTcgc > 1:81117/1‑100 (MQ=255) cTGTAGGAGCGGGTCTACCCGAGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTcgc < 2:81117/100‑1 (MQ=255) gCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGcc > 1:536712/1‑150 (MQ=255) gTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGACGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGAc > 2:636175/1‑150 (MQ=255) aCCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCAc < 2:536712/150‑1 (MQ=255) ccGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACgg < 2:605978/150‑1 (MQ=255) gccGGGGGCCTGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTg > 1:458589/2‑96 (MQ=255) gccGGGGGCCTGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTg < 2:458589/95‑1 (MQ=255) ccAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTc < 1:49516/150‑1 (MQ=255) cAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTcc > 1:851644/1‑150 (MQ=255) cAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTcc < 1:509233/150‑1 (MQ=255) | GACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGACTCACGGCCACCGCCGGTTTCC > NC_002947/6002329‑6002622 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |