Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,800,887 | C→T | 26.9% | intergenic (+54/+25) | PP_3360 → / ← PP_3361 | membrane protein/aminotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,800,887 | 0 | C | T | 26.9% | 27.7 / 7.8 | 22 | intergenic (+54/+25) | PP_3360/PP_3361 | membrane protein/aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (10/6); new base T (3/3); total (13/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 |
CAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGC > NC_002947/3800775‑3801030 | cAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCGCATGTGCACGCCCTGcggtg < 1:293036/128‑5 (MQ=255) cAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCGCATGTGCACGCCCTGcggtg > 2:293036/1‑124 (MQ=255) tCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCaa > 1:515483/1‑150 (MQ=255) ggggCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTg > 1:216596/1‑147 (MQ=255) ggggCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTg < 2:216596/147‑1 (MQ=255) ggggCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCCCGCTGAAACGCTATTTCAGCGCCGCGGCGAAGCTGTCTTTCAGCACCTGAAGAAGGAAGTCCGCATTGGcc > 1:376534/1‑150 (MQ=255) cTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGGCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGCCTTTCAGCACCTCAAGAAGGAAGTCGGCATTGCCCCGGATg > 2:129797/1‑150 (MQ=255) aTGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTAAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa > 1:76145/1‑150 (MQ=255) tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa < 2:301473/149‑1 (MQ=255) tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa > 1:301473/1‑149 (MQ=255) cTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACACGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAatca < 1:539463/150‑1 (MQ=255) ggCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAatcatc < 1:591144/150‑1 (MQ=255) gCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCg > 2:145704/1‑150 (MQ=255) cccGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGa > 2:47919/1‑85 (MQ=255) cccGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGa < 1:47919/85‑1 (MQ=255) cGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCa < 1:485972/150‑1 (MQ=255) aGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCAc > 1:89466/1‑150 (MQ=255) cGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACACTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGAATTTTCAGCACGTTGTCATGCgg > 1:162434/1‑150 (MQ=255) cAGGGCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGc > 1:132281/1‑150 (MQ=255) cgcaTGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGc < 2:132281/148‑1 (MQ=255) gcaTGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaaggaa < 1:350262/77‑1 (MQ=255) gcaTGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaaggaa > 2:350262/2‑78 (MQ=255) | CAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGC > NC_002947/3800775‑3801030 |
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: ‑ |