Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 512,643 | 0 | A | G | 48.7% | 7.7 / 41.2 | 39 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (14/6); new base G (7/12); total (21/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.62e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA > USA300TCH1516_ALE/512506‑512783 | agaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTCTTGACTTACATCTttata < 2:25062‑M1/27‑1 (MQ=255) ttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg < 2:220694‑M1/32‑1 (MQ=255) tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtata > 1:108747/1‑141 (MQ=255) cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta < 1:83162‑M1/40‑1 (MQ=255) tgttAAGCAATTTGAAAGCCTGTATAGTAATAATTTGTTGAGACAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaat < 2:222490/141‑1 (MQ=255) gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATATTAATAAAATGTATAATTAATTATTGt > 2:145976/1‑141 (MQ=255) atggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCg < 1:1772‑M1/50‑1 (MQ=255) gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg < 1:13485/141‑1 (MQ=255) gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg < 2:56376/141‑1 (MQ=255) gtcgaacttacgttccgctagggtagaacgttgccaggcagttttttaattaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTTTTGTAGGTAAGa < 1:185982‑M1/56‑1 (MQ=255) gtcgaacttaagttcagctagagtagaacgttgacaggcagttttttaataaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa > 2:221373‑M1/86‑141 (MQ=255) aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa > 2:188578/1‑141 (MQ=255) gaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa > 1:235913‑M1/83‑141 (MQ=255) acttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTACGAAAAAt < 2:102490‑M1/61‑1 (MQ=255) tttataGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg < 2:36407/139‑1 (MQ=255) tacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa > 1:220694‑M1/78‑141 (MQ=255) acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagaaaaAGTAAGTGATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc < 2:84344‑M1/65‑1 (MQ=255) gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt > 1:2071/1‑141 (MQ=255) gctagagtagagcgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTATGAAAAATGAACATTTaaa < 2:194394‑M1/72‑1 (MQ=255) tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa > 1:131182/1‑141 (MQ=255) tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa > 1:186682/1‑141 (MQ=255) agtagaacgttgccaggcagttttttaatcaaatttcggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa < 1:42286‑M1/77‑1 (MQ=255) gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaactggacaagataaAAAAAGTTATTGACTTAAATGTTGATACAATGTATAATTAATTCTTGACGGTCAGAAAACTGAACATTCAAAACTGCAtgc < 2:149120‑M1/81‑2 (MQ=255) taaaacaaaacaACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa < 2:232770/139‑1 (MQ=255) acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGATATTGACTTAAATGTTAATAAAATGTATAATTAATTATTATCGGTAAGAAAAATAAACATTAAAAACTGCATAACaaaacg > 2:112310/1‑137 (MQ=255) aaggcagtattttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc > 2:201607‑M1/52‑141 (MQ=255) aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 2:92458/1‑141 (MQ=255) tttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTGGCTAAAAGGTATAATTAATTCTTGTCGGTAATAAACATGACCATTGAAAACTGAATGACAATATGTCAACGTTaa < 1:32988‑M1/98‑1 (MQ=255) actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGCCGCCAGGACAAAGGAACGTTGAAAACTGAATGACAATATGTCAACGTTAATTc > 1:222809/1‑141 (MQ=38) ttAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa < 2:219097/141‑1 (MQ=255) tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTAAACGTTAATTCCaa > 1:208080/1‑141 (MQ=255) tcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTAATTTATTAAATGTTGATAAATTGTCAGAGTAATTCTTGTCGGTAAGAAAATTGTAGATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaa > 1:210671‑M1/37‑141 (MQ=255) tAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATGAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt > 2:34189/1‑141 (MQ=255) tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTGCAAAAAACGTAACTATAAGTTAc > 2:98731/1‑141 (MQ=255) aaaatggacaagataaAAAAAGTTATTGACTTAAAGGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGACCATTTAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa > 2:96651‑M1/17‑141 (MQ=255) caagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGTAGAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattat > 2:189961‑M1/9‑141 (MQ=255) attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa > 1:85300/1‑141 (MQ=255) ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg < 1:210065/141‑1 (MQ=255) atAAATGTTGATAAAATGAATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaa < 2:117702/140‑1 (MQ=255) ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACACTATTTAGTATTTATGAGCTAATCAAAcatcat > 2:222137/1‑141 (MQ=255) aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATaa > 1:230123/1‑141 (MQ=255) | TGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATAA > USA300TCH1516_ALE/512506‑512783 |
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: ‑ |