Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,510,928 | A→T | 52.0% | L177Q (CTG→CAG) | fecE ← | ferric citrate ABC transporter ATP binding subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,510,928 | 0 | A | T | 52.0% | 0.3 / 30.8 | 25 | L177Q (CTG→CAG) | fecE | ferric citrate ABC transporter ATP binding subunit |
Reads supporting (aligned to +/- strand): ref base A (7/5); new base T (7/6); total (14/11) | |||||||||||
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.40e-01 |
AATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCT > NC_000913/4510687‑4511177 | aaTCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATACTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTCTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGCCCACCTGGTGATTGATATCAAgtaagg > 1:128872/1‑267 (MQ=255) gCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAAATAGGTGGCTGGCCCATCAAGTaaaaaa > 2:60587/1‑268 (MQ=255) aTACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa < 1:60587/272‑1 (MQ=255) cGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGAGCGTATTCTGGGcc > 1:185575/1‑272 (MQ=255) tGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTgg > 1:133370/1‑192 (MQ=255) ggggTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGc < 2:133370/192‑1 (MQ=255) aCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGt > 2:215377/1‑206 (MQ=255) ccTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa < 1:215377/206‑1 (MQ=255) gcgcCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCCACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGAtgattg < 1:13709/272‑1 (MQ=255) cATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCt > 2:170481/1‑161 (MQ=255) aTAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtc < 1:170481/161‑1 (MQ=255) cATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGt > 2:28173/1‑178 (MQ=255) aTGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGGCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa < 1:28173/178‑1 (MQ=255) tGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGGTAGGTGGTTGGCCCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCCGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCCCGGGTAACCGACGAACGGGAAGATGATTGATCCGGGTATGGGTCATGGCGACATTAACTCGTGc > 1:196029/1‑272 (MQ=255) cGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCt > 1:15335/1‑126 (MQ=255) gCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCtc < 2:15335/126‑1 (MQ=255) cGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggcgct > 1:103867/1‑181 (MQ=255) ggCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCgcgctggcgctg < 2:103867/181‑1 (MQ=255) gATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGAGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGCCATTAACACGTGCATTGTCTTAAGCGGAAAGACGc > 1:174028/1‑272 (MQ=255) gccacacagACCCCCGTTTTCCCCTGGGTCCGGAGTGCGCCCATCAACCGCATCAGGTCCACCTGGTGGTTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGGATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAg < 2:196029/264‑1 (MQ=255) agcccagcccacaccGTTTTCCCCTGAGTCCGGGGTTCGCCCTTCAACCGCAGCTGGTCCACCTGGTGATGGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAg < 2:185575/261‑1 (MQ=255) agcaCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGAAAg > 1:57470/1‑272 (MQ=255) agcaCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAACGCTCGGTTAACCGACGAACGGCAAGCTGATTGATCCGGGTCTGGTTCATGGCGGCCTGAACTCGTGCATTGGCTTCAGCGGAGAGACGCCCCCAGATGGACaa > 2:221334/1‑271 (MQ=255) gTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGAAGCCCCCAGAGTGACAGCCAGGGATTACGAcc > 1:231036/1‑272 (MQ=255) gAGTTCGCCCATCAACCGCATCATGTCCGCCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCt < 1:221334/272‑1 (MQ=255) | AATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCT > NC_000913/4510687‑4511177 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |