| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,558,426 | T→G | 11.8% | K293T (AAG→ACG) | purA ← | adenylosuccinate synthetase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,558,426 | 0 | T | G | 11.8% | 77.0 / 5.6 | 34 | K293T (AAG→ACG) | purA | adenylosuccinate synthetase |
| Reads supporting (aligned to +/- strand): ref base T (15/15); new base G (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 = 3.23e-01 | |||||||||||
TTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTG > NC_002947/5558283‑5558572 | tttCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCa > 2:363374/1‑150 (MQ=255) tCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCGTGGCCAAg < 2:746960/150‑1 (MQ=255) cGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgt < 1:540988/150‑1 (MQ=255) gTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtc > 2:764308/1‑150 (MQ=255) gTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGAAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtc > 2:124149/1‑150 (MQ=255) aGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAAc < 2:534733/150‑1 (MQ=255) gCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACa > 2:369464/1‑150 (MQ=255) tGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGtt < 1:369464/150‑1 (MQ=255) gCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtcgg > 2:59889/1‑150 (MQ=255) cGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACATTTCGGTCGGGAACGGAcc < 2:686391/150‑1 (MQ=255) tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGc < 1:163450/97‑1 (MQ=255) tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGc > 2:163450/1‑97 (MQ=255) gACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCAcgcg > 2:283365/1‑150 (MQ=255) gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt < 2:699005/150‑1 (MQ=255) gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt > 2:469556/1‑150 (MQ=255) acGAGGGAGCGAGACAGAATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTAGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg < 2:563820/149‑1 (MQ=255) gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt > 1:339902/1‑150 (MQ=255) gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg > 2:423878/1‑150 (MQ=255) tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt > 2:276798/1‑133 (MQ=255) tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt < 1:276798/133‑1 (MQ=255) tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATAc < 1:469556/150‑1 (MQ=255) ggCATCGACACAGCGGCTACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACGGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAg < 2:197946/150‑1 (MQ=255) aaCCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAg > 2:550922/1‑150 (MQ=255) gCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAg < 1:423878/150‑1 (MQ=255) aaCGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACa > 2:525422/1‑150 (MQ=255) acgGGCACGGCAGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCgg < 2:283507/150‑1 (MQ=255) ggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGCCGAGGTACATCGGGCCGACGCCGGAACCGGt < 1:124149/150‑1 (MQ=255) ccGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGCCGCGTCGGGTCGGAACCGTACCGTAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACCTCGGGCCGACGCCGGAACCGGTGGaa > 1:16270/1‑150 (MQ=255) cGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGaaa < 1:764308/150‑1 (MQ=255) cttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccac > 1:214770/3‑150 (MQ=255) cGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccac < 1:525422/150‑1 (MQ=255) gCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTGGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGc > 1:1064/1‑150 (MQ=255) cACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAg > 2:726285/1‑150 (MQ=255) cGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtg < 2:339902/150‑1 (MQ=255) | TTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTG > NC_002947/5558283‑5558572 |
| 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: ‑ |