Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 3,459,849 | A→C | 10.8% | I79S (ATT→AGT) | bfr ← | bacterioferritin, iron storage and detoxification protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 3,459,849 | 0 | A | C | 10.8% | 97.9 / 8.2 | 37 | I79S (ATT→AGT) | bfr | bacterioferritin, iron storage and detoxification protein |
Reads supporting (aligned to +/- strand): ref base A (15/18); new base C (2/2); total (17/20) | |||||||||||
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.20e-01 |
TCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAA > CP009273/3459615‑3460114 | tcttcGCGGATCTGTGCTTGCAGATGATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGACATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACCCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCATCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTc < 1:676239/275‑1 (MQ=255) ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACCTCTTCCCCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAACCCTTcc > 1:30123/1‑274 (MQ=255) tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTGTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCATGTAAGTTTGGAAGACCTTCCAGAAAAAGaa > 1:149716/1‑274 (MQ=255) ttgacgttatgtttGAAGGCCCATCTTCTGAATCAAATCAAGTTCCGTTTCCAGCCAGTCGATATTGCCTTCTTCATCACGCAAAATTTGTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt < 2:149715/264‑1 (MQ=255) ttgaccaTAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt < 2:30123/269‑1 (MQ=255) tGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTg > 1:182113/1‑192 (MQ=255) tGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTg > 1:606949/1‑192 (MQ=255) aaTCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCANATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatcccat > 1:600881/1‑275 (MQ=255) tCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTGtgaa < 2:606940/192‑4 (MQ=255) tCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTGtgaa < 2:182110/192‑4 (MQ=255) gTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGaa > 1:206889/1‑212 (MQ=255) atatGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg < 2:206885/212‑1 (MQ=255) cttcttCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAg > 2:517075/1‑193 (MQ=255) ttcATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGaaaa < 1:517082/193‑1 (MQ=255) aaaaagCAAAATTTATATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGAcccc < 2:600872/270‑1 (MQ=255) tcaCGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg > 2:318841/1‑195 (MQ=255) gCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt < 1:318847/195‑1 (MQ=255) catATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTc > 2:628797/1‑127 (MQ=255) tCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGtt < 1:628806/127‑1 (MQ=255) gAACGCTATCGGCATTACCAAGTGTCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg > 1:492653/1‑152 (MQ=255) aaCGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAAt > 1:167313/1‑189 (MQ=255) aCGCTATCGGCATAACCAATTTCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTNAACATCTTNACCAATGTTCATATTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACCTTCGGGCCTGGAGAAAGTACTGATTGATTGCGAc < 1:464300/275‑1 (MQ=255) cGCTATCGGCATAACCAATTGCCGCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATTCGCTCAATATAACGATCGGCGTGTCTCATCTCATNAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGTGCATGGAGAAAGTACTGATTGATTGCGAc > 1:487930/1‑274 (MQ=255) gCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt < 2:492646/152‑1 (MQ=255) cTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGAt < 2:167310/189‑1 (MQ=255) ggCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATg < 2:135646/91‑1 (MQ=255) aTTGCCGCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTCCCTCAACATCTTCACCAATGTTCAgat > 2:138798/1‑86 (MQ=255) ccTCACGCAAATCCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACCTCTTCACCAATGTTCAGTTTGcc < 1:138799/88‑1 (MQ=255) tCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTttgt > 1:691632/1‑274 (MQ=255) gCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATATTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTgaga < 2:691623/275‑1 (MQ=255) aTTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGtt > 1:210993/1‑74 (MQ=255) ttcGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGcc < 2:210989/74‑1 (MQ=255) agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg > 1:249370/1‑78 (MQ=255) tttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt < 2:249365/78‑1 (MQ=255) ttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCt > 1:506451/1‑41 (MQ=255) tCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTaa < 2:506444/41‑1 (MQ=255) gtCACAAATGTTTAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTACAGAAATAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGGTTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACCGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGaa < 2:487923/272‑1 (MQ=255) | TCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAA > CP009273/3459615‑3460114 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |