| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 955,600 | G→T | 18.2% | R80L (CGC→CTC) | alaC → | aminotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 955,600 | 0 | G | T | 18.2% | 55.0 / 2.4 | 22 | R80L (CGC→CTC) | alaC | aminotransferase |
| Reads supporting (aligned to +/- strand): ref base G (13/5); new base T (2/2); total (15/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01 | |||||||||||
ATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACC > NC_002947/955476‑955722 | atcGACCTGAGCATGGGCAACACCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCAtct > 2:270905/1‑136 (MQ=255) tcGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCAtctc < 1:372272/136‑1 (MQ=255) aTGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACa > 2:272142/1‑136 (MQ=255) aTGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGACCAGCGTGAAGACGCCCACGGCTACTCCACGTCGCGCGGCCTCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACa > 2:271696/1‑136 (MQ=255) ccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTa > 2:29761/1‑136 (MQ=255) gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCTCCGGGCCAtctc < 1:59638/109‑1 (MQ=255) gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCTCCGGGCCAtctc > 2:125156/1‑109 (MQ=255) gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCTCCGGGCCAtctc > 2:59638/1‑109 (MQ=255) gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCTCCGGGCCAtctc < 1:125156/109‑1 (MQ=255) gcgcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGt > 1:196255/1‑136 (MQ=255) cacaTCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCgg > 1:163840/1‑136 (MQ=255) gTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCa < 1:192935/73‑1 (MQ=255) gTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCa > 2:192935/1‑73 (MQ=255) tACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCGACGTGGCGCGGCATCCCTCGCCTGCGCCGGGCCATCTCCATCTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCAc < 1:56110/136‑1 (MQ=255) cAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCa > 2:10827/1‑136 (MQ=255) cAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCa > 2:11333/1‑136 (MQ=255) cACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGccc < 1:60163/136‑1 (MQ=255) aCGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCa > 1:148386/1‑136 (MQ=255) ccACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCt > 1:8970/1‑136 (MQ=255) cgcgcgGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCAc > 1:349968/1‑136 (MQ=255) cgcgcgGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCAc > 2:406991/1‑136 (MQ=255) aTCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGAcc < 1:270905/136‑1 (MQ=255) | ATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACC > NC_002947/955476‑955722 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |