Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 320,471 | G→A | 12.1% | intergenic (‑140/‑138) | rmlC ← / → aguA | dTDP‑4‑dehydrorhamnose 3,5‑epimerase/agmatine deiminase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 320,471 | 0 | G | A | 12.1% | 82.3 / 3.9 | 33 | intergenic (‑140/‑138) | rmlC/aguA | dTDP‑4‑dehydrorhamnose 3,5‑epimerase/agmatine deiminase |
Reads supporting (aligned to +/- strand): ref base G (16/13); new base A (2/2); total (18/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01 |
GATGTTCATAGCGGCTCGCTGGCCCGTCGTTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATCACAGAGGTCACCCATGAAAACCCT > NC_002947/320323‑320619 | tattttCATAGCCGCTCTCTGGCCCGACGTTCAAGTTTTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCACGTAGCTGGAACAGGATTCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTGTCACCAGACCAGCTCTGATAGgt < 2:94903/146‑1 (MQ=255) tgTTCATAGCGGCTCGCTGGCCCGTCGTTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGgtgt < 2:722642/150‑1 (MQ=255) tAGCGGCTCGCTGGCCCGTCGTTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTg > 1:717747/1‑150 (MQ=255) cTCGCTGGCCCGTCGTTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCg < 2:686902/150‑1 (MQ=255) cgtcgtTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTc < 2:340039/150‑1 (MQ=255) aaGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGATGTCGATTGTTATCGATTTAAATTCGCCGAATg > 1:486051/1‑150 (MQ=255) ttGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGcgcg < 2:104731/150‑1 (MQ=255) aCCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGaa > 1:154247/1‑150 (MQ=255) ccTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGAttt > 1:980964/1‑150 (MQ=255) cTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGgt > 1:193536/1‑98 (MQ=255) cTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGgt < 2:193536/98‑1 (MQ=255) cGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGc > 1:691619/1‑150 (MQ=255) cGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGc < 1:735730/150‑1 (MQ=255) tCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAata < 1:234173/150‑1 (MQ=255) gcgcgTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCg > 2:784451/1‑150 (MQ=255) gTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGatat < 2:154247/150‑1 (MQ=255) tGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGa > 2:1038589/1‑78 (MQ=255) tGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGACCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGa < 1:1038589/78‑1 (MQ=255) tAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATc > 2:192588/1‑138 (MQ=255) tAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATc < 1:192588/138‑1 (MQ=255) aCAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGtt > 2:775499/1‑150 (MQ=255) aTGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATaaa < 2:717747/150‑1 (MQ=255) cATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTg > 2:460586/1‑150 (MQ=255) tAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCTTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGCTTTTCAGGCTTCAAATATCGCTATATATCGGTTATAAATTGAc > 2:172156/1‑150 (MQ=255) aGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAAt > 1:660947/1‑88 (MQ=255) aGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAAt < 2:660947/88‑1 (MQ=255) gtgGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGATGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGc < 2:486051/150‑1 (MQ=255) tATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGATGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGtt < 1:316325/129‑1 (MQ=255) tATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGATGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGtt > 2:316325/1‑129 (MQ=255) tAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAAttcttc < 2:345758/89‑1 (MQ=255) tAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAAttcttc > 1:345758/1‑89 (MQ=255) aGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATAcc < 1:775499/150‑1 (MQ=255) gTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGcc < 2:980964/150‑1 (MQ=255) tttCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATcac > 2:1194963/1‑150 (MQ=255) tttCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATcac > 2:129252/1‑150 (MQ=255) tttCACCAGACAAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATcac > 2:933900/1‑150 (MQ=255) aTAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATCACAGAGGTCACCCATGAAAAc > 2:99544/1‑150 (MQ=255) ggTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATCACAGAGGTCACCCATGAAAACCCt < 1:99544/150‑1 (MQ=255) | GATGTTCATAGCGGCTCGCTGGCCCGTCGTTCAAGTATTGACCCTAGTCGGCCTTGCCGGAAATTTCCGCGCGTAGCTGTAACAGGATGCATAAGGTGTGGCTATGGTAGACAGCAGGCCGTTCCTTTCACCAGACCAGCGCTGATAGGTGTCGATTGTTATCGATTTAAATTCGCCGAATGGCGCGGAATTCTTCGATTTTCAGGCTTCAAATATCGATATATATCGGTTATAAATTGACAAACGCCCCAGCATACCCGGCTCAGGGCCATCACAGAGGTCACCCATGAAAACCCT > NC_002947/320323‑320619 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |