Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 512,643 | 0 | A | G | 39.7% | 45.4 / 64.4 | 63 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (20/18); new base G (11/14); total (31/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA > USA300TCH1516_ALE/512501‑512782 | aacacagaagttaagctccttagcgtcgatggtagtcgaacttacgcaccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa < 1:278577‑M1/29‑1 (MQ=255) cacagaagttaagctccgtagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAt > 1:179315‑M1/118‑148 (MQ=255) atgtgtAATTTTTACTGTCTTTAAGCAGTTTGTAATCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTGTAAAAAGTTATTGACTTAAATGTTAATAAAAt < 1:403798/146‑1 (MQ=255) acagaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaagatcaaatggccaagttaaAAAAAGTTATTGACTTAAATGTTGATAAAATg < 1:84259‑M1/32‑1 (MQ=255) gtgtAATTTTTAATGTTGTTGAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAACAATAAATTTAAAAAGTTATTGACTTAAGTGTTAATAAAATGt > 1:447107/1‑148 (MQ=255) cagaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt < 1:589640‑M1/33‑1 (MQ=255) aaTTTTTACTGTTGTGAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa < 1:202804/148‑1 (MQ=255) ttttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta < 1:478804/148‑1 (MQ=255) ttttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta < 1:82742/148‑1 (MQ=255) tccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTATATGTTGATAAAATGTATAATTAATTCt > 1:335916‑M1/104‑148 (MQ=255) tccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt > 1:519838‑M1/104‑148 (MQ=255) cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt < 1:621442‑M1/52‑1 (MQ=255) aGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAGAACAACTCAACTTAAGAAATCACTTGAATTACTAACGAAAATTAATTTTAAAAAGATATTGACTTAGATGTTAATAAAATGTATAATTAATTGTTTTCGGTaa > 1:181158/1‑148 (MQ=255) tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa < 1:482047/148‑1 (MQ=255) ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACGTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa > 1:175610/1‑148 (MQ=255) gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTGATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa < 1:163510‑M1/60‑1 (MQ=255) tcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa < 1:267858‑M1/64‑1 (MQ=255) tcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa > 1:135039‑M1/85‑148 (MQ=255) acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa > 1:512050‑M1/77‑148 (MQ=255) atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa > 1:400008/1‑148 (MQ=255) tttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa > 1:530189/1‑148 (MQ=255) atttGTTGAGGCAAACAAAACAACTCAACTTAAGACATAACTTGAATTACTAAGGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAGGAGAAATGAACATTGAAAACTGAATGa < 1:354200/148‑1 (MQ=255) atttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa < 1:437347/148‑1 (MQ=255) tttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATgcc > 1:164669/1‑146 (MQ=255) ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa < 1:618264/148‑1 (MQ=255) tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa < 1:247974/148‑1 (MQ=255) agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa > 1:404400‑M1/65‑148 (MQ=255) ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata > 1:57951/1‑148 (MQ=255) gCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACCATATGtaa > 1:187467/1‑146 (MQ=255) gCAAACAAAACAACTCAACTTAAGAAAAAACTTGAATTACTAACGAGAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 1:354149/1‑148 (MQ=255) aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg > 1:106397/1‑148 (MQ=255) acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt < 1:617864/148‑1 (MQ=255) caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGAAAGAAGAAGCAACAGTGAAAACTGAATGACAATATGTAAAAGTTa < 1:488552‑M1/97‑1 (MQ=255) ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTACGCCAACGGAACATGGAAAACTGAATGACACTATGTCAacgttcat < 1:15693‑M1/99‑4 (MQ=255) ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt < 1:255012‑M1/99‑1 (MQ=255) tttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATTTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACAACTGAATGACAATAAGTCAACGTTAATTACaaa < 1:28459‑M1/105‑1 (MQ=255) tttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGGGTTGAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa < 1:300849‑M1/105‑1 (MQ=255) tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa > 1:308053/1‑148 (MQ=255) tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa > 1:43591‑M1/42‑148 (MQ=255) aaGAAATAACTTGATTTACTAACGAAAATTAATTTTGAAAAGATATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAGAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATCCCAAAAAACGTa < 1:446432/148‑1 (MQ=38) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTTAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa > 1:628115‑M1/36‑148 (MQ=255) aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa < 1:203794/148‑1 (MQ=255) gAAAGAATTGGAATTACTAACGAAAATTAATTTTAAAAAGTTATCGACTTAAATGTTAATAAAATGTATAATTAACTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:356040/148‑1 (MQ=37) aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACta > 1:177697/1‑148 (MQ=255) ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt < 1:474361/148‑1 (MQ=255) ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTTAATTCCAAAAAACGTAACTATAAGtt > 1:427623/1‑148 (MQ=255) aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAGGTTACaa > 1:493444/1‑148 (MQ=255) tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc < 1:416644/148‑1 (MQ=255) aataaaatggacagaagaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAGATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa > 1:110901‑M1/20‑148 (MQ=255) cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACat > 1:322342/1‑148 (MQ=255) gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt > 1:28627/1‑148 (MQ=255) aattaattTTAAAAAGTTATTGACTTAAATGCTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGCACATTGAAAACTGAATGACAATATGTCAACGTTAATTATAAAAAACGTAACTATAAGTTCCAAACAATAATTaat > 1:214684/1‑146 (MQ=25) ttttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTGTAAGTTACAAACATTATTTAGTATTTAt > 1:229557/1‑148 (MQ=255) aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAGTTAGTATTTATGa > 1:27807‑M1/3‑148 (MQ=255) aaaGTTATTGACTTAAATGTTGATAAAATGTATGATTAATTCTTGTGGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATCAGGTACCAACATTATTTAGTATTTATGAGcag > 1:111761/1‑146 (MQ=21) aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa < 1:299262/148‑1 (MQ=255) aaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa > 1:580482/1‑148 (MQ=255) gTTATTGACATAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATc < 1:548055/148‑1 (MQ=255) tATTGACTTAAATGTTAAGAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa < 1:16044/148‑1 (MQ=255) aCTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTTCAATCATTATTTAGTATTTCTGAGCTAATCAAAcatc < 1:465830/148‑1 (MQ=255) cTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca < 1:67646/148‑1 (MQ=255) ttAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat < 1:632551/148‑1 (MQ=255) tAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTAATGAGCTAATCAAAcatcaca > 1:71354/1‑146 (MQ=255) | AAAGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA > USA300TCH1516_ALE/512501‑512782 |
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: ‑ |