| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,548,571 | 0 | G | T | 31.2% | 15.7 / 10.5 | 16 | T500T (ACC→ACA) | malQ | 4‑alpha‑glucanotransferase (amylomaltase) |
| Reads supporting (aligned to +/- strand): ref base G (10/1); new base T (0/5); total (10/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.37e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAA > NC_000913/3548445‑3548706 | cGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCAc < 7:137890/143‑1 (MQ=255)cGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCAc > 8:137890/1‑143 (MQ=255) tttCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCAc > 6:161613/1‑119 (MQ=255) gTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAAACACCATACAGCGATGACGTTTACTTTCGAGTGCCAGCATCGAGagaaga > 7:55063/1‑148 (MQ=255) tCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGATGCGCAGCTTACCGACAATCTCTACCGGTACGGTAACGAGATCTTCACCAATCACCATACAGCGATGCAGTTTACTTTCGAGTGCCAGAATCTAGAGAAGatcagg > 7:262/1‑147 (MQ=255) aTAGAGCACTTTGTAAGAGTACACACCTCTGCTGCGCAGATTTTTGTTTTTATTTATAGGGAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGt < 8:55063/149‑1 (MQ=255) tttGTAAGAGTACACACCGCTGCTGCGCAGCTTATTTTTAATATTTATAGGTAGTGTACCGAGATCTTCACCAATCAc < 2:69639/78‑1 (MQ=255) tttGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCAc > 1:69639/1‑78 (MQ=255) aaGAGTACACACCGCTGCTGCGCAGCTTATTTTTTTTATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCg < 2:28962/82‑1 (MQ=255) aaGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCg > 1:28962/1‑82 (MQ=255) acacCGCTGCTGCGCAGCTTATTTTTAATATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGGCCTGATCTGc < 4:10851/149‑1 (MQ=255) acacCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTAAACAATCACAATACAGAGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATACACCGGATAGTGAACATACGCTCCCTGATCAGc > 3:10851/1‑149 (MQ=255) gctgctGCGCAGCTTATTTATTATATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc < 8:124489/114‑1 (MQ=255) gctgctGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc > 7:124489/1‑114 (MQ=255) tctcTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACaaa > 2:84326/1‑149 (MQ=255) tACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTc > 5:133805/1‑57 (MQ=255) | CGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAA > NC_000913/3548445‑3548706 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |