| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,256,423 | A→C | 12.8% | E53A (GAG→GCG) ‡ | PP_4633 → | DUF1244 domain‑containing protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,256,423 | 0 | A | C | 12.8% | 87.0 / 9.2 | 39 | E53A (GAG→GCG) ‡ | PP_4633 | DUF1244 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base A (18/16); new base C (3/2); total (21/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01 | |||||||||||
TCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGAAAAAGGAAAACCTCGTGACTGA > NC_002947/5256285‑5256562 | tcgagctcgagGCCGCCCCGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTgg < 1:568028/150‑1 (MQ=255) cgccgcGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGga < 2:316521/150‑1 (MQ=255) gccgcGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGgat > 1:641541/1‑150 (MQ=255) ccgcGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGgatg > 1:269439/1‑150 (MQ=255) ggCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGAc > 2:738825/1‑128 (MQ=255) ggCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGAc < 1:738825/128‑1 (MQ=255) cacCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCACGTGGTACAAGGCCGCCGCGGGCGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGc < 2:506211/150‑1 (MQ=255) acCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGc < 2:702127/102‑1 (MQ=255) acCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGc > 1:702127/1‑102 (MQ=255) acCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCa > 2:579376/1‑150 (MQ=255) cATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATa < 1:579376/150‑1 (MQ=255) aTCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCa < 1:337500/146‑1 (MQ=255) aTCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCa > 2:337500/1‑146 (MQ=255) tCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATAcc > 2:185604/1‑150 (MQ=255) aCCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGCTCGTCAg > 1:304554/1‑74 (MQ=255) aCCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGCTCGTCAg < 2:304554/74‑1 (MQ=255) atgaACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGaaa < 2:42921/132‑1 (MQ=255) atgaACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGaaa > 1:42921/1‑132 (MQ=255) tgaACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGt < 2:460577/126‑1 (MQ=255) tgaACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGt > 1:460577/1‑126 (MQ=255) gaACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGc < 1:185604/150‑1 (MQ=255) aCCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCa < 2:641541/150‑1 (MQ=255) cGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGa > 2:281494/1‑150 (MQ=255) cGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGa > 1:625180/1‑150 (MQ=255) ggTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGaa < 2:269439/150‑1 (MQ=255) ggTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGaa > 1:251657/1‑150 (MQ=255) tGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAAcaacaa > 2:479524/1‑150 (MQ=255) cAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGa > 2:587657/1‑144 (MQ=255) cAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGa < 1:587657/144‑1 (MQ=255) aaGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCgg > 1:4926/1‑67 (MQ=255) aaGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCgg < 2:4926/67‑1 (MQ=255) gTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGga < 2:128402/48‑1 (MQ=255) gTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGga > 1:128402/1‑48 (MQ=255) aGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAg > 1:15400/1‑54 (MQ=255) aGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAg < 2:15400/54‑1 (MQ=255) ccgccgCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAATTGGCGTACGAAAAAGGAACACCACGtaact > 2:325391/1‑150 (MQ=255) tccagcTGACGCTCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGAAAAAGGAAAACCTCGtgactg > 2:328014/5‑150 (MQ=255) tccagcTGACGCTCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGAAAAAGGAAAACCTCGtgactg > 2:86862/5‑150 (MQ=255) ccagcTGACGCTCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGAAAAAGGAAAACCTCGtgactga < 1:86862/147‑1 (MQ=255) | TCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACACCTGCAAAAGCGCACCGATGTACAGAACATCGACCTGATGAACCTGTCCGGTTTCTGCCGCAACTGCCTGTCCAAGTGGTACAAGGCCGCCGCTGACGAGCGTCAGCTGGACCTGAGCCTGGATGATGCCCGCGAAGCGGTGTACGGCATGCCATACGCCGAGTGGAAAGCCCAATACCAGAAAGAAGCCAGCGCCGAACAACAAGCGGCGTTCGAAAAAGGAAAACCTCGTGACTGA > NC_002947/5256285‑5256562 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |