Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 277,427 | 0 | G | C | 43.5% | 6.6 / 30.4 | 23 | A84A (GCG→GCC) | USA300HOU_RS01215 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (9/4); new base C (3/7); total (12/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.95e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGACATATC > USA300TCH1516_ALE/277296‑277561 | tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat > 2:242143/1‑141 (MQ=255) tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat > 1:395218/1‑141 (MQ=255) tCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAatgat > 1:287036/1‑141 (MQ=255) tCCGACATACATAATGAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAGATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAatgat > 1:103630/1‑141 (MQ=255) tCCGACATACATAATGAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAatgat > 1:318830/1‑141 (MQ=255) cataATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTa > 2:333137/1‑141 (MQ=255) gAATCTATCAATAAACAAATGGAAGCCAAAGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCtt < 2:41229/141‑1 (MQ=255) aTCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTaa > 1:33555/1‑141 (MQ=255) aTCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTaa > 2:184759/1‑141 (MQ=255) aCAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGatat < 2:137276/141‑1 (MQ=255) caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt < 2:144451/140‑1 (MQ=255) caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt < 1:174003/140‑1 (MQ=255) caaaGCGCATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGt < 1:173985/140‑1 (MQ=255) gcgcATGAAACAGCGCAAAATGCAGATTTAAAAACCGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTAttt > 1:115885/1‑141 (MQ=255) gATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCa < 1:242143/141‑1 (MQ=255) gATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCa < 2:395218/141‑1 (MQ=255) gCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACg > 1:37373/1‑141 (MQ=255) gCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACg > 2:39228/1‑141 (MQ=255) tcaatcGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAAtt > 1:346582/1‑141 (MQ=255) tcGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTtgtg < 1:290786/141‑1 (MQ=255) aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGGTATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGATATATc < 2:263930/141‑1 (MQ=255) aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGATATATc < 2:267730/141‑1 (MQ=255) aCTAGCCGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGTTTTCTGAAGTATTTAAGCCACATACTAAAAACGAAGCAGATGAAATATTTATAGAAAAAAAAAAAAAAACTACGCCAGCAATTTGTGATATATc < 2:1399/141‑1 (MQ=255) | TCCGACATACATAATGAATCTATCGATAAACAAATGGAAGCTAAAGCGCATGAAACAGCGCAAAATACAGATTTAAAAAACGAAGCAAGAAGTTTATTTGATAATGCAACCAAATCAATCGGTAGACTAGCGGGCAATGATGAAAGCTTAAATCTTAATTTAAAAGATATGCTTTCTGAAGTATTTAAGCCGCATACTAAAAACGAAGCAGATGAAATATTTATAGCGGGTACTGCTAAAACTACGCCAGCAATTTGTGACATATC > USA300TCH1516_ALE/277296‑277561 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |