Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP000730 | 195,286 | T→A | intergenic (‑97/‑82) | USA300HOU_0183 ← / → USA300HOU_0184 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 195,286 | 0 | T | A | 82.2% | 116.7 / 15.4 | 45 | intergenic (‑97/‑82) | USA300HOU_0183/USA300HOU_0184 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (5/3); new base A (15/22); total (20/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.35e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
GTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAG > CP000730/195151‑195422 | gTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTaaaaaa < 1:307639/141‑1 (MQ=255) tCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAat > 2:351906/1‑141 (MQ=255) cTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCtt < 2:323511/141‑1 (MQ=255) atAAACTTGCTACGAAAGATAATATATAAAAAAAGATAGTCCTCTATATATATGTTTATTAACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACATAATTTTATCAAAAAAAAAAAAAAAATATCtt < 2:289544/140‑1 (MQ=255) tAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACAAGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTa < 1:246848/141‑1 (MQ=255) tAAACTTGCTACGAAAGATTAATTAAAAATAATGATAGCCTCATATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGTATTTTATCTAAAAATAAAAAAAAATATCTTa < 2:352499/141‑1 (MQ=255) aaCTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAat < 1:304595/141‑1 (MQ=255) aCTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTAACTAAAAAATAAAAAAAATATCTTAata < 2:300226/141‑1 (MQ=255) aCAGGTTCTTTATCCAGAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAAt < 2:312703/141‑1 (MQ=255) aTTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATc > 2:307456/1‑129 (MQ=255) aTTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATc < 1:307456/129‑1 (MQ=255) ttCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGAt > 2:195917/1‑141 (MQ=255) tCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATAAAAAAAAAAATATCTTAATATCCTTGTAATCCGATa < 2:67328/141‑1 (MQ=255) ttATCCATAATGATAGCCGCCTATTTATATGGTTTTTTACTTATACCATAACATGCTTTTTATATTCTTTGAAAATATATTTTAAAGATCTTTATCTACATTTTAAAAAAAATCTGCTAATATAATTGTAAGCAGATaata > 2:53619/1‑139 (MQ=255) tATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGaa > 2:78615/1‑141 (MQ=255) tATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGaa > 2:166471/1‑141 (MQ=255) aaTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGt > 2:38000/1‑141 (MQ=255) aTAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAata < 1:391092/141‑1 (MQ=255) aGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATa < 2:263500/124‑1 (MQ=255) aGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATa > 1:263500/1‑124 (MQ=255) gCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATAttt > 2:8277/1‑141 (MQ=255) tatatataTGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCaa > 1:64537/1‑141 (MQ=255) atatatatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAAc > 1:86109/1‑141 (MQ=255) atatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAtt < 1:390705/141‑1 (MQ=255) tttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTtat < 1:258694/141‑1 (MQ=255) cTTATACCCTAACATGATTTTTATATTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc > 1:200792/1‑141 (MQ=255) ttATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCt < 1:24046/141‑1 (MQ=255) ataCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCtta > 2:400899/1‑141 (MQ=255) aTTTTTATATTTTTTTAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACattat > 2:35004/1‑141 (MQ=255) aTTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACattat > 1:395918/1‑141 (MQ=255) tttATACTCTTTGAAATTATATTTTACAGAATTTTATCTAAATATTTAAGAAAGTATGTTTATATCATTGTGATCCGAAAAGAATTATAGTAATATTTTTTCAATCATCGTTAGAGGCGGTATTAATAATGTCATTAgtta < 2:117374/141‑2 (MQ=255) tttATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTAtttt < 1:276956/141‑1 (MQ=255) tataCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTtta < 2:160269/141‑1 (MQ=255) taCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttatt > 2:372552/1‑141 (MQ=255) tctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAACATATCTTAATATCCTCATATGCGGATTAAAAATAGGGTAATAGATTGTCAACCATGTTAATAGGAGGGCTTATTAATTACATTATTTTTATTAGa > 1:143293/1‑141 (MQ=255) atatTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCCACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATc > 1:235789/1‑141 (MQ=255) atatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATc > 1:155102/1‑141 (MQ=255) tatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCt < 2:168359/141‑1 (MQ=255) tACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAtt < 2:86109/141‑1 (MQ=255) aCAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAttt < 2:38256/141‑1 (MQ=255) aCAGAATTTTATATAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAttt < 2:38251/141‑1 (MQ=255) aaTTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCa < 1:306787/141‑1 (MQ=255) aTCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACg < 1:351906/141‑1 (MQ=255) tCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa > 2:327619/1‑141 (MQ=255) aaTATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTTATTTTGCATCAACGAaaga < 1:200633/141‑1 (MQ=255) tatTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAaagaag < 1:372552/141‑1 (MQ=255) | GTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAG > CP000730/195151‑195422 |
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: ‑ |