| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 6,095,227 | C→A | 15.4% | D34Y (GAC→TAC) | pycB ← | pyruvate carboxylase subunit B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 6,095,227 | 0 | C | A | 15.4% | 61.0 / 3.0 | 26 | D34Y (GAC→TAC) | pycB | pyruvate carboxylase subunit B |
| Reads supporting (aligned to +/- strand): ref base C (14/8); new base A (2/2); total (16/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.25e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01 | |||||||||||
CAGCATTTGCAGGCGGGTGTTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAGGCCTGCGTGGGCGGCGATGGCGGCGGCGATGG > NC_002947/6095081‑6095375 | cagcaTTTGCAGGCGGGTGTTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCa < 1:737186/150‑1 (MQ=255) aTTTGCAGGCGGGTGTTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATc < 2:75925/150‑1 (MQ=255) tttGCAGGCGGGTGTTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCg > 2:147180/1‑150 (MQ=255) gtTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTc > 2:429183/1‑150 (MQ=255) gcgcTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTg > 2:126875/1‑150 (MQ=255) tCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGa > 1:119833/1‑150 (MQ=255) ccACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTAGCa < 2:71015/96‑1 (MQ=255) ccACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTAGCa < 2:71002/96‑1 (MQ=255) ccACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTAGCa > 1:71002/1‑96 (MQ=255) ccACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTAGCa > 1:71015/1‑96 (MQ=255) tCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCgcgc > 1:85729/1‑150 (MQ=255) ttGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATc < 1:147180/150‑1 (MQ=255) gAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTcgc < 1:72100/137‑1 (MQ=255) gAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTcgc > 2:72100/1‑137 (MQ=255) gAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTc < 1:503719/150‑1 (MQ=255) gAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGAATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCtt > 1:5493/1‑150 (MQ=255) ccaccCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATtata > 1:340419/1‑150 (MQ=255) caccCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAg > 1:519918/1‑150 (MQ=255) tCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTcaca < 2:85729/150‑1 (MQ=255) cAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAgg > 2:223885/1‑150 (MQ=255) gCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAGGCCTGCGTGggcg > 2:214570/1‑150 (MQ=255) ccTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCTGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAGGCCTGCGTGggcggcgat > 2:487016/1‑150 (MQ=255) gAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAAtttt > 2:280924/1‑99 (MQ=255) gAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAAtttt < 1:280924/99‑1 (MQ=255) tGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGTCATTATAGGGTTCCTCACAGGCCTGCGTGggcggcgatggcggcggcgat < 1:429183/150‑1 (MQ=255) tCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGGGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAGGCCTGCGTGggcggcgatggcggcggcgatgg > 1:427647/1‑150 (MQ=255) | CAGCATTTGCAGGCGGGTGTTTGGCAGCGCTGCGCGCAGCTTGCGCAGGCGCTCCCACGGGTCTTCCTTGAGGAAGCGCACGCAGGCGTCGAAGGTTGCACCACCCCAGACTTCCAGCGACCAGTAGCCAACCTTGTCGAGCTTGTCGCAGATCGGCAGCATGTCTTCGGTACGCATGCGGGTAGCCAGCAGGGACTGGTGAGCGTCGCGCAGGATCGTGTCGGTTACGTGAATTTTCTTGGACATTATAGGGTTCCTCACAGGCCTGCGTGGGCGGCGATGGCGGCGGCGATGG > NC_002947/6095081‑6095375 |
| 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: ‑ |