| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 3,564,985 | T→C | 11.1% | F21L (TTT→CTT) | yihS → | predicted glucosamine isomerase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 3,564,985 | 0 | T | C | 11.1% | 163.6 / 12.4 | 63 | F21L (TTT→CTT) | yihS | predicted glucosamine isomerase |
| Reads supporting (aligned to +/- strand): ref base T (28/28); new base C (2/5); total (30/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.29e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01 | |||||||||||
ATAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA > AP009048/3564916‑3565051 | aTAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCtt > 1:4293772/1‑71 (MQ=255) tAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCtt < 1:4429998/70‑1 (MQ=255) aGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAtt > 1:823562/1‑71 (MQ=255) aGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAtt < 1:1961512/71‑1 (MQ=255) ggAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt < 1:1289491/71‑1 (MQ=255) aTGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTgg > 1:2683154/1‑71 (MQ=255) aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTa < 1:712962/70‑1 (MQ=255) aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTa > 1:180208/1‑70 (MQ=255) aaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa < 1:2079373/71‑1 (MQ=255) aTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa > 1:3071793/1‑70 (MQ=255) tttAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAtt < 1:3839123/70‑1 (MQ=255) ttAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTaa > 1:71471/1‑64 (MQ=255) tAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTgg > 1:1789263/1‑60 (MQ=255) aaCACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCg < 1:3259557/70‑1 (MQ=255) acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGCCCGCATCTTTg > 1:1939196/1‑51 (MQ=255) acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaaa > 1:3784561/1‑63 (MQ=255) acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGt < 1:1992138/70‑1 (MQ=255) acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGt < 1:2737687/70‑1 (MQ=255) cccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt > 1:2540093/1‑52 (MQ=255) ccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTc < 1:2476993/63‑1 (MQ=255) ccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGc > 1:2970742/1‑71 (MQ=255) ccTAAGCCACAACC‑TTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAg < 1:4409477/67‑1 (MQ=255) cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc < 1:4180512/71‑1 (MQ=255) tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt < 1:3723819/49‑1 (MQ=255) tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt > 1:501572/1‑49 (MQ=255) tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc > 1:2638092/1‑70 (MQ=255) aGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAgaga > 1:1660167/1‑64 (MQ=255) acaACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGt < 1:923471/47‑1 (MQ=255) acaACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa < 1:1872682/50‑1 (MQ=255) acaACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGt < 1:1425666/71‑1 (MQ=255) caACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtt < 1:1337454/71‑1 (MQ=255) aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAt < 1:2737976/51‑1 (MQ=255) aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTcc > 1:682267/1‑54 (MQ=255) aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt < 1:1679833/71‑1 (MQ=255) aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt < 1:4223309/71‑1 (MQ=255) aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtt < 1:3578276/70‑1 (MQ=255) aCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt > 1:1105966/1‑40 (MQ=255) aCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCTTAGTGCCGACTGGtttt < 1:762792/71‑1 (MQ=255) cGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTgg < 1:606733/71‑1 (MQ=255) cGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAg < 1:1634492/55‑1 (MQ=255) ttGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa < 1:4499125/52‑1 (MQ=255) ggCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg > 1:3151819/1‑71 (MQ=255) ggCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg > 1:4535636/1‑71 (MQ=255) ggCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg > 1:3454736/1‑71 (MQ=255) cTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGt < 1:1934764/47‑1 (MQ=255) tGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGt < 1:4593622/46‑1 (MQ=255) gAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAgg > 1:4248318/1‑71 (MQ=255) aCAGGAAACCGTCCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAg < 1:4395735/68‑1 (MQ=255) aCAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGc > 1:217913/1‑61 (MQ=255) cAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa > 1:2423005/1‑42 (MQ=255) aGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGc > 1:653636/1‑45 (MQ=255) ggAAACCGACAGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAg > 1:3008781/1‑65 (MQ=255) gAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaa < 1:1870827/71‑1 (MQ=255) aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGa > 1:2922118/1‑45 (MQ=255) aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGa > 1:917744/1‑45 (MQ=255) aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaaa < 1:4478676/71‑1 (MQ=255) aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaaa < 1:894624/71‑1 (MQ=255) aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa < 1:630018/69‑1 (MQ=255) aaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAg < 1:2633438/71‑1 (MQ=255) gaccgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGc > 1:2300993/1‑70 (MQ=255) accgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCa > 1:1724915/1‑70 (MQ=255) cgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAAt > 1:4220883/1‑71 (MQ=255) cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt > 1:449899/1‑42 (MQ=255) cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa < 1:4207656/71‑1 (MQ=255) tCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa > 1:1303636/1‑69 (MQ=255) | ATAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA > AP009048/3564916‑3565051 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |