| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 512,643 | 0 | A | G | 70.3% | 44.9 / 25.9 | 37 | intergenic (+711/‑150) | recR/rrs1 | recombination protein RecR/16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (7/4); new base G (11/15); total (18/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAA > CP000730/512508‑512775 | aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa < 2:100718‑M1/29‑1 (MQ=255) ttaagctccttagcgtaaatggtagtcgaaattacgttaagctagagtagaacgttgccaggcagttttttaataaaattttggttaaaaaataaaatgaaaaaaataaAAAAAGTTATTGACTTAAATGTTGATAAAATg > 2:234058‑M1/110‑141 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata < 1:247214‑M1/36‑1 (MQ=255) aCTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa > 2:180219/1‑141 (MQ=255) gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta > 2:225629/1‑141 (MQ=255) gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGATATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa < 2:113151‑M1/56‑1 (MQ=255) aaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa > 2:42676/1‑141 (MQ=255) cttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑Tg > 2:217192‑M1/80‑141 (MQ=255) ttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGa < 1:279693‑M1/63‑1 (MQ=255) gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATGACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACAtt > 2:56987/1‑141 (MQ=255) gagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTg < 2:111537‑M1/76‑1 (MQ=255) gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTGATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGa < 2:209622‑M1/81‑1 (MQ=255) gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAATGAACATTGAAAACTGAATGACAatat < 1:92651‑M1/88‑1 (MQ=255) aggcagatttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGATGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCa < 2:248219‑M1/91‑1 (MQ=255) acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAAc < 1:131047/141‑1 (MQ=255) tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTaa > 2:100905/1‑141 (MQ=255) tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt > 1:100718‑M1/42‑141 (MQ=255) tAAGAAATAACTTGAATTACTAACGAAAATTAATTATAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa < 2:253982/141‑1 (MQ=255) ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa > 1:173587‑M1/32‑90 (MQ=255) ttttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa < 2:173587‑M1/59‑1 (MQ=255) aaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTACAAAAAACGTAACtata > 2:146849‑M1/24‑141 (MQ=255) aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa < 1:180219/141‑1 (MQ=255) aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa > 2:212569/1‑141 (MQ=255) tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa < 2:4552‑M1/129‑1 (MQ=255) ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat < 2:240130‑M1/130‑1 (MQ=255) ggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat < 2:226176‑M1/130‑1 (MQ=255) gacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt < 2:267722‑M1/131‑1 (MQ=255) gacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt > 1:226176‑M1/11‑141 (MQ=255) aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt < 1:100905/141‑1 (MQ=255) ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg > 2:279693‑M1/5‑141 (MQ=255) taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt > 2:108658‑M1/4‑141 (MQ=255) aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt > 2:173644/1‑141 (MQ=255) aaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt < 1:234115/141‑1 (MQ=255) aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt > 2:93405/1‑141 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt > 1:156296/1‑139 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt < 2:156296/139‑1 (MQ=255) tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaaa > 2:144749/1‑141 (MQ=255) | TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAA‑TGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAA > CP000730/512508‑512775 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |