| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,216,481 | (G)6→5 | coding (11/765 nt) | nfeF ← | NADPH‑dependent ferric chelate reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,216,476 | 0 | G | . | 93.5% | 113.6 / 2.1 | 31 | coding (16/765 nt) | nfeF | NADPH‑dependent ferric chelate reductase |
| Reads supporting (aligned to +/- strand): ref base G (1/1); new base . (12/17); total (13/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCT > NC_000913/3216344‑3216597 | aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCg > 1:15392/1‑133 (MQ=255)aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCCCAGATCATTGCGAAAGCGCTGCGGGTAGCg > 1:281241/1‑133 (MQ=255) gCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCgg < 2:15392/133‑1 (MQ=255) cACGCGATGTAAAACCGTCCAGCGCCTCACCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCgggg < 2:592065/133‑1 (MQ=255) gcgATGTAAAACCGTCCAGCGCCTCGTCTCCGAGGACAATGGGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCgggggtg < 2:281241/133‑3 (MQ=255) gcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCgggggtg < 2:584810/133‑3 (MQ=255) tGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGAGCAGATCAATGCTAACGCGCTGCGGATAGC‑GGGGGTGTTAt > 2:225372/1‑133 (MQ=255) aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCa < 2:321468/133‑1 (MQ=255) gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTa < 2:368090/133‑1 (MQ=255) tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACTCAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTAt > 1:610879/1‑133 (MQ=255) tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGG‑GTTATTCATTGTTAt > 2:223121/1‑133 (MQ=255) tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTAt > 2:387825/1‑133 (MQ=255) tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTAt > 2:311588/1‑133 (MQ=255) tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTAt > 1:550748/1‑133 (MQ=255) aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATcac < 2:114689/133‑3 (MQ=255) cTCGCCGCGGAGGACTATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCg < 1:265824/133‑1 (MQ=255) cgccgcCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTg < 2:27698/133‑1 (MQ=255) gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTaa < 1:311588/133‑1 (MQ=255) aGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt < 1:75499/133‑1 (MQ=255) tGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCta > 2:565045/1‑133 (MQ=255) cTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc < 1:387825/133‑1 (MQ=255) tGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa < 1:312241/133‑1 (MQ=255) aCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatga < 2:357728/133‑1 (MQ=255) tcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa > 1:421639/1‑133 (MQ=255) agtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATa < 2:421639/133‑1 (MQ=255) cgcAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTGGCAAATGAGAATGATTGTTAATCATGATAAATGCAAACGAt > 2:168908/1‑133 (MQ=255) cgcAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAt > 2:326286/1‑133 (MQ=255) cgcAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAt > 2:546509/1‑133 (MQ=255) aGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTg < 1:565045/133‑1 (MQ=255) aTCATTGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa < 1:546509/133‑1 (MQ=255) tCATTGCGCACGCGCGGCTGGTGGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATAGTTAATCATGATAAATGCAAGCGATTTGTAg < 1:223121/133‑1 (MQ=255) tGCGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACt > 1:550057/1‑133 (MQ=255) cGAACGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGa < 2:550748/133‑1 (MQ=255) cGAACGCGATGCGGGTAG‑GGGGGGTGTTATCCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGa < 1:225372/133‑1 (MQ=255) aCGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat < 1:591452/131‑1 (MQ=255) aCGCGCTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat > 2:591452/1‑131 (MQ=255) gcTGCGGGTAGC‑GGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCt < 1:326286/133‑1 (MQ=255) | AGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCT > NC_000913/3216344‑3216597 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |