Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,136,732 | A→T | 9.5% | L21Q (CTG→CAG) | fadE ← | medium‑long chain acyl‑CoA dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,136,732 | 0 | A | T | 9.5% | 120.0 / 4.9 | 42 | L21Q (CTG→CAG) | fadE | medium‑long chain acyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (19/19); new base T (2/2); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01 |
GTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAG > NC_002947/2136583‑2136880 | gTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCa < 1:874444/150‑1 (MQ=255) tCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCg > 1:350790/1‑150 (MQ=255) cACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGgcgc < 1:1205566/150‑1 (MQ=255) cACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGc > 2:802354/1‑150 (MQ=255) aCTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCTGGCGccctt < 2:838126/147‑4 (MQ=255) aCTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCTGGCGccctt > 1:838126/1‑144 (MQ=255) aCTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCg > 1:1184700/1‑150 (MQ=255) gCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGa < 1:387130/150‑1 (MQ=255) cGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAg < 2:883673/150‑1 (MQ=255) cGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAg < 1:181558/150‑1 (MQ=255) ggCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGAt < 1:673332/150‑1 (MQ=255) ttCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGc > 1:218527/1‑150 (MQ=255) ccAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCAccacca < 2:1065445/150‑1 (MQ=255) cAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATAAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAAGTACGCCGCGCCGATCACCAGCACAaccaa > 2:293353/1‑150 (MQ=255) gATCAGCGCCAGTAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCaa < 1:167945/150‑1 (MQ=255) aTCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAAc > 1:457934/1‑150 (MQ=255) tCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaac < 2:1014593/143‑1 (MQ=255) tCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCacaac > 1:1014593/1‑143 (MQ=255) ccagccGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGtcctcc > 1:184130/1‑150 (MQ=255) agccGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCtt < 2:1184700/150‑1 (MQ=255) ccGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTgg < 2:184130/150‑1 (MQ=255) ggggCACTACTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAAc > 2:38052/1‑150 (MQ=255) ctgAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCtt < 1:802354/150‑1 (MQ=255) gAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGa > 2:495909/1‑150 (MQ=255) aaaaTGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACcaca > 1:1088105/1‑109 (MQ=255) aaaaTGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACcaca < 2:1088105/109‑1 (MQ=255) catgagccATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGc < 1:67122/144‑1 (MQ=255) aaaTGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCCACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGc < 2:350790/150‑1 (MQ=255) cATCAACAGCACCTAGGCCGCCACCATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCa < 2:478302/150‑1 (MQ=255) aTCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCAcc > 2:92053/1‑84 (MQ=255) aTCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCAcc < 1:92053/84‑1 (MQ=255) aaCAGCACATAGGCCGCCACGATGCCGAGAAGCTCCAAGGGCGCCAGGCGCCAGTGCGTGAGGTACGTCGCGATGAACACCAGCACCAGCAACCACAACAACAACAGAATCGGTCCTCCTAGTAACCAAAGGCTTGACCCATCCcacaca < 2:229273/150‑1 (MQ=255) gATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGa > 1:244693/1‑150 (MQ=255) cGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGc > 2:313802/1‑150 (MQ=255) tCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTc > 1:244451/1‑150 (MQ=255) cAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGaa > 1:1177503/1‑150 (MQ=255) accggCGCCTGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGa < 1:1157402/98‑1 (MQ=255) accggCGCCTGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGa > 2:1157402/4‑101 (MQ=255) gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGc > 1:514817/1‑101 (MQ=255) gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGc < 2:514817/101‑1 (MQ=255) gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAg < 1:313802/150‑1 (MQ=255) gggCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAg > 1:167038/1‑150 (MQ=255) cAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAg > 2:57673/1‑150 (MQ=255) | GTTTGCGTCGCCACTGCGGGAGCGCAACCAGGGCGATCTTCAGGGCCAGCACCACCCAGATCAGCGCCAGCAGCCAGCCGGGGGCACTGCTGAAAATGCCCATCAACAGCACATAGGCCGCCACGATGCCGAGAATCTGCAAGGGCGCCAGGCGCCGGTGCGTGAGGTACGCCGCGCCGATCACCAGCACCACCAACCACAACAGCAACATAATCGGTCCTCCTTGGAACCAAGGGCTTGAGCCATTCCAGACAGCTTAGACGCAATGCTGTGAACGGCATGTTCAGAACAGTGACAG > NC_002947/2136583‑2136880 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |