| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,680,185 | 1 | . | T | 21.1% | 45.0 / 12.7 | 19 | intergenic (+237/+295) | cbpD/USA300HOU_2537 | possible secretory antigen/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base . (9/6); new base T (2/2); total (11/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTAATTGCTTTAATA > CP000730/2680055‑2680309 | ttGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAAttttt > 2:43527/1‑141 (MQ=255) tataTTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGtt < 1:61565/141‑1 (MQ=255) attcctatctatcATAAGCTATATATAATATTTTATGGTAATAGTGAAATAATACGGTATAATTTAACGTACGCTTTCGATGTAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCAATAGTAaaaaaa < 2:312704/129‑1 (MQ=255) cATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAaata > 1:20193/1‑141 (MQ=255) tatataAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTga < 1:192526/141‑1 (MQ=255) taAAATTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGAT‑TTTTTTCATCATCATCAGTTaaaaaatcag < 2:69854/122‑5 (MQ=255) aaaaTTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTG‑ATTTTTTTCATCATCATCAGTTaaaaaatcaga > 2:192923/1‑117 (MQ=255) aaaaTTTTATGTTTTTGTTTAAGCAATCCACCATAACTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTG‑ATTTTTTTCATCATCATCAGTTaaaaaatcaga < 1:192923/122‑6 (MQ=255) tattatTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGCCCTGTTTAAATtgat > 1:282684/1‑141 (MQ=255) aataatTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTG‑TTAAATTGATGACTCAAAAc > 1:28713/1‑141 (MQ=255) aatTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCt < 2:300768/141‑1 (MQ=255) taAGTTAACATACGTTTTCGATATAAACCTTGTCCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTGTTCGTAATTAGGTGAGAAAAGAATCGATTGGCTGTTGTGAATCCTGTTTAGATGGATGTCTCAAAACGTTTATcgaac > 2:5101/1‑137 (MQ=255) taACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTT‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCt > 2:116598/1‑141 (MQ=255) aaCTTAACATACGTTTTCGATACAAACCTTTTTCTAAATCACTTACAATTTCTTGAT‑TTTTTTCATCATCATCAGTTaaaaaat > 1:91095/1‑83 (MQ=255) gTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGttt > 1:194662/1‑141 (MQ=255) ttttCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTa > 1:284356/1‑141 (MQ=255) ttCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTAAt < 1:218491/141‑1 (MQ=255) cGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTAATTg > 1:273688/1‑141 (MQ=255) aCCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTAATTGCTTTAata < 2:235264/141‑1 (MQ=255) | TTGAACGTATATTTCTTCTTCAACTATTATTTCATCTTTAGCATAATCTATATATAAAATTTTATGCTATTATTTAAATAATTCGCTATAACTTAACATACGTTTTCGATATAAACCTTGTTCTAAATCTC‑AATAATTTTTTG‑C‑TGTTTTCATCGTCATTAGTTAAAAAAATAATTTAACTGAGTTTTTGAGACCTGTTTAAATTGATGACTCAAAACCTTTATCCCCTTTTCACTCGGTTTAATTGCTTTAATA > CP000730/2680055‑2680309 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |