Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 302,135 | G→A | Q147* (CAG→TAG) | paoA ← | PaoABC aldehyde oxidoreductase, 2Fe‑2S subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 302,135 | 0 | G | A | 98.1% | 164.6 / ‑4.3 | 53 | Q147* (CAG→TAG) | paoA | PaoABC aldehyde oxidoreductase, 2Fe‑2S subunit |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (21/31); total (21/32) | |||||||||||
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.88e-01 |
GGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCC > NC_000913/301999‑302270 | ggAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCTGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACa < 2:421315/139‑1 (MQ=255) ggAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCa < 2:557645/139‑1 (MQ=255) ggAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACa < 2:1294301/139‑1 (MQ=255) gAGCGGAAACCAAATCGACCGTGACGTGACTAGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACAt < 1:52312/139‑1 (MQ=255) gAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGt < 2:624149/139‑1 (MQ=255) cgtgacgtgaCTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGc < 2:944741/139‑1 (MQ=255) acgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAc > 1:1363830/1‑139 (MQ=255) cgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc < 1:1156920/139‑1 (MQ=255) cgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc < 2:57974/139‑1 (MQ=255) gacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAg < 2:916816/139‑1 (MQ=255) gacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAg < 2:387126/139‑1 (MQ=255) gacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAg < 1:572697/139‑1 (MQ=255) ggAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGt < 2:803704/112‑1 (MQ=255) ggAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGt > 1:803704/1‑112 (MQ=255) aaTGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCa > 2:9285/1‑139 (MQ=255) aaTGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCa < 1:9285/139‑1 (MQ=255) cTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGgtggtg > 1:90619/1‑139 (MQ=255) tGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGa < 1:1198920/139‑1 (MQ=255) aTCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGAtct < 1:1428565/139‑1 (MQ=255) ctTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGc < 1:25069/139‑1 (MQ=255) tAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGccccc < 2:90619/139‑1 (MQ=255) cTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATggg < 2:650204/79‑1 (MQ=255) cTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATggg > 1:650204/1‑79 (MQ=255) atgaGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGc < 2:509617/89‑1 (MQ=255) atgaGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGc > 1:509617/1‑89 (MQ=255) aGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAAtggtgg < 1:98684/108‑1 (MQ=255) aGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAAtggtgg > 2:98684/1‑108 (MQ=255) gCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAg < 1:252709/85‑1 (MQ=255) gCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAg > 2:252709/1‑85 (MQ=255) cAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCACCTACATGGGGTGAAGATTATCTGGCGAACCCAGGcc < 1:1120354/77‑1 (MQ=255) cAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCACCTACATGGGGTGAAGATTATCTGGCGAACCCAGGcc > 2:1120354/1‑77 (MQ=255) gCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATgggg < 2:36967/46‑1 (MQ=255) gCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATgggg > 1:36967/1‑46 (MQ=255) cGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATTGTGGTGATCTCGGccccc < 2:794909/95‑1 (MQ=255) cGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGccccc > 1:794909/1‑95 (MQ=255) cACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGccccc < 2:186717/93‑1 (MQ=255) cACTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGccccc > 1:186717/1‑93 (MQ=255) aCTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc > 1:1008014/1‑61 (MQ=255) aCTGGAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc < 2:1008014/61‑1 (MQ=255) ggAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc > 2:287095/1‑58 (MQ=255) ggAAGCCATCATGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc < 1:287095/58‑1 (MQ=255) catGCTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGt > 1:249746/1‑139 (MQ=255) cTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc > 1:1226345/1‑45 (MQ=255) cTTGATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAAcc < 2:1226345/45‑1 (MQ=255) gATAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCt > 2:1068059/1‑139 (MQ=255) aTAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGa < 1:752495/61‑1 (MQ=255) aTAAAGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGa > 2:752495/1‑61 (MQ=255) aGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAAtg > 2:577233/1‑50 (MQ=255) aGGCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAAtg < 1:577233/50‑1 (MQ=255) ggCCGCCTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCAc < 2:249746/139‑1 (MQ=255) gccTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCGTGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTAc < 2:596962/124‑1 (MQ=255) gccTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTAc > 1:596962/1‑124 (MQ=255) ccTACATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTcc > 1:336220/1‑139 (MQ=255) | GGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCC > NC_000913/301999‑302270 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |