Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,559,994 | T→A | 11.8% | Q185L (CAG→CTG) | hisZ ← | ATP phosphoribosyltransferase regulatory subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,559,994 | 0 | T | A | 11.8% | 89.0 / 4.7 | 34 | Q185L (CAG→CTG) | hisZ | ATP phosphoribosyltransferase regulatory subunit |
Reads supporting (aligned to +/- strand): ref base T (15/15); new base A (2/2); total (17/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01 |
GGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTA > NC_002947/5559860‑5560137 | gggCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCgcgc > 2:198991/1‑150 (MQ=255) gCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCg < 2:158205/150‑1 (MQ=255) ttCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATAGCGTCGAACAGCTGCTGTTCGACCGCGCCGGa > 2:516537/1‑150 (MQ=255) tCAGCCAGTACTTCCCGACCACCACCCAGCTCGACCAGGGCTCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGAc < 1:274214/150‑1 (MQ=255) ttCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCAAGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCggcggc > 2:474323/1‑150 (MQ=255) accacACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAg > 1:363964/1‑150 (MQ=255) ccacACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAgg > 1:260142/1‑150 (MQ=255) cacacaGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGc > 1:420870/1‑150 (MQ=255) gcgcAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATg > 1:543572/1‑150 (MQ=255) cagcaTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGACTGCACTTCATCGACAGCCTTGCGTTGAATGGAGTGGAACAGCTGCTGTTAAACAGCACTTGCTGGACCGGCGTGACTGGCCAGGCCGCGGTAGATGCCGACATGGcc < 1:322995/150‑1 (MQ=255) caTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAg < 2:272273/150‑1 (MQ=255) caTGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAg < 1:821089/150‑1 (MQ=255) tgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCa > 1:492492/1‑117 (MQ=255) tgcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCa < 2:492492/117‑1 (MQ=255) gcCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAt < 2:363964/150‑1 (MQ=255) ggTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGaa < 2:673576/150‑1 (MQ=255) cGGCAGGTCGGCCGGCAGCGCCTGCACTTCATCGACAGCCGTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCGGACATGGCCCAGGTCCATGTGAACATCCgg < 1:111561/150‑1 (MQ=255) gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt > 1:373991/1‑138 (MQ=255) gCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt < 2:373991/138‑1 (MQ=255) cgcCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAg < 1:590963/150‑1 (MQ=255) ccTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCg > 1:423467/1‑150 (MQ=255) ttCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAg < 2:260142/150‑1 (MQ=255) cATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCAAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCa < 1:198991/150‑1 (MQ=255) gACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAg > 2:434340/1‑150 (MQ=255) gTTGCACGGCGTTGCACAGCTGGTGTTTGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACtt < 2:266249/150‑1 (MQ=255) gCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGa < 1:434340/150‑1 (MQ=255) aTGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGCCAGACCGGCGGCCCGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACa > 2:223126/1‑150 (MQ=255) ggCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATc > 2:111915/1‑150 (MQ=255) gTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga > 2:847864/1‑134 (MQ=255) gTCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCtga < 1:847864/134‑1 (MQ=255) tCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt > 2:535201/1‑85 (MQ=255) tCGAACAGCAGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCAtgt < 1:535201/85‑1 (MQ=255) aaCAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTa > 1:815589/1‑150 (MQ=255) aaCAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGACAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTa > 1:57333/1‑150 (MQ=255) | GGGCTTCAGCCAGTACTTCCCGACCACCACACAGCTCGACCAGGGCGCGCAGCATGCTGCCCAGGTCTTTCGGCAGGTCGGCCGTCAGCGCCTGCACTTCATCGACAGCCTTGCGTTGCATGGCGTCGAACAGCTGCTGTTCGACCGCGCCGGACAGACCGGCGGCACGGGCCAGGCCGCGGTAGATGCCGACATGGCCCAGGTCCATGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTA > NC_002947/5559860‑5560137 |
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: ‑ |