| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 512,643 | A→G | 53.5% | intergenic (+711/‑150) | recR → / → rrs1 | recombination protein RecR/16S ribosomal RNA |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 512,643 | 0 | A | G | 53.5% | 12.1 / 68.7 | 58 | intergenic (+711/‑150) | recR/rrs1 | recombination protein RecR/16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (11/16); new base G (21/10); total (32/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.34e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT > CP000730/512503‑512781 | aGGTGTAATTTTTACTGATGTTCAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta < 1:204391/141‑1 (MQ=255)aGGTGTAATTTTTAATACTGTAACTCAGTGTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATAACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta < 1:72606/141‑1 (MQ=255) agaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttactcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACATAAATGTtgata < 2:243147‑M1/27‑1 (MQ=255) atttttatttttTTAAAACAGTTTGAAATTCTGTATATGATTTACTTTTTGAGTTAAAAAAAACAACTTAACTTAAGAATTAACTTAAATTACTAACTAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt < 2:234870/141‑1 (MQ=255) aattttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt < 1:41690/139‑1 (MQ=255) aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt < 1:41643/141‑1 (MQ=255) gtaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg > 1:225301‑M1/110‑141 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata < 1:119524‑M1/36‑1 (MQ=255) aCTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa > 2:65718/1‑141 (MQ=255) tgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaat < 2:240686/141‑1 (MQ=255) gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt < 1:38002/141‑1 (MQ=255) gcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc > 1:169395‑M1/98‑141 (MQ=255) gcgtccattgtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc > 1:169397‑M1/98‑141 (MQ=255) cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatacaatggataagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt > 1:90039‑M1/97‑141 (MQ=255) gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTTAACTTAATAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCtttt > 2:97261/1‑139 (MQ=255) ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATCAAATGTATATTTAATTCTTGTCGGt > 1:136293‑M1/90‑141 (MQ=255) agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg > 1:12042‑M1/87‑141 (MQ=255) aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGGTAAGa > 1:166811/1‑141 (MQ=255) aGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaa > 2:55393/1‑141 (MQ=255) ccTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa < 1:240660/141‑1 (MQ=255) gTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTACGAAAAATGa > 2:22048/1‑141 (MQ=255) tacgttacgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa > 2:241790‑M1/78‑141 (MQ=255) cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGtcggtaagaaacaaaaatg < 2:133167‑M1/66‑9 (MQ=255) ccgtccgctagagtagaacgttgacaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa > 2:2385‑M1/76‑141 (MQ=255) actatttataTGTTGATGCAAACAAAACAATTTAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt < 2:254066/139‑1 (MQ=255) tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg < 1:119641‑M1/69‑1 (MQ=255) ctagagtagagcgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa > 1:25376‑M1/69‑141 (MQ=255) tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGGAAAc > 2:187266‑M1/68‑141 (MQ=255) ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg > 2:183011/1‑141 (MQ=255) caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt < 1:55393/141‑1 (MQ=255) ggcagttatttaaacaaattttggtcaaaaaataaaaaaaaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa > 2:248256‑M1/50‑141 (MQ=255) gttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 2:49935‑M1/46‑136 (MQ=255) gttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa < 1:49935‑M1/91‑1 (MQ=255) tttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa < 2:46556‑M1/98‑1 (MQ=255) ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACAACTGAATGACAATATGTCAACGTTAAt < 1:36618‑M1/99‑1 (MQ=255) aaactTAAGAAATAACTTGAATTACTAACGCAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt < 2:39577/140‑1 (MQ=255) aaGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGGTATTGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTAAATAAAAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa < 2:78414/141‑1 (MQ=37) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa > 2:234347‑M1/36‑141 (MQ=255) aaattttggttgaaaaagaaaatggacaagataaAACAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTTTTTTCGTTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa < 2:193162‑M1/107‑1 (MQ=255) aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc < 1:38482/141‑1 (MQ=255) aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata > 2:33779‑M1/24‑141 (MQ=255) aTTACTAACAAAAATTAATTTAAAAAAGTTTTTGACGTAAATGTTAATAAAATGTACAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa < 2:220145/141‑1 (MQ=37) taatAACCAAAATTAATTTTAAAAAGTTATTGACTTAAAGGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt < 2:125056/138‑1 (MQ=255) ataaaatgaataagataaAAAAAGTTATTGACTAAAATGTTGATAAAATGTATGAGTAATTCTTGTCTGTAAGAAAAATGACCATTGAAAACTGAATTACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa > 2:148805‑M1/19‑141 (MQ=255) aaaaaatgagacagataaAAAAAGTTATTGATTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCATAAAACGTAACTATAAGTTa > 2:207742‑M1/19‑141 (MQ=255) tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc < 1:261041/141‑1 (MQ=255) aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGCATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa > 1:185732‑M1/17‑141 (MQ=255) aaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa > 1:14307‑M1/17‑141 (MQ=255) aaaaatgacaagataaAAAAAGTTATTGACGTAAATGTTGATAAAATATAGAATTAACTCATGACGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATTTCAACGTTAAGACCAAAAAACGTAACTATAAGTTACa > 2:154149‑M1/17‑141 (MQ=255) aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa > 1:230104‑M1/15‑141 (MQ=255) aatggaaaagacaaAAAAAGGTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGGCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa > 1:250812‑M1/15‑141 (MQ=255) taattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattact > 2:145334/1‑139 (MQ=255) aagagaaaaaaaaTTATTGAATTAAATGTTGAAAAAAGGAAAAATTAATTCTTGCAGGAAAGAAAAATGAACATTGAAAATTGAATGACAATATGTAAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACattatt < 2:110712/134‑1 (MQ=12) tttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt < 1:183011/141‑1 (MQ=255) aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAAATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAAAACGTAACTATAAGTTACAAACATAATTTAGTAt < 1:190994‑M1/140‑1 (MQ=255) aaaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt > 2:8687/1‑141 (MQ=255) gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa > 1:12001/1‑141 (MQ=255) gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa > 1:12044/1‑141 (MQ=255) gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa > 1:54612/1‑141 (MQ=255) ttGACTTAAATGTTGATAAAAGGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAGTGAAAACTGAATGACGATATGTCAACGTTAATTCCAAAGAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa < 1:188819/141‑1 (MQ=21) ttAATAAAATGTATAATTAATTCTTGACGCTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat < 2:193885/141‑1 (MQ=17) | AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT > CP000730/512503‑512781 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |