| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,929,867 | 0 | G | T | 65.4% | 46.1 / 56.7 | 52 | intergenic (‑950/‑218) | USA300HOU_1794/USA300HOU_1796 | hypothetical protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (8/10); new base T (19/15); total (27/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTA > CP000730/1929727‑1929999 | tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑g‑g < 2:251992/141‑1 (MQ=255) tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTg < 1:266673/141‑1 (MQ=255) aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGa > 1:14015/1‑141 (MQ=255) aaTCTATCAGCTAGAATCATATAACTTATCAGAAGGAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAttt < 2:24253/141‑1 (MQ=255) aTCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGtt < 2:259050/141‑1 (MQ=255) tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGc < 1:1325/141‑1 (MQ=255) cTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATCGTAGAGTCGCCTATCTCTCAGGCGTGAATTGAGACTCAGAGAGCAGGTG‑TA‑TAA‑G‑GGGATGCTTATTTTCGTTCACATCATAGCACCAGtca > 1:195005/1‑141 (MQ=255) cTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtca < 1:85138/141‑1 (MQ=255) cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg > 1:81655/13‑141 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCa > 2:218013/11‑125 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCa < 1:218013/115‑1 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg < 2:156497/131‑1 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg > 1:201363/11‑141 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGGATGTG‑CA‑TTACT‑GTGCTAATCATTTTCGTTCACATCATATCACCAGTGATCACTGGCTGTGctg > 2:181905/8‑139 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGTGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg > 2:252388/8‑139 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTgcgctg > 2:42269/8‑139 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg > 2:221164/8‑139 (MQ=255) tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt < 1:258287/135‑4 (MQ=255) tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt < 2:17244/135‑4 (MQ=255) tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt < 1:117423/135‑4 (MQ=255) tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt < 1:515/135‑4 (MQ=255) agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt > 1:207261/6‑137 (MQ=255) aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt > 1:185923/1‑141 (MQ=255) atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt < 2:62870/140‑1 (MQ=255) tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTTTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta > 2:242101/1‑140 (MQ=255) ggCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt < 2:103188/140‑1 (MQ=255) gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat < 2:179284/141‑3 (MQ=255) cTTGTATATTTGCGGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTGTCGTTCACATCATGGCACCAGTCATCAGTGGCTGTGCCATGGCGttt > 2:57351/1‑141 (MQ=255) cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt < 1:219500/141‑4 (MQ=255) cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt > 2:128042/1‑141 (MQ=255) ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta < 1:147782/141‑5 (MQ=255) ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta < 1:133551/141‑5 (MQ=255) cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGCCAATT‑‑GACGAAGAGAGGAGGTGCCT‑TAC‑G‑TTGCTAATCAGTTTAGTTCACATCAACGCACCAGTCATTATTGGCTGTGCTGTTGCGTATTATACTTATTg > 2:114549/5‑141 (MQ=255) acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg < 1:260840/138‑1 (MQ=255) taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa < 1:172826/138‑1 (MQ=255) atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac > 1:73171/2‑140 (MQ=255) agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa > 1:204639/1‑138 (MQ=255) gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata < 2:65554/141‑7 (MQ=255) tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa > 1:58534/1‑134 (MQ=255) cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat > 1:234277/1‑133 (MQ=255) cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat > 1:38588/1‑133 (MQ=255) cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat > 2:101639/1‑133 (MQ=255) gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt < 2:73171/141‑10 (MQ=255) aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata > 2:69785/1‑128 (MQ=255) aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTGCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata > 2:52756/1‑128 (MQ=255) tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag > 2:134766/1‑127 (MQ=255) ctaaattgACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc > 1:135896/8‑86 (MQ=255) ctaaattgACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc < 2:135896/79‑1 (MQ=255) aGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt < 2:164721/141‑1 (MQ=255) ggTG‑TA‑TGA‑G‑GTGATGGGTATGTTCGCTCAGATGTGAGCACCAGTGAGCGGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc < 2:213846/141‑1 (MQ=255) tg‑tA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa > 2:227049/1‑141 (MQ=255) aa‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc < 1:169303/68‑1 (MQ=255) aa‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc > 2:169303/1‑68 (MQ=255) | TGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTA > CP000730/1929727‑1929999 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |