Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,778,934 | (T)6→5 | intergenic (‑564/+187) | BW25113_RS25705 ← / ← BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,778,929 | 0 | T | . | 97.6% | 168.8 / ‑1.5 | 42 | intergenic (‑559/+192) | BW25113_RS25705/BW25113_RS13825 | DUF5507 domain‑containing protein/tRNA‑Ile |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (20/21); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01 |
AAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC > NZ_CP009273/2778840‑2779012 | aaGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAt < 1:832784/90‑1 (MQ=255) ataaATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGt < 2:332634/90‑1 (MQ=255) ataaATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGt < 1:712396/90‑1 (MQ=255) aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACt > 1:157742/1‑90 (MQ=255) aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACt > 2:999900/1‑90 (MQ=255) aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACt > 2:5132/1‑90 (MQ=255) aTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACtt > 1:410790/1‑90 (MQ=255) tGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACttt > 2:765641/1‑90 (MQ=255) gATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACtttt < 2:50657/90‑1 (MQ=255) ataataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAt < 2:837881/73‑1 (MQ=255) ataataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAt > 1:837881/1‑73 (MQ=255) aataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGaa > 1:443294/1‑90 (MQ=255) taCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATg < 2:206534/90‑1 (MQ=255) taCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATg < 2:157742/90‑1 (MQ=255) ctctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAg > 2:824966/1‑90 (MQ=255) tctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGt > 2:298271/1‑90 (MQ=255) tctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGt < 1:49938/90‑1 (MQ=255) ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGt > 2:907654/1‑90 (MQ=255) ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTa > 2:178448/1‑90 (MQ=255) gCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTg < 2:735207/90‑1 (MQ=255) tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGAt < 2:443294/90‑1 (MQ=255) tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATgg > 1:896339/1‑90 (MQ=255) tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATgg > 1:727402/1‑90 (MQ=255) cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa < 1:242518/90‑1 (MQ=255) cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa < 1:768687/90‑1 (MQ=255) cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa < 1:418326/90‑1 (MQ=255) ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa > 1:497127/1‑90 (MQ=255) ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATaaa > 1:351216/1‑90 (MQ=255) aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc < 2:11431/90‑1 (MQ=255) ccACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGc > 2:286057/1‑90 (MQ=255) acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa > 1:361401/1‑72 (MQ=255) acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa < 2:361401/72‑1 (MQ=255) tGTAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGcc > 2:696315/1‑90 (MQ=255) tAAGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc < 2:443265/90‑1 (MQ=255) aaGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGAt > 2:10318/1‑61 (MQ=255) aaGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGAt < 1:10318/61‑1 (MQ=255) aaGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc > 2:512474/1‑71 (MQ=255) aaGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc < 1:512474/71‑1 (MQ=255) aaGCAATTTGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACg > 2:725817/1‑90 (MQ=255) tttGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaa < 1:298271/90‑1 (MQ=255) ttGCCCGCTTGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa < 2:41947/90‑1 (MQ=255) tGGCATAGCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGcc < 1:174894/90‑1 (MQ=255) gCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTc > 1:171225/1‑90 (MQ=255) gCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTc < 2:564972/90‑1 (MQ=255) gCAGGCA‑TTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTc < 1:999900/90‑1 (MQ=255) | AAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC > NZ_CP009273/2778840‑2779012 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |