Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,173,361 | Δ2 bp | 100% | pseudogene (917‑918/1358 nt) | gatC ← | galactitol‑specific PTS enzyme IIC component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,173,361 | 0 | C | . | 81.8% | 22.1 / 2.9 | 11 | pseudogene (918/1358 nt) | gatC | galactitol‑specific PTS enzyme IIC component |
Reads supporting (aligned to +/- strand): ref base C (0/2); new base . (5/4); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
* | NC_000913 | 2,173,362 | 0 | C | . | 81.8% | 24.2 / 0.8 | 11 | pseudogene (917/1358 nt) | gatC | galactitol‑specific PTS enzyme IIC component |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (5/4); minor base G (0/1); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA > NC_000913/2173225‑2173491 || ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCTAAACGGTAGCATCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt < 1:237824/143‑6 (MQ=255) cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa < 2:123050/143‑1 (MQ=255) gcGAGGAAGGCGATGGTGGCGAGATCGCCAAACGGCGGCACCTGATTACCGGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGTTACCA‑‑CCGCCGTATGTCCTAGCAGCAACGCCGGATCAAg < 2:125188/142‑1 (MQ=255) gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATAGTGAGAGGAATAAAAAACAGGATTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGACGGAGCAAGGCCAATCAGGAACTCCTg > 1:40774/1‑143 (MQ=255) tGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGTGATAAAATTCAGTTTTTCTTTTAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg < 1:410030/143‑1 (MQ=255) tGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGACGAATAACCCAAAAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg > 1:286225/1‑143 (MQ=255) tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGTTCCCA‑‑CCTCCGTATGTCCCAGCAGCAACGGCGGATCAAGGCAAATAAGGAACTCCTAACCGCCGAACTTCGACTGTAAAca > 1:299439/1‑142 (MQ=255) cccGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGAGAGCAGAAACCA‑‑CCGAAGTATGGCCCAGCAGCAACGCGGGATCAAGGCGACTCAAGCAATGCTCGCCGCCGAACCGACCCTGTAAAcacgtacg > 2:347517/1‑136 (MQ=255) cccacaAACTGCAATTAAAATGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt < 2:163795/141‑1 (MQ=255) tcgggggATAAAAAGCAGGCTTGGCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAAGGACGGCGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat < 2:49751/139‑1 (MQ=255) tgtCAGGCGTGCCGATGCCA‑‑CCGCCGTAGGTCCCAGCAGCAACGGCGGATCGATGTGAATCAGGAGCTCCGGACCGCCGAAGTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCATGATGGGtt < 2:312975/141‑1 (MQ=255) gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGATAAACCATCCATGATGGGTTTAATCa > 2:166584/1‑143 (MQ=255) || CCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA > NC_000913/2173225‑2173491 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |