Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,975,678 | T→A | 10.6% | intergenic (+622/‑522) | PP_5727 → / → yifB | hypothetical protein/ATP‑dependent chelatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,975,678 | 0 | T | A | 10.6% | 118.5 / 5.5 | 47 | intergenic (+622/‑522) | PP_5727/yifB | hypothetical protein/ATP‑dependent chelatase |
Reads supporting (aligned to +/- strand): ref base T (17/25); new base A (3/2); total (20/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01 |
TTTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAT > NC_002947/5975533‑5975825 | tttGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATaa < 5:58577/148‑1 (MQ=255) tttGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATaa < 6:159445/148‑1 (MQ=255) aaaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAGAAGTGAGGAATTATAAATATTGAACCTGGAAAAGATACGATGCGTATAAGCGGGCCATACGAATACCTGCAAAACATATATAatattt > 1:98434/1‑147 (MQ=255) aaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt < 2:236487/149‑1 (MQ=255) gCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTaa < 2:225320/149‑1 (MQ=255) cGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAg < 6:65608/149‑1 (MQ=255) ggaggaTTCTCGAGAGAGTAAAATAAATTCAGCATGCTCCATCAATTCAAAAGTGA‑‑AATTATACATATTGAACCTGGAAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt < 3:3907/149‑1 (MQ=255) ggaggaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt < 1:104670/149‑1 (MQ=255) gaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCt > 4:301011/1‑149 (MQ=255) gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGt < 7:141026/148‑1 (MQ=255) gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg > 2:198260/1‑149 (MQ=255) agagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTgg < 2:120732/149‑1 (MQ=255) ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg > 6:199581/1‑126 (MQ=255) ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg < 5:199581/127‑2 (MQ=255) ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa < 2:157339/149‑1 (MQ=255) ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa > 3:73641/1‑149 (MQ=255) aataaatTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca < 1:198260/149‑1 (MQ=255) aataaatTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca > 3:176782/1‑149 (MQ=255) ttCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt > 4:201276/1‑112 (MQ=255) ttCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt < 3:201276/112‑1 (MQ=255) gCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGcc > 8:65769/1‑149 (MQ=255) aTGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCg > 6:97750/1‑148 (MQ=255) aaTGCACAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATCCGTATAAGCGGGCCATAAGAAGACCTGCAAAACATATATAAGATATTTTGAATTTTAGTGCCTCGTAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatc < 2:98434/148‑1 (MQ=255) tGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc < 4:176782/149‑1 (MQ=255) gCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcat > 2:250155/1‑147 (MQ=255) cATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTg < 3:45394/148‑1 (MQ=255) atTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAACACGAGTGCCa > 6:7190/1‑148 (MQ=255) tGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc > 8:157260/1‑120 (MQ=255) tGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc < 7:157260/120‑1 (MQ=255) gCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagct < 1:250155/149‑1 (MQ=255) gCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagc > 5:275016/1‑148 (MQ=255) aaaGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagctag < 3:301011/149‑1 (MQ=255) gATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg < 8:284870/54‑1 (MQ=255) gATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg > 7:284870/1‑54 (MQ=255) aTGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACt > 6:17155/1‑149 (MQ=255) cGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGAt < 5:97750/149‑1 (MQ=255) cGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAACACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGAt < 5:7190/149‑1 (MQ=255) cGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgt < 6:275016/149‑1 (MQ=255) gCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg > 3:258801/1‑147 (MQ=255) gCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg < 4:258801/147‑1 (MQ=255) ccATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc > 5:234379/1‑149 (MQ=255) aCGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCgg > 8:148530/1‑148 (MQ=255) ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta > 1:70729/5‑116 (MQ=38) ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta < 2:70729/112‑1 (MQ=38) tGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACg > 5:280327/1‑149 (MQ=255) aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTa < 7:65769/149‑1 (MQ=255) aCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTaa < 7:148530/149‑1 (MQ=255) ataAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt > 7:267000/1‑149 (MQ=255) | TTTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAT > NC_002947/5975533‑5975825 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |