| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,800,887 | C→T | 26.3% | 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.3% | 25.1 / 7.8 | 19 | intergenic (+54/+25) | PP_3360/PP_3361 | membrane protein/aminotransferase |
| Reads supporting (aligned to +/- strand): ref base C (5/9); new base T (3/2); total (8/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.03e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCGTGCCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGA > NC_002947/3800739‑3801033 | ccGTGCCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCCGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTTGCCAGATAGGTGCACCGCAGGGCGTGCACATgcg > 1:84615/1‑150 (MQ=255) gTGCCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCa > 1:205604/1‑150 (MQ=255) tGCCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCAc < 2:113489/150‑1 (MQ=255) gcgcgcGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGGTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGc < 2:171768/150‑1 (MQ=255) gATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCgcgc < 1:411390/150‑1 (MQ=255) tGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAgc < 1:279961/150‑1 (MQ=255) tCCGGCCCGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTTGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGaa < 2:84615/150‑1 (MQ=255) gAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCTCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCaagaag < 2:205604/150‑1 (MQ=255) gAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCaagaag < 1:404358/150‑1 (MQ=255) cTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCAtttt < 2:355429/150‑1 (MQ=255) ccAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTc < 1:460389/63‑1 (MQ=255) ccAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTc > 2:460389/1‑63 (MQ=255) tGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCa > 1:33136/1‑150 (MQ=255) gCAGGGCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCggg > 2:154576/1‑150 (MQ=255) gCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTAGTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCggg > 1:393491/1‑150 (MQ=255) ggCGTGCGCATGTGCACGCCCTGCGCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCtt > 1:497534/1‑150 (MQ=255) gcaTGTGCACGCCCTGCGCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGcc < 2:497534/149‑1 (MQ=255) aTGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGc > 1:509367/1‑93 (MQ=255) aTGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGc < 2:509367/93‑1 (MQ=255) aTGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGa > 2:432290/1‑150 (MQ=255) | CCGTGCCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGA > NC_002947/3800739‑3801033 |
| 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: ‑ |