Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,178,335 | G→A | 10.2% | P404P (CCG→CCA) | guaA → | glutamine‑hydrolyzing GMP synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,178,335 | 0 | G | A | 10.2% | 101.2 / 2.7 | 39 | P404P (CCG→CCA) | guaA | glutamine‑hydrolyzing GMP synthase |
Reads supporting (aligned to +/- strand): ref base G (14/21); new base A (2/2); total (16/23) | |||||||||||
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.94e-01 |
GTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGGTGTTCCAGC > NC_002947/1178191‑1178480 | gTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTcccgggccc < 2:803201/150‑1 (MQ=255) cGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTcccgggcccgg < 1:1103420/150‑1 (MQ=255) aaCGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgc > 2:137551/1‑150 (MQ=255) gagtgtgcCTGCCGGAGGAGATGAACCTCCAGCTAGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCACGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgc < 2:790086/144‑1 (MQ=255) cTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCAGGCCTGGGCGTGCGTATCCTc > 1:428518/1‑150 (MQ=255) ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgcgt > 1:203634/1‑141 (MQ=255) ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGgcgtgcgt < 2:203634/141‑1 (MQ=255) ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaacgg > 2:582452/1‑134 (MQ=255) ccGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaacgg < 1:582452/139‑6 (MQ=255) aggagATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaag < 2:76121/150‑1 (MQ=255) aaCCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtgccgttt < 1:588094/150‑7 (MQ=255) ccTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGa > 1:244988/1‑150 (MQ=255) ccTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGa < 2:325978/150‑1 (MQ=255) cTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGcc > 2:392141/1‑150 (MQ=255) ccGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTg < 2:729217/150‑1 (MQ=255) tGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCgcgc > 2:1049433/1‑150 (MQ=255) cTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAgg > 1:1057202/1‑121 (MQ=255) cTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAgg < 2:1057202/121‑1 (MQ=255) ttCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACcac > 1:489183/1‑150 (MQ=255) gACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaac < 2:666173/86‑4 (MQ=255) gACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCAGGcccgggcccgggaac > 1:666173/1‑83 (MQ=255) gAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCg < 1:353403/118‑1 (MQ=255) gAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCg > 2:353403/1‑118 (MQ=255) cgcAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAagag > 1:7984/1‑150 (MQ=255) cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTc > 2:1043992/1‑85 (MQ=255) cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTc < 1:1043992/85‑1 (MQ=255) cAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCt < 2:179321/150‑1 (MQ=255) aGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTgc < 1:650536/150‑1 (MQ=255) cGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCaa < 2:7984/150‑1 (MQ=255) ccTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCAGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGc < 2:428518/150‑1 (MQ=255) aGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGTAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCGTCGAAGAGCTGCGCAAGGCCGAc < 2:554807/149‑1 (MQ=255) aGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGAc > 1:554807/1‑149 (MQ=255) tGGGCATGTGATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGt < 1:326153/150‑1 (MQ=255) tCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtg < 1:1071259/56‑1 (MQ=255) tCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGgtgaagtg > 2:1071259/1‑56 (MQ=255) cTACCCCCATCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGtt < 1:46854/150‑1 (MQ=255) gCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCgtgg < 1:137551/150‑1 (MQ=255) aCCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGgtgt < 2:489183/150‑1 (MQ=255) cccGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGc > 1:893433/1‑138 (MQ=255) cccGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGc < 2:893433/138‑1 (MQ=255) tCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAATGCCGACTGGTACCTCAAGACCATCCAGGCGTTCGTGGTGTTCCAGc < 2:244988/150‑1 (MQ=255) | GTCGCACCACAACGTAGGCGGCCTGCCGGAGGAGATGAACCTCAAGCTGGTCGAGCCGCTGCGTGAACTGTTCAAGGACGAAGTGCGCAAGATCGGCCTGGAGCTGGGCCTGCCATACGACATGGTCTACCGCCACCCGTTCCCGGGCCCGGGCCTGGGCGTGCGTATCCTCGGTGAAGTGAAGAAGGAATACGCCGACATCCTGCGTCGCGCCGACCACATCTTCATCGAAGAGCTGCGCAAGGCCGACTGGTACCACAAGACCAGCCAGGCGTTCGTGGTGTTCCAGC > NC_002947/1178191‑1178480 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |