Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 2,343,438 | C→T | 100% | intergenic (+114/+250) | USA300HOU_2200 → / ← alsD | hypothetical protein/acetolactate decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 2,343,438 | 0 | C | T | 87.0% | 62.9 / 4.7 | 31 | intergenic (+114/+250) | USA300HOU_2200/alsD | hypothetical protein/acetolactate decarboxylase |
Reads supporting (aligned to +/- strand): ref base C (1/3); new base T (13/14); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.07e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ATACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTAA > CP000730/2343304‑2343572 | aTACAAATTGTAGGTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATTGt > 2:331616/1‑141 (MQ=255) tACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTa > 2:5393/1‑141 (MQ=255) aCAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATGGTaa > 2:328013/1‑141 (MQ=255) aTTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAAttt < 1:250567/141‑1 (MQ=255) ataatTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTa < 2:242770/141‑1 (MQ=255) aatTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTAcc > 1:318308/1‑141 (MQ=255) aaTAGTTTTATAGGGCAAGCGCTGAAAAATTTGGGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTta < 2:287876/141‑1 (MQ=255) aaTAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTta > 2:103770/1‑141 (MQ=255) tAGTTTTATAGGGCAAGAACCTCAAAGGGTGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTAGCACAATTtata > 2:322822/1‑141 (MQ=255) tataGGGCAAGCGCTGAAAAATTGGCGTTTTTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTaa < 1:160120/141‑1 (MQ=255) gggCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACAt < 1:284284/141‑1 (MQ=255) aaaaaTTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAAtct < 2:218750/141‑1 (MQ=255) ggCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTATCGTAGTAACAGATTGGCCTATGTAATAGCAAATTTTATAGGAGTGTGGCACCCAAATAAAAATCGAAACATACCAGCCAATCTCagaaaa > 2:229408/1‑136 (MQ=255) tttGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTgg > 1:98374/1‑141 (MQ=255) tttttACTTATTGACCTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTACATGAGTGTGTACCACAATTTATAATCTAAACATAACCGCAAATCTCAAATATCTGGAACCta < 2:154730/141‑1 (MQ=255) tcttAGGTGTTGAACTTGTAAACTATTGTTCTTCCCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAACTTTATATGAATGTCCACCACAATTTATAATCTAAACATAACAGCAAACCTCAAATATCTGAAACCtaa < 2:273761/139‑1 (MQ=255) ttAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCtaata > 1:15214/1‑141 (MQ=255) gTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAAt < 1:333869/141‑1 (MQ=255) aTTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATACATccc > 2:277692/1‑141 (MQ=255) tGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCt < 1:5393/141‑1 (MQ=255) gAACTTGTAAATTATTGTTGTTTTCTTCTCGTAGGAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAAGCCCCTc < 2:22855/141‑1 (MQ=255) gAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTc < 2:130907/141‑1 (MQ=255) ttgttTGCTTCTCGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGa > 1:134715/1‑141 (MQ=255) tGCTTCTCGTAGTAACCGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGa < 1:307396/141‑1 (MQ=255) tctcGTAGTAACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGattta < 2:98374/141‑1 (MQ=255) aagtaACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGAtttatttta < 1:104374/140‑1 (MQ=255) taACAGATTGGCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTttatta > 1:56694/1‑141 (MQ=255) cAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACtt < 1:1510/141‑1 (MQ=255) aTTGGCCTATATAATACTACAGTTTATATGAGTGTGTACCACACTTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTCTTTTATTACTTAGa > 1:171310/1‑141 (MQ=255) gCCTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTa < 1:264257/141‑1 (MQ=255) ccTATATAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTaa > 1:34007/1‑141 (MQ=255) | ATACAAATTGTAGTTAAATAATAATTTGAAAATAGTTTTATAGGGCAAGCGCTGAAAAATTGGCGTTTGTCTTTTTAGTTATTGAACTTGTAAATTATTGTTGTTTGCTTCTCGTAGTAACAGATTGGCCTATACAATAGTAAATTTTATATGAGTGTGTACCACAATTTATAATCTAAACATAACAGCAAATCTCAAATATCTGGAACCTAATAAATCCCCTCACTACTGCTAATAGTGAAGGGATTTATTTTATTACTTAGACGTAA > CP000730/2343304‑2343572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |