Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 512,643 | 0 | A | G | 51.7% | ‑6.2 / 81.3 | 60 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (18/11); new base G (14/17); total (32/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA > USA300TCH1516_ALE/512507‑512782 | gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa > 1:11587/1‑141 (MQ=255) agttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaaccttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa > 1:199119‑M1/112‑141 (MQ=255) aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACGTGAATTACGAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaaaa > 1:65930/1‑140 (MQ=255) aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCAATAAAGGATGTAGTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt > 1:328972/1‑141 (MQ=255) ttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATg > 1:11639‑M1/110‑141 (MQ=255) cgccttagcgtcgatggtagtcgaacttacataccgctagagtagaaagttgacaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaa > 2:344098‑M1/105‑141 (MQ=255) cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCt < 1:299744‑M1/45‑1 (MQ=255) agagggtaggcgaacttacgtaccgctagagtagaacattgaatggctgtattgtaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGt > 2:38785‑M1/94‑141 (MQ=255) gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTc < 1:332985‑M1/49‑1 (MQ=255) gTTTGACAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTctg > 2:277467/1‑139 (MQ=255) ttagacgaaaatactttcctatagagtaaaactttgaaatgtagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa > 2:315845‑M1/89‑141 (MQ=255) ttaaaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa < 2:322073/138‑1 (MQ=255) ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTa > 2:226414/1‑141 (MQ=255) agtagaacttacgttacgagagagtagaaagttgacaggcagttttttaataaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATCATTAATTCTTGTCGGTAAg < 2:218606‑M1/55‑1 (MQ=255) cgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaa > 1:170918‑M1/84‑141 (MQ=255) ccTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa > 1:301687/1‑141 (MQ=255) taGTATTGATTTATTGAAGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa < 2:270664/141‑1 (MQ=255) tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa > 1:8358/1‑141 (MQ=255) tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa < 2:11587/141‑1 (MQ=255) atagagtagaaggtagccaggcagttttttaataaaaatttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATGATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa > 2:286304‑M1/69‑141 (MQ=255) ttgttgAGGTAAACAAAACAACTCAACTTAAGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg < 2:61671/141‑1 (MQ=255) agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatcaaatggacaacatacAAAAAGTTCTTGCCTTATATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATTAACATTGAAAACTGa < 2:323448‑M1/77‑1 (MQ=255) agtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacacgataaAAAAAGTTATTGACTTACATGTTGATATAATGTATAATTAATTCTTGTCGGTAAGACAACTGAACATTGAAAactca < 2:287213‑M1/77‑3 (MQ=255) gtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa < 1:183311‑M1/78‑1 (MQ=255) ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt > 1:278056/1‑141 (MQ=255) ttgAGGCAAAAAAAACAACTAAACTTAAGAAATAACTTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt < 2:65930/141‑1 (MQ=255) gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTTAACACTGAATGa < 2:98150‑M1/81‑1 (MQ=255) aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc > 1:66530/1‑141 (MQ=255) gttgccaggcagttttttaatccaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACAACTGAATGACAat < 1:329714‑M1/85‑1 (MQ=255) gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATCATCATCTAAAACTGTATGAcaac < 2:243039‑M1/85‑2 (MQ=255) ttgccaggcagttttttaatcaaattttggttaaaaaataaaatggagaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata < 2:178008‑M1/86‑1 (MQ=255) aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat > 1:140692/1‑141 (MQ=255) gccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg < 1:170401‑M1/88‑1 (MQ=255) acaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATCATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTCAATGCCACTATg > 2:189359/1‑141 (MQ=255) caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc > 1:41107‑M1/52‑141 (MQ=255) aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 1:43996‑M1/51‑141 (MQ=255) aggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa < 1:335765‑M1/91‑1 (MQ=255) aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 1:320388/1‑141 (MQ=255) aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGATAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 2:255370/1‑141 (MQ=255) aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTAATGGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATAACAATATGACAAcat > 2:41402/1‑139 (MQ=255) actcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt < 1:138514/141‑1 (MQ=255) aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGATCACCAGTATGTCAACGTTAAtt > 1:216236/1‑141 (MQ=255) aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAAAACTGAATGACAATATGTAAACGTTAAtt > 2:198331/1‑141 (MQ=255) actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATAACAATATGTCAACATTAATTc > 2:272966/1‑141 (MQ=255) atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa > 1:284281‑M1/38‑141 (MQ=255) gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATAGAATACTTAATGACAATATGTCAACGTTAATTCCaaaaa > 1:64583/1‑141 (MQ=255) aaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTTACTTATATTTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATCAACATTGAAAACTGAATGACAATATGTCAACGTTACTTCCaaaaa < 2:89140‑M1/107‑1 (MQ=255) ttggtgaaaaaataaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAACATGTCAACGTTAATTCCAAAAAACGTa > 2:231337‑M1/30‑141 (MQ=255) tggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACCATCTGTCAACTGTAATTCCAAAAAACGTaa < 2:138643‑M1/113‑1 (MQ=255) aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata < 1:226414/141‑1 (MQ=255) aattaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAAAACGTAACTATAAGTTATAAACatt > 2:212587/1‑141 (MQ=255) ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTTCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa < 2:313437/141‑1 (MQ=255) aaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa < 2:320388/141‑1 (MQ=255) aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg > 1:275170/1‑141 (MQ=255) gTTATTCACACAAATGGTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa < 1:129202/141‑1 (MQ=18) ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg < 1:178008/141‑1 (MQ=255) ttATTGAATTAAAGGTTAAACGAAGGTATAATTAAGTCTTGTCGGGAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg < 2:297575/141‑1 (MQ=17) ttGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa < 2:278056/141‑1 (MQ=255) tGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTaa > 1:249486/1‑141 (MQ=255) gTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca > 1:301611/1‑141 (MQ=255) gTTGATAAAATGTATAATTAAGTCTTGTCGGTAAGAAAAATGAACATTGAAAGCTGGATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatca < 2:190399/141‑1 (MQ=255) taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa < 1:260605/141‑1 (MQ=255) taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa < 2:8358/141‑1 (MQ=255) | GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA > USA300TCH1516_ALE/512507‑512782 |
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: ‑ |