| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,957,951 | T→A | Q61L (CAG→CTG) | torY ← | TMAO reductase III (TorYZ), cytochrome c‑type subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,957,951 | 0 | T | A | 96.4% | 185.1 / ‑1.3 | 57 | Q61L (CAG→CTG) | torY | TMAO reductase III (TorYZ), cytochrome c‑type subunit |
| Reads supporting (aligned to +/- strand): ref base T (2/0); new base A (22/31); total (26/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.86e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01 | |||||||||||
ATGATAAATATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTA > NC_000913/1957854‑1958044 | aTGATAAATATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAgg < 2:891545/100‑1 (MQ=255) tGATAAATATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGt < 1:1098879/100‑1 (MQ=255) aTAAATATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTc < 2:1212962/100‑1 (MQ=255) aaTATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGaa > 2:1248846/1‑100 (MQ=255) tCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGt < 1:556550/100‑1 (MQ=255) cTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTg > 1:566792/1‑100 (MQ=255) cTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTg > 2:1248807/1‑100 (MQ=255) cTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTg < 2:1261749/100‑1 (MQ=255) tttAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTgg < 1:210202/100‑1 (MQ=255) tttAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTgg < 2:59297/100‑1 (MQ=255) ttAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTggn < 2:467322/100‑2 (MQ=255) ttccgatcTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGAc < 1:736576/93‑1 (MQ=255) tAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGAc > 1:1498876/1‑100 (MQ=255) gCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGtt > 1:703825/1‑99 (MQ=255) ccTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTcc > 1:308068/1‑100 (MQ=255) tttAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCt < 1:568110/100‑1 (MQ=255) tttAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCt < 1:1247962/100‑1 (MQ=255) ttAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGATcggaa > 2:736574/1‑96 (MQ=255) aGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTcc < 1:823683/100‑1 (MQ=255) aGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTcc < 1:5874/100‑1 (MQ=255) gCAAATAAATAATCCATCCCTGACTTTGGAATATGACAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCt > 2:827615/1‑100 (MQ=255) cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCCGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTc > 2:1401321/1‑100 (MQ=255) cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTc < 1:685276/100‑1 (MQ=255) cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTc < 2:139512/100‑1 (MQ=255) cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTc > 1:799435/1‑100 (MQ=255) aaataaatAATCCATCCCTGACTTTAGAATATGACAATCGGCACATTACACACATATCCCTTTAAGATACAAAAAGTAGACAGTTCCCTTATATTtctc > 2:143860/1‑99 (MQ=255) aataaataATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCCGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCaa > 1:1111601/1‑100 (MQ=255) aaataaTCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAg > 2:1220726/1‑100 (MQ=255) aataatCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAgg > 1:764249/1‑100 (MQ=255) ataatCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGt > 1:996262/1‑100 (MQ=255) taatCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGtt < 1:854777/100‑1 (MQ=255) taatCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTGCCCTGATATTCCTCAAAAGGtt < 2:934399/100‑1 (MQ=255) aatCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGttt < 1:1395447/100‑1 (MQ=255) atccatccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTa > 1:1254477/1‑100 (MQ=255) catccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTc > 2:1460461/1‑100 (MQ=255) ccTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCg < 2:1411229/100‑1 (MQ=255) gACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGt > 2:160503/1‑100 (MQ=255) tttGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGc < 1:1376776/100‑1 (MQ=255) ttGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCa > 1:477267/1‑100 (MQ=255) ttGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCa < 2:307076/100‑1 (MQ=255) ttGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCa < 2:1407143/100‑1 (MQ=255) ggAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAgg > 1:877213/1‑100 (MQ=255) ggAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAgg < 1:962976/100‑1 (MQ=255) tatGGCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAg < 1:1460469/100‑1 (MQ=255) atGGCAATCGGCACGTTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAgg < 2:53886/100‑1 (MQ=255) ggCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCa > 2:195951/1‑100 (MQ=255) gCAATCGGCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCaa > 2:199472/1‑100 (MQ=255) ggCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCt > 1:484460/1‑100 (MQ=255) gCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTTGACAGTTCCCTGATATTCCTCAAAGGGTTTACTCATCGCGTGGCAGGAAAGGCAAAATGCt > 1:189418/1‑99 (MQ=255) gCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTg < 1:701545/100‑1 (MQ=255) gCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTg < 2:70373/100‑1 (MQ=255) gCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTg < 1:1248851/100‑1 (MQ=255) gCACATTCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTg < 1:904826/100‑1 (MQ=255) ttCCGCACGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATcc > 2:1096383/1‑100 (MQ=255) cGTATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGttt < 1:160503/100‑1 (MQ=255) tATCCCTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTa < 1:1248812/100‑1 (MQ=255) ccTTTCAGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTa > 1:193814/1‑100 (MQ=255) | ATGATAAATATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTA > NC_000913/1957854‑1958044 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |