Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 91,323 | G→A | E98K (GAA→AAA) | ftsL → | cell division protein FtsL |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 91,323 | 0 | G | A | 96.8% | 102.2 / ‑4.2 | 31 | E98K (GAA→AAA) | ftsL | cell division protein FtsL |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (13/17); total (13/18) | |||||||||||
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.81e-01 |
CTGGAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGA > NC_000913/91245‑91400 | cTGGAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGaaa < 2:20025/81‑1 (MQ=255) gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAgg < 1:20026/81‑1 (MQ=255) gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAgg < 2:20027/81‑1 (MQ=255) gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAg > 2:20074/1‑80 (MQ=255) cGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATc < 2:20028/81‑1 (MQ=255) gagaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCg > 1:20088/1‑81 (MQ=255) agaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGc > 2:20090/1‑81 (MQ=255) gaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGcc < 1:20029/81‑1 (MQ=255) gaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGcc < 1:20030/81‑1 (MQ=255) gaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGcc < 1:20031/81‑1 (MQ=255) aTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCa < 2:20032/81‑1 (MQ=255) gCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACg > 1:20069/1‑81 (MQ=255) tAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAG‑GGGGTGAAAAGGATCGCCACGGaa < 2:20033/79‑1 (MQ=255) aGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGaaa < 1:20034/80‑1 (MQ=255) gACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAg < 1:20035/81‑1 (MQ=255) aTTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTg > 2:20094/1‑81 (MQ=255) gAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCa > 1:20077/1‑80 (MQ=255) tcgcAACCTGATCCTTGACGAGAATGCGCTCAGCGACCACAGCCGGGTGAAAAGGATCGCCGCGGAAAAGCTGCAAATgca < 2:20036/80‑1 (MQ=255) cgcAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATgcag > 1:20071/1‑81 (MQ=255) cTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGtt < 2:20037/81‑1 (MQ=255) tGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGtt > 1:20096/1‑80 (MQ=255) tCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATc < 1:20039/81‑1 (MQ=255) tCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATc < 2:20038/81‑1 (MQ=255) agagAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTcaca < 1:20040/81‑1 (MQ=255) gaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGa < 2:20041/81‑1 (MQ=255) aTGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGaaa < 1:20042/81‑1 (MQ=255) aTGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGaaa < 2:20043/81‑1 (MQ=255) gcgcTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCTAATGCAGCATGTTGATCCGTCACAAGaaaa > 2:20089/1‑80 (MQ=255) gACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTg > 2:656269/1‑81 (MQ=255) cATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCa > 2:20080/1‑80 (MQ=255) cATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCa > 1:20085/1‑80 (MQ=255) gTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGa > 1:20093/1‑81 (MQ=255) | CTGGAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGA > NC_000913/91245‑91400 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |