| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,093,803 | (A)6→5 | coding (928/1359 nt) | ycdT → | diguanylate cyclase, membrane‑anchored |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,093,798 | 0 | A | . | 98.0% | 221.6 / ‑2.2 | 50 | coding (923/1359 nt) | ycdT | diguanylate cyclase, membrane‑anchored |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (21/28); total (22/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01 | |||||||||||
ACGAGTAACAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGAT > NC_000913/1093700‑1093888 | aCGAGTAACAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaa > 2:1056462/1‑100 (MQ=255) gAGTAACAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaa < 1:772126/100‑1 (MQ=255) aGTAACAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaa < 1:1065866/100‑1 (MQ=255) aCAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaacgcc < 1:474432/100‑5 (MQ=255) aCAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaacgcc > 1:475413/1‑96 (MQ=255) aCAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaacgcc > 2:597851/1‑96 (MQ=255) cAAAGAACATTGCACATCGCGATCCCTTAACCTATATATGTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCt < 2:916113/100‑1 (MQ=255) cAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCt > 2:334957/1‑100 (MQ=255) aaCATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTg < 2:1274879/100‑1 (MQ=255) aaCATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTg < 1:511884/100‑1 (MQ=255) ttGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTc < 2:105978/100‑1 (MQ=255) gCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCAt > 2:1233544/1‑100 (MQ=255) cacaTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATg > 2:1012865/1‑100 (MQ=255) aTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGAtt > 2:282739/1‑100 (MQ=255) aTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGAt < 2:110291/99‑1 (MQ=255) gcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATg < 2:1275094/100‑1 (MQ=255) gcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATg > 2:719352/1‑100 (MQ=255) gcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATg > 2:860260/1‑100 (MQ=255) gcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATg < 2:1275131/100‑1 (MQ=255) gcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATg > 1:185422/1‑100 (MQ=255) gATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGa < 1:479470/100‑1 (MQ=255) gATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGa < 1:195303/100‑1 (MQ=255) gATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGa < 1:975476/100‑1 (MQ=255) gATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGa < 1:148270/100‑1 (MQ=255) aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGat < 1:816806/100‑1 (MQ=255) aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATAAGCATCAGCCC‑AAAAACGCTTTATTGCGTCATGATTATGGat < 2:757974/100‑1 (MQ=255) ccTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATc > 1:1086730/1‑100 (MQ=255) tAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGAc > 2:949955/1‑100 (MQ=255) aaCCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGAcc < 1:1172536/100‑1 (MQ=255) tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCa < 2:523685/100‑1 (MQ=255) tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCa < 2:431096/100‑1 (MQ=255) tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAg > 2:1212140/1‑100 (MQ=255) tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAg < 1:1212151/100‑1 (MQ=255) aaCAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGt > 1:1206668/1‑100 (MQ=255) aCAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTc < 2:675126/100‑1 (MQ=255) gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCTTCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAAc < 2:877600/100‑1 (MQ=255) gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAAc > 2:1069051/1‑100 (MQ=255) gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAAc < 1:646106/100‑1 (MQ=255) gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAAc < 1:334960/100‑1 (MQ=255) attTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACAcc > 2:1119386/1‑100 (MQ=255) tttttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTTCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCt > 1:174158/1‑100 (MQ=255) tttttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCt < 2:1264997/100‑1 (MQ=255) ttttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTg < 1:74225/100‑1 (MQ=255) tttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTgg > 1:239236/1‑100 (MQ=255) tttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTgg > 2:1235857/1‑100 (MQ=255) ttttAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTggg < 1:479975/100‑1 (MQ=255) tttAATGAACTTACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTgggg < 1:107415/100‑1 (MQ=255) tgattGAACTGTCAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTgggg > 1:572008/5‑99 (MQ=255) tAATGAACTGACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGc > 1:197945/1‑100 (MQ=255) gACAGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTg > 1:709929/1‑100 (MQ=255) aGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCg < 2:197943/100‑1 (MQ=255) aGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCg < 1:597855/100‑1 (MQ=255) aGTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCg < 2:1086720/100‑1 (MQ=255) gTTCAATCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGa > 2:432663/1‑100 (MQ=255) aaTCAGCATCAGCCC‑AAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAg > 1:133837/1‑100 (MQ=255) catcagcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGAt > 2:149188/1‑100 (MQ=255) | ACGAGTAACAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGAT > NC_000913/1093700‑1093888 |
| 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: ‑ |