Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,059,436 | (G)6→7 | coding (181/2547 nt) | torA → | trimethylamine N‑oxide reductase 1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,059,430 | 1 | . | G | 92.6% | 85.7 / 0.5 | 27 | coding (175/2547 nt) | torA | trimethylamine N‑oxide reductase 1 |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (13/12); total (14/13) | |||||||||||
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.94e-01 |
CTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑T‑G‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCG > NC_000913/1059315‑1059557 | cTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCcgcgc > 2:151370/1‑133 (MQ=255) ttAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCcgcgcg < 2:89879/133‑1 (MQ=255) tgaCCGTCGCTGGGATGCTGGGGGCGGCCTTGGTAACGCCGAGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGGCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGa < 2:291695/131‑1 (MQ=255) gCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCGTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGTGACGGTGAAg > 2:253596/1‑133 (MQ=255) gCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgt < 2:16794/133‑1 (MQ=255) ccGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgtg > 1:66644/1‑133 (MQ=255) ccGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgtg < 1:33765/133‑1 (MQ=255) gTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTggc > 2:497780/1‑133 (MQ=255) cATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTggcgg < 1:151370/133‑1 (MQ=255) cATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑T‑G‑GG‑GCGGCTATCCGCGCGACGGTGAAGGATTGTCGGTTTGTggcgg > 2:161903/1‑133 (MQ=255) cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg > 1:395797/1‑133 (MQ=255) cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg > 2:366624/1‑133 (MQ=255) cgcgACGTACGACTGAGGTGCAAGCGGCGAATGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg < 1:161903/133‑1 (MQ=255) cgACGTGCCACTGCGTCGCAAGTGGCGACTCACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa < 1:173386/133‑1 (MQ=255) tGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGGCGCGTTGGGGCGGCAAGACCGGTCGAACTGGa > 2:173386/1‑133 (MQ=255) gACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgat > 2:524501/1‑133 (MQ=255) gACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgat < 1:313952/133‑1 (MQ=255) gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccgga > 1:440087/1‑133 (MQ=255) tCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccggattgc > 1:14243/1‑133 (MQ=255) tctcGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg < 1:514829/133‑1 (MQ=255) tctcGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg < 1:366624/133‑1 (MQ=255) tctcGAAAGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg < 2:66644/133‑1 (MQ=255) tctcGAAAGAGGGCATTCTTACCGGGTCGCAA‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg < 2:440087/133‑1 (MQ=255) tcGAAAGAGGACATTTTTACCGGGCCGCCC‑T‑GTGG‑G‑GGCTAACCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATtgccggattgccggat < 2:185129/133‑1 (MQ=255) aaaGAGGGCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCAc > 1:314642/1‑133 (MQ=255) gCATTCTTACCGGGTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTAcaca > 2:112585/1‑133 (MQ=255) ttCTTACCGGGTCGCACTT‑G‑GG‑C‑GGCTATCGGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAAcg < 1:212278/133‑1 (MQ=255) gggTCGCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGt > 1:355416/1‑133 (MQ=255) gCAC‑TGG‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg < 2:314642/133‑1 (MQ=255) gCAC‑T‑G‑GGCG‑GGCTATCCACGCCACGGTGACGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg < 1:175703/133‑1 (MQ=255) | CTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCAC‑T‑G‑GG‑G‑GGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCG > NC_000913/1059315‑1059557 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |