| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,852,005 | G→A | 100% | R350H (CGT→CAT) | relA → | ATP:GTP 3'‑pyrophosphotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,852,005 | 0 | G | A | 94.5% | 103.8 / ‑1.8 | 36 | R350H (CGT→CAT) | relA | ATP:GTP 3'‑pyrophosphotransferase |
| Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (19/15); total (20/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGC > NC_002947/1851869‑1852150 | ggCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATg > 2:30591/1‑150 (MQ=255) gCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGc < 1:30591/150‑1 (MQ=255) aTCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTTCATATCCATACCCACGGCATGCAc < 2:814329/150‑1 (MQ=255) gcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTgg > 2:305173/1‑150 (MQ=255) ttCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgcg > 1:554229/1‑150 (MQ=255) cAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCa < 2:294864/150‑1 (MQ=255) gacgacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTAc > 2:637041/1‑150 (MQ=255) aTCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACAGTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAaga > 1:812003/1‑107 (MQ=255) aTCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAaga < 2:812003/107‑1 (MQ=255) ccccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCACCCCCGGCCCGCAGCTGGAGGTGCAGAGCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCaa > 1:105760/1‑150 (MQ=255) cccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACt > 1:845789/1‑107 (MQ=255) cccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACt < 2:845789/107‑1 (MQ=255) cccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAg > 1:606129/1‑150 (MQ=255) cAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc > 1:329434/1‑150 (MQ=255) aaGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACTGCAGTCGTTCACGAAGAGGCCGAACTTGGCGTATGCGCCCAATGGAGGCACACGTGGAAAGACGACAAGccc > 1:53356/1‑150 (MQ=255) aGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCa > 2:705491/1‑150 (MQ=255) aGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCa > 1:35516/1‑150 (MQ=255) gagaACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGc > 1:447151/1‑150 (MQ=255) tACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc < 2:554229/150‑1 (MQ=255) cTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTAt > 2:243239/1‑95 (MQ=255) cTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTAt < 1:243239/95‑1 (MQ=255) acacTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGa > 1:80993/1‑150 (MQ=255) acTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGAtt > 1:634042/1‑150 (MQ=255) cTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGACACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGAttt > 2:670751/1‑150 (MQ=255) aTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTg > 2:788661/1‑150 (MQ=255) cccGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAg < 2:80993/150‑1 (MQ=255) cAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGa > 2:33132/1‑150 (MQ=255) aGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGt < 2:105760/150‑1 (MQ=255) aGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGt < 1:305173/150‑1 (MQ=255) gTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTgg < 2:329434/150‑1 (MQ=255) cTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCAc < 2:447151/150‑1 (MQ=255) cTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCAc < 1:637041/150‑1 (MQ=255) ggTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAaga < 1:788661/150‑1 (MQ=255) gTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGaa > 1:760488/1‑147 (MQ=255) gTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGaa < 2:760488/147‑1 (MQ=255) aTCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGc < 2:634042/150‑1 (MQ=255) | GGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGC > NC_002947/1851869‑1852150 |
| 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: ‑ |