Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 512,643 | 0 | A | G | 69.6% | 22.7 / 16.5 | 23 | intergenic (+711/‑150) | recR/rrs1 | recombination protein RecR/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (4/3); new base G (12/4); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA > CP000730/512494‑512783 | aGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATATTTTTTATTTGTTAAGTTTAACAAATTAATTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa < 1:42269/151‑1 (MQ=255) cacagaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGta < 2:252941‑M1/34‑1 (MQ=255) cagaagttaagctccttagcgtcgatggtagtcgaacttacgttgcgcttgagtagaacgttgccaggcagttttttagtcaaattgtgggtaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTAGATAACATGtata < 2:452861‑M1/36‑1 (MQ=255) aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTc > 1:116530‑M1/108‑151 (MQ=255) ttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCgg > 2:353046‑M1/101‑151 (MQ=255) gttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGAAAACAAAACAACTCAACTTAAGAAATAACTTGACTTACTAAAGAACATTATTTTTAAAAAATTATTGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa > 1:294221/1‑151 (MQ=255) tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg < 2:111289/151‑1 (MQ=255) tacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc > 2:158682‑M1/78‑151 (MQ=255) ccgctagagtagaacgttgccaggcagtattttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTCTTGACGTAAATGTTGATAAAATGTATAATTAATTCTAGTCGGTAAGAAAAATGAACATTTAAAACTGAATg < 2:279499‑M1/80‑1 (MQ=255) caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt > 1:308347/1‑151 (MQ=255) caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAATTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt > 1:385071‑M1/52‑151 (MQ=255) ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt < 1:186666‑M1/111‑1 (MQ=255) gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTAAAAAACAAAAAACGTAACtat < 2:308347/151‑1 (MQ=255) aaaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata < 1:162470/151‑1 (MQ=255) taaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa > 1:314355‑M1/25‑151 (MQ=255) caagataaAAAAAGTTATTGGCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAATATTGAAAACTGAATGATAATATGTCAACGTTAATTACAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa > 2:195245‑M1/9‑151 (MQ=255) attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa > 2:102924/1‑151 (MQ=255) taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc > 1:280564‑M1/4‑151 (MQ=255) tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAca > 2:106286/1‑151 (MQ=255) naTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAca > 2:434143/2‑151 (MQ=255) gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTCTGAGCTAATCAAAcatcat > 2:224224/1‑151 (MQ=255) gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat > 2:332257/1‑151 (MQ=255) aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa > 2:301664/1‑151 (MQ=255) cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATaa > 2:306256/1‑151 (MQ=255) | AGATGACAAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA > CP000730/512494‑512783 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |