Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 445,105 | A→T | 19.4% | intergenic (+78/‑210) | 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,105 | 0 | A | T | 19.4% | 56.5 / 10.5 | 31 | intergenic (+78/‑210) | PP_0367/PP_0368 | hypothetical protein/acyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (15/10); new base T (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 = 4.68e-01 |
CGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACT > NC_002947/444967‑445253 | 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 > 4:37213/1‑129 (MQ=255) gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc < 3:37213/129‑1 (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) acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt > 2:295711/1‑107 (MQ=255) acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa > 1:257373/1‑148 (MQ=255) acaacaATTTTCTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt < 1:295711/107‑1 (MQ=255) acaacaATTTTCCCGTGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGCCATGGCATGACTGATTTCCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa < 8:221795/149‑1 (MQ=255) aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa > 7:97998/1‑129 (MQ=255) aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa < 8:97998/129‑1 (MQ=255) aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa < 6:354344/140‑1 (MQ=255) aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa > 5:354344/1‑140 (MQ=255) aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa > 3:2099/1‑148 (MQ=255) aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac > 6:129744/1‑149 (MQ=255) aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac < 5:155680/149‑1 (MQ=255) aTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGCTGAGTGGCTTGGCATGACTGATTTCCCATTTTTCGTGGCGTAAATGGCACCAAAGTCCATCCTTGACAGAAACGTGAGCTAAAAGTAAGTTTCAAACAACt > 4:314428/1‑149 (MQ=255) | CGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACT > NC_002947/444967‑445253 |
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: ‑ |