Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 2,606,786 | C→A | 11.9% | H85N (CAC→AAC) | hmpA → | fused nitric oxide dioxygenase/dihydropteridine reductase 2 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,606,786 | 0 | C | A | 11.9% | 88.0 / 2.6 | 32 | H85N (CAC→AAC) | hmpA | fused nitric oxide dioxygenase/dihydropteridine reductase 2 |
Reads supporting (aligned to +/- strand): ref base C (10/18); new base A (2/2); total (12/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.20e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 |
ACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCA > REL606/2606640‑2606930 | aCCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGcaca > 1:578532/1‑150 (MQ=255) aCCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGcaca > 1:15871/1‑150 (MQ=255) aTGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGAACACCAGTTTCCAGATCAAACCGGAACAg > 1:627773/1‑150 (MQ=255) aTGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGAACACCAGTTTCCAGATCAAACCGGAACAg > 1:181799/1‑150 (MQ=255) aGTAACCAGAGTAATGGCGATCAACGTGAAGACCTGGTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAATAGTac < 1:144774/150‑1 (MQ=255) ccAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGGAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAt < 1:456323/150‑1 (MQ=255) ccAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAt < 2:296985/150‑1 (MQ=255) ccAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAt < 1:222546/150‑1 (MQ=255) ccAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAt < 2:580228/150‑1 (MQ=255) ccAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAt < 2:15922/150‑1 (MQ=255) gCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACAtcgt > 2:497671/1‑150 (MQ=255) tgtataagagacagcTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACAcc < 1:479365/136‑1 (MQ=255) aGCCCTGTTTAACGCTAATGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAAc < 1:421123/150‑1 (MQ=255) cTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACAccctgtctcttataca > 2:480724/1‑136 (MQ=255) gTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGa > 2:330439/1‑150 (MQ=255) tttAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGAc < 2:276641/150‑1 (MQ=255) aTTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTc > 2:139209/1‑150 (MQ=255) aTTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTc > 2:61963/1‑150 (MQ=255) aTTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTc > 2:288767/1‑150 (MQ=255) aTTGCCGCCTACGCCAGTAATATTGAAAACCAGCCTGCGCGGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCCATGCTGGACGAAATgata > 1:431241/1‑147 (MQ=255) cgccTAAGCCAGTAATATTGAAAACCTGCCTGCGCTTCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCACATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGccc < 1:339139/150‑1 (MQ=255) gccTACGCCAGTAATATTGAAAACCTTCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCg > 1:249905/1‑150 (MQ=255) cgcttGACAACCTGCCCTGGCAGCGGCCCGCGATCGAAAAAATCGTGCATATGCCAACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGATGAACACCTGTTGGCAACGCTGGACGAAATGTTTAGCCCGGGGCAGGAAGTGc < 1:178096/147‑1 (MQ=255) tatTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAACCCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGAGCt > 1:95287/1‑150 (MQ=255) ttGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGAACAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTgg > 1:336702/1‑150 (MQ=255) gAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGAACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc < 2:629622/150‑1 (MQ=255) gAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGAACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc < 2:182319/150‑1 (MQ=255) cgcTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAg < 2:188769/150‑1 (MQ=255) gctgctGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGc < 2:258215/150‑1 (MQ=255) gcCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATgg < 2:454832/150‑1 (MQ=255) cGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTAc < 2:456050/150‑1 (MQ=255) cGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTAc < 1:301700/150‑1 (MQ=255) gAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCt < 2:432467/150‑1 (MQ=255) aaTCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGt < 1:496263/150‑1 (MQ=255) cAGAAGTATACCAGTTTCCAGATCAAACCGGAACAGTACATCTTCGTCGGTGAACACCTGTTGGCAACGCTTGACGAAATGTTCAGCCCGGGGCAGGAAGTGCGGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATc < 1:287948/150‑1 (MQ=255) aGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCa < 1:138813/150‑1 (MQ=255) | ACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGTTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCA > REL606/2606640‑2606930 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |