Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,191,676 | C→T | 27.4% | H366H (CAC→CAT) | icd → | e14 prophage; isocitrate dehydrogenase, specific for NADP+ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 1,191,676 | 0 | C | T | 27.4% | 25.4 / 9.0 | 22 | H366H (CAC→CAT) | icd | e14 prophage; isocitrate dehydrogenase, specific for NADP+ |
Reads supporting (aligned to +/- strand): ref base C (12/4); new base T (3/3); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.34e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01 |
ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATC > CP009273/1191536‑1191814 | atcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGGCGCCACATga > 1:370486/1‑144 (MQ=255) gtcccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGActg > 1:283408/3‑143 (MQ=255) aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACACAGCAAATCCTGGCTCTATTATTCTCCCCGCTGAGATCATGCTGCGCCACATGGGTTGGCCCGAAACGGccgcgc > 1:28896/1‑140 (MQ=255) gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg > 2:394352/1‑145 (MQ=255) gcgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt < 1:285087/145‑1 (MQ=255) tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCCTCTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAAGTATGGAAGtcg > 2:123017/1‑145 (MQ=37) aGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGGTTGTTAAAGGTATGGAAGGCGCAATc > 2:404699/1‑145 (MQ=25) cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTACGGAAGGCGCAATCAAtgc < 2:256749/145‑1 (MQ=25) cccacGGTACTGCGCCGAAATATGCCGGAAAGGACAAAGTAAATCATGGCTCTAATATACACTCCGATGAGATGATGCTGTGCCACATGGGTTGGACCGAAGCAGCTGGCTGAATTGTTAAAGGTAAGGAAGGGGCAATCAAcgc > 1:84701/1‑145 (MQ=255) gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGACGCAATCAACGCGAAAACCGTAAcc > 2:42219/1‑145 (MQ=255) atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCCACGCGAAAACCGTAACCTATGGCtt > 1:404552/1‑145 (MQ=255) tGCTGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCCCTCCGCTGGGATGATGCTGCGCCATATGGGTGGGACGGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc < 2:178446/145‑1 (MQ=14) tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc < 2:175399/145‑1 (MQ=14) tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGATAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTc > 2:124860/1‑145 (MQ=255) gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga > 2:249060/1‑145 (MQ=255) gTAAATCCTGGCTCTATTATTCTCCCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggatg < 1:265115/145‑1 (MQ=255) aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGgc < 1:394352/145‑1 (MQ=255) ctgttattTTCTCCGCTGATATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgct > 1:419871/4‑145 (MQ=255) ctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGACGGCCCAATCATTGGTACAACCAAAACCTATAAATTCGAAAATCTGATGGATGGCGCTAAAgtgct > 2:185727/1‑145 (MQ=255) ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg > 2:418634/1‑145 (MQ=255) tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGCCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt > 2:176029/1‑145 (MQ=255) aggcggcgcgtcaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg < 2:28724/134‑1 (MQ=255) gccccACAAGGGTTGCACCGAAGCGGCCGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatc < 2:91034/142‑1 (MQ=255) | ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATC > CP009273/1191536‑1191814 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |