Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 445,095 | T→C | 19.3% | intergenic (+68/‑220) | PP_0367 → / → PP_0368 | hypothetical protein/acyl‑CoA dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 445,095 | 0 | T | C | 19.3% | 58.8 / 8.7 | 31 | intergenic (+68/‑220) | PP_0367/PP_0368 | hypothetical protein/acyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (15/10); new base C (3/3); total (18/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01 |
GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAG > NC_002947/444949‑445241 | gACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa < 1:326045/149‑1 (MQ=255) ccGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa < 5:115800/149‑1 (MQ=255) aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGaca > 5:333242/1‑148 (MQ=255) aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGacaaaattgt > 1:156885/1‑143 (MQ=255) cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTga < 2:156885/149‑2 (MQ=255) gggcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc < 7:247930/146‑1 (MQ=255) gctcggcAAGGCATCCGACGAAATGCTTCGCCGCATCATCGACACCTACGCCTGAATACTTGCGGCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc < 8:298974/147‑1 (MQ=255) aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc < 6:333242/149‑1 (MQ=255) gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc < 3:37213/129‑1 (MQ=255) gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc > 4:37213/1‑129 (MQ=255) tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCGCATTCGTACGCTAAAACAACAATTTTTTAAGGTCGGGCCTCTTAGCGGGGC‑AGCCCGCGCCAAACAGGGGTtg > 7:56890/1‑149 (MQ=255) cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGCATGCGTCCGCTCAGACAACAATTTTTTGTGGTGGGGCCTAATAGCCGGCA‑AGCCCGCGCCGACAAGGGGTTGTGGGATGGCATGAGTGATTTTGCATTTTTAGTgg > 6:155680/1‑149 (MQ=255) cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc < 1:70690/86‑1 (MQ=255) cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc > 2:70690/1‑86 (MQ=255) accttCAGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGTCCTACCACAAAACAATTTTGTCTGGCCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGt < 7:334441/146‑1 (MQ=38) gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGAAAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt < 1:228687/142‑1 (MQ=255) gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCGATCCTGGa > 7:221795/1‑146 (MQ=255) cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc > 7:150484/1‑149 (MQ=255) gcgcgAATGCGTACGATCAGACAACAATTTTGTATGGTCGGGCCTCTTCGCGTGCA‑AGCGCGCTCCTTCAAGGGTTTGTTGCTTGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGTCTCCTTAGTCCATCCTTGa < 4:1761/139‑1 (MQ=255) ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAGAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc > 1:98631/1‑149 (MQ=255) ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGTCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGACCATCATTGACAGATAgg > 3:1761/1‑147 (MQ=255) cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt > 5:28086/1‑148 (MQ=255) cGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc > 8:38742/1‑149 (MQ=255) gcccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt < 5:252000/86‑1 (MQ=18) gcacGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt > 6:252000/4‑88 (MQ=18) gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg > 8:71125/1‑149 (MQ=255) ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg > 2:46324/1‑140 (MQ=255) ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg < 1:46324/140‑1 (MQ=255) ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGTCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt > 6:106845/1‑148 (MQ=38) gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg > 4:275588/1‑149 (MQ=255) gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATCGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg > 6:364900/1‑149 (MQ=255) | GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAG > NC_002947/444949‑445241 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |