Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,975,888 | A→C | 9.1% | N274H (AAC→CAC) ‡ | pheA → | bifunctional chorismate mutase/prephenate dehydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,975,888 | 0 | A | C | 9.1% | 160.9 / 7.8 | 55 | N274H (AAC→CAC) ‡ | pheA | bifunctional chorismate mutase/prephenate dehydratase |
Reads supporting (aligned to +/- strand): ref base A (26/24); new base C (3/2); total (29/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAG > NC_002947/1975742‑1976033 | aTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc < 3:199750/149‑1 (MQ=255) tGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAcc < 2:10945/149‑1 (MQ=255) tGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc > 1:47770/1‑148 (MQ=255) cGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAgg < 7:288161/149‑1 (MQ=255) gCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAcc > 7:100438/1‑149 (MQ=255) gTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCgacgac < 1:248219/149‑1 (MQ=255) tgTGCAACTCGGCGCCTTTCGCGGGCTTTATGGCGGCCACCCTGTCCGGCCTGACCCGTCTGGCCTATAATATCCAAGACCGCCCTCACAACTCCCCGCGCTTCCTGTTTAGCGTTAACCAGCAAGTACCGCCGACCGGCGACGACaa < 8:100438/147‑1 (MQ=255) tGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGAc < 2:47770/149‑1 (MQ=255) gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt > 4:164326/1‑149 (MQ=255) gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGAACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCAGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt > 4:33845/1‑149 (MQ=255) ggcggcGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGatc < 5:241256/148‑1 (MQ=255) ccGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATg < 7:8838/149‑1 (MQ=255) cGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt < 6:127976/147‑1 (MQ=255) cGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt > 5:127976/1‑147 (MQ=255) gCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGAAACGATGatc < 5:198496/149‑3 (MQ=255) atatGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGCCGACGACAAGACCTCGATCATCGTTTCGATGAGcaa > 5:249743/1‑148 (MQ=255) atcttctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCg > 6:257674/8‑149 (MQ=255) cTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGtt > 2:93993/1‑149 (MQ=255) ctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTg < 5:257674/145‑1 (MQ=255) ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGTAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGt < 3:33845/149‑1 (MQ=255) ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGGAACCAGGAAGTACCGCCGACCGGCGCCGACAAGAGCTCGATCATCGTTTGGATGAGCGACAAGCCAGGTGGGTTGCATGAGt > 6:102416/1‑149 (MQ=255) tGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGGTTACCGATCATCAGGAAgc > 5:247114/1‑75 (MQ=25) tGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGGTTACCGATCATCAGGAAgc < 6:247114/76‑2 (MQ=25) aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCt > 5:158174/1‑149 (MQ=255) ccGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTgg > 8:293437/1‑149 (MQ=255) tCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGc > 6:118150/1‑149 (MQ=255) ccGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc > 6:134148/1‑149 (MQ=255) gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATc < 6:101531/149‑1 (MQ=255) gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCAGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt > 6:671/1‑122 (MQ=255) gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCAGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt < 5:671/122‑1 (MQ=255) aagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATTATCGGTAACCAGGAATTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCa < 5:271445/149‑1 (MQ=255) aTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa < 5:134148/149‑1 (MQ=255) tCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGCCCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGATTTGCTGGTGCCGTTCTATCAGAAc > 2:153933/1‑149 (MQ=255) aaGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCa < 1:121273/107‑1 (MQ=255) aaGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCa > 2:121273/1‑107 (MQ=255) aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc < 3:160963/114‑1 (MQ=255) aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc > 4:160963/1‑114 (MQ=255) cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGCCCTCGATCATCGTTTCGAAGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAcc < 5:102416/149‑1 (MQ=255) ccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt < 3:164326/149‑1 (MQ=255) tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt > 8:54034/1‑149 (MQ=255) cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc < 6:249743/148‑1 (MQ=255) cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc < 5:118150/148‑1 (MQ=255) cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCg > 7:120611/1‑149 (MQ=255) cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAACACCTCGATCATCGTTTCGATGAGCAACAAGCCAGCTGCGTTGc > 4:104361/1‑98 (MQ=255) cggtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGAGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc > 1:62039/6‑147 (MQ=255) gtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc > 7:235462/4‑115 (MQ=38) gtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc < 8:235462/112‑1 (MQ=38) gcTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcc > 7:110341/1‑149 (MQ=255) cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTg > 1:57448/1‑102 (MQ=255) cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTg < 2:57448/102‑1 (MQ=255) cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgccc > 7:87984/1‑149 (MQ=255) cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca > 1:252814/1‑121 (MQ=255) cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca < 2:252814/124‑4 (MQ=255) tCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgcg > 4:249930/1‑148 (MQ=255) ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg > 8:132460/1‑149 (MQ=255) | ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAG > NC_002947/1975742‑1976033 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |