Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,082,821 | T→G | 12.9% | intergenic (+150/+56) | loiP → / ← speB | Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,082,821 | 0 | T | G | 12.9% | 88.3 / 7.2 | 31 | intergenic (+150/+56) | loiP/speB | Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase |
Reads supporting (aligned to +/- strand): ref base T (14/13); new base G (2/2); total (16/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 = 6.88e-01 |
GACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGT > NC_000913/3082721‑3082962 | gACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgccttat > 1:286442/1‑107 (MQ=37) gACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgccttat < 2:286442/111‑5 (MQ=37) aCATGACTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc < 2:97604/124‑1 (MQ=255) aCATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc > 1:97604/1‑124 (MQ=255) cATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCa < 1:310664/101‑1 (MQ=255) cATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCa > 2:310664/1‑101 (MQ=255) caTTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt < 2:471964/117‑1 (MQ=255) caTTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt > 1:471964/1‑117 (MQ=255) ttGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgcctta > 2:775055/1‑83 (MQ=21) ttGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCATACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgcctta < 1:775055/86‑4 (MQ=21) atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt > 3:518272/1‑87 (MQ=255) atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt < 4:518272/87‑1 (MQ=255) atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCAtt > 2:945979/1‑148 (MQ=255) ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat < 3:671290/133‑1 (MQ=255) ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat > 4:671290/1‑133 (MQ=255) cgcATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAata > 3:502437/1‑132 (MQ=255) cgcATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAata < 4:502437/132‑1 (MQ=255) cATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagcgcc < 1:945979/148‑1 (MQ=255) tGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGc > 3:300992/1‑148 (MQ=255) cGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGAtt < 4:300992/148‑1 (MQ=255) tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat < 4:876377/96‑1 (MQ=255) tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAat > 3:876377/1‑96 (MQ=255) tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTc < 2:819405/108‑1 (MQ=255) tCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTc > 1:819405/1‑108 (MQ=255) aTAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAAtata > 2:562831/1‑95 (MQ=255) aTAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAAtata < 1:562831/95‑1 (MQ=255) tAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATc > 3:288406/1‑148 (MQ=255) cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGc < 4:72603/67‑1 (MQ=255) cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGc > 3:72603/1‑67 (MQ=255) cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc < 2:706588/81‑1 (MQ=255) cGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc > 1:706588/1‑81 (MQ=255) ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTgg < 1:301093/54‑1 (MQ=255) ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTgg > 2:301093/1‑54 (MQ=255) ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGt < 2:424144/143‑1 (MQ=255) ttCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGt > 1:424144/1‑143 (MQ=255) tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc > 2:57416/1‑131 (MQ=255) tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTc < 1:57416/131‑1 (MQ=255) | GACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGT > NC_000913/3082721‑3082962 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |