Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 512,643 | 0 | A | G | 41.5% | 47.9 / 56.0 | 65 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (22/16); new base G (12/15); total (34/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT > USA300TCH1516_ALE/512503‑512781 | aGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta > 1:225213/1‑141 (MQ=255) aGGTGTAATTTTTACTGTAGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta > 1:146577/1‑141 (MQ=255) ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa < 1:104781/141‑1 (MQ=255) agttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa < 1:292610‑M1/30‑1 (MQ=255) aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatcaaatggacaagatacCAAAATTTATTGACTTTCATGTTGATAAAATGta < 2:242845‑M1/34‑1 (MQ=255) aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGta < 2:188760‑M1/34‑1 (MQ=255) ttATTGGTGTGAAGAAGTATAAAAACATGTAAAGTATTTATTTGTTGAGGCAAACAAAACAACTCAAATTAAGAAATAAACTGTTATATTTTTTTAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat < 2:51326/141‑1 (MQ=255) ttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat < 1:75594/141‑1 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata < 2:34661‑M1/36‑1 (MQ=255) aCTGTTGTTAAGCAGTGTGAAAGCCTGTATAGTATTTATTTGTTGAGGAAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa < 2:246691/141‑1 (MQ=255) ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggccgctgctttatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt > 1:221452‑M1/103‑141 (MQ=255) cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatacaatgggcaagataaAAAAAGTTATTGACTTAAATGTTGATAGAATGTAtgatta < 1:11098‑M1/40‑1 (MQ=255) gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt < 1:63396‑M1/46‑1 (MQ=255) gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTgtt < 2:4371‑M1/49‑2 (MQ=255) aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg > 1:144719/1‑141 (MQ=255) ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa > 2:174371/1‑141 (MQ=255) ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa > 1:14517/1‑141 (MQ=255) tagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGCCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTaa < 2:86820‑M1/54‑1 (MQ=255) aGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGACTTTGCGTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaa < 1:70070/141‑1 (MQ=255) tataGTATTTATTTGTTGAGGCAAACAACACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa > 1:64705/1‑141 (MQ=255) acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc < 1:37193‑M1/65‑1 (MQ=255) ttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt < 1:33831‑M1/68‑1 (MQ=255) gTATTTATTTGTTGAGACAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt < 2:188277/141‑1 (MQ=255) atttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATGGACTTAAATGTTAATAAAATGTATAATTAATTTTTGTCGGTAAGAAAAATGAACATTGa > 2:88637/1‑141 (MQ=255) ttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa < 1:174038/141‑1 (MQ=255) tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa > 2:38584/1‑141 (MQ=255) atgtgttgAGGCAAACAAAACAACTCAACTTAAAGAATAACTTGAAGTAGTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCGTGTCGGTAAGAAAAATGAACATTGAAAAc < 1:10784/138‑1 (MQ=255) agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa > 1:34661‑M1/65‑141 (MQ=255) gttgAGGCAAACAAAACAACTCAACTTAAGAAATCACTTGAATTACTAACGAACAATACTAGTAAAAAGTTATTGACTTAAATGTAAATAAAATGTATAATTAATTCTTCTCGGTAAGAAAAATAAACATCGAAAACTGaa > 2:92871/1‑141 (MQ=255) agaacgttgacaggaagttttttaatcaaatgttcgttaaaaaataaaatagagaggataaGAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg > 2:261101‑M1/62‑141 (MQ=255) gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAAGAAAATGAACATTGAAAACTGCATAacagt < 2:30927‑M1/85‑3 (MQ=255) caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAat < 1:65766/141‑1 (MQ=255) aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat > 1:9706/1‑141 (MQ=255) aaaagaagttttataataaaattatggttaaaaaatattttgaacaaaataaAAAAAGTTATTGACTTAAATGTAGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt > 2:214114‑M1/53‑141 (MQ=255) aacaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa < 2:225213/141‑1 (MQ=255) attAAGATATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc < 2:100278/140‑1 (MQ=255) tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt > 1:160779‑M1/31‑141 (MQ=255) aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa > 1:155673/1‑141 (MQ=255) cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc > 2:225602/1‑141 (MQ=255) aCTAACGAAAATTAATTTTAAAAAGTTATTGACATAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt > 1:122664/1‑141 (MQ=255) cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGCACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAATtatacatta > 2:81276/1‑136 (MQ=255) cTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCTAAAAAACGTAACTATAAGTTa > 2:88751/1‑141 (MQ=255) aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa < 2:26786/141‑1 (MQ=255) cGAAAATTAATTTTAAAAAGTTATGGACTTAAATGTTAATAAAATGTATAATTAATTCTTGACGATAATAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa < 1:142989/141‑1 (MQ=255) tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa > 1:120961‑M1/13‑141 (MQ=255) attaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta < 1:174371/141‑1 (MQ=255) ttaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACATTAATTCCAAAAAACGTAACTATAAGTTACAAACattat > 2:169467/1‑141 (MQ=255) agataaaaaaaaATATTAAATTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAAGAACATTAAAAACTGAATGACAATATGTAAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTAttt < 2:121528/135‑1 (MQ=19) attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa > 1:276567/1‑141 (MQ=255) attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa > 1:208774/1‑141 (MQ=255) attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACATAACTATAAGTTACAAACCTTATTTa > 2:83530/1‑141 (MQ=255) atatAAAAAGTTATTTACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg < 1:250175/138‑1 (MQ=255) tataaaaGTGTTTGACTTAAATGTTAATAAAAGTTATAATAAATACTTGTCGATAAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa < 2:254237/138‑1 (MQ=18) tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt < 2:134198/141‑1 (MQ=255) ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATAAg > 2:82803/1‑141 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa < 1:206543/141‑1 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa < 1:252949/141‑1 (MQ=255) tGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTAACAAACATTAGTTAGAAATTATGAGCTaa > 2:242809/1‑141 (MQ=255) gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt > 1:200800/1‑141 (MQ=255) aacgtAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt < 1:82027/137‑1 (MQ=255) tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGCATAACAATATGTCAACGTTAATTACAAAAAACGTAACTATAACTTCCAACCATTATTTAATATATATGAGCTAATGaaa > 2:93147/1‑141 (MQ=18) tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaaa > 2:20580/1‑141 (MQ=255) aaaTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATAACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAGTTATAAACTAATCAAAc > 2:201698/1‑141 (MQ=18) aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat > 2:86254/1‑141 (MQ=255) aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat > 1:292361/1‑141 (MQ=255) aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat > 2:158124/1‑141 (MQ=255) ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat > 2:234590/1‑141 (MQ=255) ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat > 1:99969/1‑141 (MQ=255) | AGGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT > USA300TCH1516_ALE/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: ‑ |