| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,661,866 | G→T | 42.3% | A89D (GCC→GAC) | iscR ← | DNA‑binding transcriptional dual regulator IscR |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,661,866 | 0 | G | T | 42.3% | 1.7 / 13.2 | 19 | A89D (GCC→GAC) | iscR | DNA‑binding transcriptional dual regulator IscR |
| Reads supporting (aligned to +/- strand): ref base G (5/6); new base T (4/4); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01 | |||||||||||
AGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTAC > NC_000913/2661725‑2661983 | aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTAGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCACGCATTTATCACCGCCCTGACAGCCGCCTTTACCCGGACAACAGGTAGc > 1:257540/1‑143 (MQ=255) ctCCTTGTTACTCCCCAGTTCGTCTCAAGTACTCTTGTTCAGAATGCCGTTGATATGGTCGCGCATCTCACGCCTCCGCGCGTGAGCGTGTGTCAGGCATTTATCGCCGCCCTGTCAGCCGCCTTTACCCTGACAACGGGTgg < 2:317271/142‑1 (MQ=255) tattaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCg < 2:354095/143‑1 (MQ=255) gTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGgcg < 1:412114/143‑1 (MQ=255) tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTa > 1:272404/1‑143 (MQ=255) aaTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCATATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGGCGCTAATTACTTCGc < 1:392800/143‑1 (MQ=255) gttgttGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGGCGCTAATTACTTTGCCaa > 2:160300/1‑143 (MQ=255) gaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATCTATCGCCGCCCTCGCAGCCGCCTTCACCCTGACAACGGGTGGCATCTACAGATGCGTCAACGGCACTAATCACCTCGACAAGGGCGAtc > 2:413987/1‑142 (MQ=255) aaaaCCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGGCGCTAATTACTTCCCTACCGGCGAtgc > 1:158886/1‑143 (MQ=255) gggTTAGACGGTCGCTCAAATCACGTCACATCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCGTGACAACGGGTGGCATCTACATATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGctgctg < 2:75886/142‑1 (MQ=255) gTGAGACGGTCGCTCAAATCCCGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTGCAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGc < 1:218529/143‑1 (MQ=255) tGAGACGGGCGGTCAAACCAAGCACCATCGCGTG‑‑‑‑‑G‑GGCAGGCCTTCATCGCCGCCCTGGCAGCCGCCTTTACCCGGACAACGGGTGGCACCTACAGAGTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCa < 1:115864/143‑1 (MQ=255) gTCGCTCCAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGCCGATGCTGCTGGCATCTTTGc < 1:332170/143‑1 (MQ=255) tCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGCCAACGGCGCTACTTACTTCGCCACCAGCGATGCTTCTGGCATCTTTGcc > 2:321265/1‑143 (MQ=255) acaGCGCGTA‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCg > 1:154634/1‑143 (MQ=255) gttCAGGCATTTATCGCCGCACTGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAGCCACCGCCTGGTCCAc < 1:395898/141‑1 (MQ=255) aTTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTg < 1:267478/143‑1 (MQ=255) tCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAc < 1:389735/143‑1 (MQ=255) tGGCAGCCGCCTTTACCCTGACAACGGGTGTCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGCTACAAGACCAtt > 2:199726/1‑143 (MQ=255) gccgccTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTAc > 1:304001/1‑143 (MQ=255) | AGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTG‑‑‑‑‑G‑GTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTAC > NC_000913/2661725‑2661983 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |