| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,413,542 | A→T | 100% | Y67N (TAT→AAT) | ydaC ← | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,413,542 | 0 | A | T | 94.5% | 50.2 / ‑3.7 | 18 | Y67N (TAT→AAT) | ydaC | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
| Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (10/7); minor base C (0/1); total (10/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT > NC_000913/1413401‑1413684 | gcctCAGGCCGTATTCAACCATCCCCTCTTCATTAAAGCTTACTTGCTCCATCGCATCACCTCTGGCAGGCGGCAATGGTAAACCACTGTTTATTCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCa < 1:681161/146‑1 (MQ=255) gCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCg > 2:438771/1‑149 (MQ=255) tCCACTCTTCATTAACGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAATTTTAATTTGAATTCAATGCAATTCTCGAGGACTGATGTTtc > 2:52039/1‑149 (MQ=255) tCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAtt > 1:131132/1‑149 (MQ=255) aTTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaa > 1:42869/1‑149 (MQ=255) ttAAAGATTACTTGTGCCAGCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa > 1:502502/1‑149 (MQ=255) tgtgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGttt < 2:184742/149‑1 (MQ=255) gtgCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAACTTTCTCGCAATTAAAATTTATCAGtttt < 2:203518/149‑1 (MQ=255) gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc > 2:573267/1‑149 (MQ=255) cAATGTTAGACTGAAATTGACGCCCGATGATGATTATAAATAATCAGCTATGAAGTTTTAATTAGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGctt > 1:448860/1‑149 (MQ=255) aaTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTgag < 2:539823/149‑1 (MQ=255) ccGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGGGAGCAttttt > 1:98642/1‑149 (MQ=255) gATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc < 1:376369/149‑1 (MQ=255) attattAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCt < 1:66041/149‑1 (MQ=255) cAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGAATGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGATCTCTGGAAACGCCTGCTTTTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGt < 2:42869/149‑1 (MQ=255) ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTCACATTTATCAGCCTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGGTTTTGa > 2:594068/1‑149 (MQ=255) ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCACTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGttttgattttga > 2:251981/1‑149 (MQ=255) ttGAATCCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCtt < 2:61954/149‑1 (MQ=255) | GCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT > NC_000913/1413401‑1413684 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |