Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 512,643 | A→G | 64.2% | 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 | 64.2% | 37.7 / 30.2 | 41 | intergenic (+711/‑150) | recR/rrs1 | recombination protein RecR/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (10/5); new base G (11/15); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.97e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 |
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA > CP000730/512503‑512773 | aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta < 1:137180/141‑1 (MQ=255) agaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgata < 2:138410‑M1/27‑1 (MQ=255) ttattttctACTTTTGTCAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGACGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAACATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa < 1:214592/137‑1 (MQ=255) gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta < 2:234767/141‑1 (MQ=255) cttagcgtcgatggtaatcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta < 1:169177‑M1/40‑1 (MQ=255) tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatccaattatggctaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTg < 2:95026‑M1/47‑1 (MQ=255) aaaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg < 2:108277/140‑1 (MQ=255) gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATGAAATGTATAATTAATTATTGTCGGTAAGa < 1:101551‑M1/56‑1 (MQ=255) gaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa > 2:248663‑M1/83‑141 (MQ=255) acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc < 1:174690‑M1/65‑1 (MQ=255) acgttccgctagagtagaacgtgaccatacactgtttccatcaaattttggttaaaaaataaaatggacaagataaAAAAAATTAATGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATgaaa < 1:67911‑M1/65‑2 (MQ=255) cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggataaagaaaaaaatgggcaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTAGGTAAGAAGAATGAACa < 2:156523‑M1/66‑1 (MQ=255) tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg > 2:235689‑M1/73‑141 (MQ=255) tgttgAAGCAAACAAAACAACTAAAATTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCAAGTAAGTAAAAAAAATGAACATTGAAAACTGa < 2:251395/141‑1 (MQ=255) gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa > 2:27283/1‑141 (MQ=255) gttgAGGAACAGAAAGCAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGac > 2:275423/1‑140 (MQ=255) ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt > 1:62452/1‑141 (MQ=255) tcacgCAAACAAAACAACTCAACTTAAGTAATCACTCGAATTACTCACAAAGATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAACATATATAATTAATTTTTGTCTGTAAGAAAGAGGCACATTAaaaaatcgatt > 2:123998/5‑134 (MQ=255) gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAACCTGAATGa < 1:167270‑M1/81‑1 (MQ=255) gttgccaggcagttttttaatcaaattttgggtaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGGCGGTAAGAAAAAGGAACATTGaaaactgaagtgcaat < 1:142383‑M1/85‑8 (MQ=255) gtagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc > 2:112949‑M1/49‑141 (MQ=255) acaacTCGACTTAAGAAATAACTTGAATTACTAACGCAAATTAAATTTAAAAAGTTATTGACTTACAGGTTAATAACATGTATAATTAATTCTCGTCGGTAAGAAAAAGGAACATTGCACACTCAATGACAACATGGCAAc > 2:49645/1‑141 (MQ=255) acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaag > 2:273691/1‑140 (MQ=255) ttttttaaaacaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt > 2:275078‑M1/46‑141 (MQ=255) caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt < 1:186232/141‑1 (MQ=255) taatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc < 1:116316‑M1/102‑1 (MQ=255) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAGTTTGTACATTGAAAACTGAATGACAATATGTCGACGTTAATTCCaaaa < 1:21814‑M1/106‑1 (MQ=255) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAGAATATGTGAACGTTCATtccacaa < 1:4753‑M1/106‑4 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:85497‑M1/114‑1 (MQ=255) aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCAATAGAATGTATAATTAATTTTTGTTGGTAATAAAAATGAACATTGAAATCTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACtata > 2:280515/1‑141 (MQ=37) aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa > 1:138410‑M1/17‑141 (MQ=255) ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGAAAAATCCAAAAAACGTAACTATAAGTTACAAACat > 1:247678‑M1/12‑141 (MQ=255) acaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACatta > 1:197078‑M1/10‑141 (MQ=255) acaagataaAAAAAGTAATTGACATAAATGTTGATAAAATGAATAATTAATTCTTGTCAGTAAGAAAAATGAACATTGAAAACTTAATGACCATATGTCAAAAAAAAAAACAAAAAACGTAACTATAAGTTACAAACatta > 2:231731‑M1/10‑141 (MQ=255) taattTTAAAAAGTTAGTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt < 2:84251/141‑1 (MQ=255) aagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattatt < 2:85290‑M1/134‑1 (MQ=255) attTTAAAATGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa > 1:208124/1‑141 (MQ=255) ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg > 2:236400‑M1/5‑141 (MQ=255) aaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt < 2:61303/141‑1 (MQ=255) ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg < 1:96879/141‑1 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATAAACTaa > 2:191636/1‑141 (MQ=255) gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt > 1:268879/1‑141 (MQ=255) gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt > 2:220583/1‑141 (MQ=255) cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa > 2:148383/1‑141 (MQ=255) | AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA > CP000730/512503‑512773 |
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: ‑ |