| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 657,797 | C→A | 8.7% | D28E (GAC→GAA) | cspE → | DNA‑binding transcriptional repressor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 657,797 | 0 | C | A | 8.7% | 292.5 / 7.5 | 91 | D28E (GAC→GAA) | cspE | DNA‑binding transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base C (46/37); new base A (6/2); total (52/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.59e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01 | |||||||||||
AACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCG > AP009048/657732‑657861 | aaCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt > 1:4278237/1‑70 (MQ=255) cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg < 1:3431917/71‑1 (MQ=255) cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg < 1:4821509/71‑1 (MQ=255) cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg < 1:4435674/71‑1 (MQ=255) gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg > 1:3369993/1‑70 (MQ=255) tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa > 1:2956402/1‑70 (MQ=255) tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa > 1:3956411/1‑70 (MQ=255) aaGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac < 1:4326880/70‑1 (MQ=255) gTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACt > 1:4092656/1‑71 (MQ=255) gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct < 1:1870735/71‑1 (MQ=255) gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct < 1:7325972/71‑1 (MQ=255) gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct < 1:1673916/71‑1 (MQ=25) gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct < 1:6192711/71‑1 (MQ=25) tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTacac > 1:3305081/1‑64 (MQ=255) tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg > 1:3708827/1‑71 (MQ=255) tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg > 1:2044733/1‑71 (MQ=255) tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg > 1:2713369/1‑71 (MQ=255) aGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc < 1:1736395/71‑1 (MQ=255) gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:5668175/1‑70 (MQ=255) gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:7035422/1‑70 (MQ=255) gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:7063905/1‑70 (MQ=255) tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt < 1:7669629/49‑1 (MQ=255) tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt < 1:4037709/49‑1 (MQ=255) tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:1791814/1‑69 (MQ=255) tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa < 1:2587352/71‑1 (MQ=255) ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:5165310/1‑68 (MQ=255) cAAAGTATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc < 1:5409624/60‑1 (MQ=255) cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTc > 1:4542572/1‑49 (MQ=255) cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc < 1:2329395/60‑1 (MQ=255) aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc > 1:7511701/1‑71 (MQ=255) aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg > 1:6676690/1‑47 (MQ=255) aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa > 1:6912680/1‑49 (MQ=255) ggATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac < 1:5762591/49‑1 (MQ=255) gATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt < 1:5627647/71‑1 (MQ=255) aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc < 1:6614943/61‑1 (MQ=255) tCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt > 1:3941755/1‑71 (MQ=255) cGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt < 1:7902918/70‑1 (MQ=255) ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct > 1:2361193/1‑49 (MQ=255) ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt < 1:3558143/66‑1 (MQ=255) gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg < 1:2895236/39‑1 (MQ=255) gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt < 1:371197/71‑1 (MQ=255) tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt > 1:4961038/1‑69 (MQ=255) ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa > 1:805729/1‑56 (MQ=255) ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa > 1:6116692/1‑56 (MQ=255) tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt > 1:1441973/1‑68 (MQ=255) tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa > 1:5787280/1‑70 (MQ=255) tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa > 1:5446529/1‑70 (MQ=255) cATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa > 1:2539389/1‑71 (MQ=255) aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca > 1:7744232/1‑38 (MQ=255) aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa < 1:3504256/70‑1 (MQ=255) aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc > 1:1308675/1‑71 (MQ=255) ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc > 1:1283595/1‑55 (MQ=255) ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct > 1:5507068/1‑71 (MQ=255) aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca > 1:4893129/1‑35 (MQ=255) aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc < 1:6520011/70‑1 (MQ=255) aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc < 1:2718746/70‑1 (MQ=255) cTCCTGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc > 1:1057902/1‑67 (MQ=255) tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc < 1:1127338/39‑1 (MQ=255) cGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCa > 1:2385637/1‑41 (MQ=255) ggAAGACGGCAGCAAAGACGTGTTCTTACACTTCTCTGCAATCCAGACTAATGGTTTTaa < 1:364742/60‑1 (MQ=255) ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa < 1:2888888/60‑1 (MQ=255) ggAAGACGGAAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg < 1:7858808/71‑1 (MQ=255) gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc > 1:2234123/1‑42 (MQ=255) gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa < 1:3585061/71‑1 (MQ=255) aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGa > 1:2249522/1‑45 (MQ=255) aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg > 1:3603404/1‑69 (MQ=255) aGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg > 1:6733109/1‑68 (MQ=255) gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa < 1:1332503/68‑1 (MQ=255) gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa < 1:6056808/68‑1 (MQ=255) gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa < 1:6065605/68‑1 (MQ=255) gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAgg > 1:1862221/1‑71 (MQ=255) aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt < 1:6885882/71‑1 (MQ=255) cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc > 1:47654/1‑42 (MQ=255) cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg > 1:5650000/1‑68 (MQ=255) cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc < 1:5249713/71‑1 (MQ=255) ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt > 1:4800274/1‑69 (MQ=255) ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCa < 1:7136193/71‑1 (MQ=255) cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTa < 1:6929864/50‑1 (MQ=255) cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc < 1:4612432/71‑1 (MQ=255) cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc < 1:5720084/71‑1 (MQ=255) agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa > 1:3387973/1‑52 (MQ=255) aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt > 1:2666809/1‑42 (MQ=255) aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt > 1:4568527/1‑42 (MQ=255) aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt > 1:5670939/1‑42 (MQ=255) aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct < 1:317109/49‑1 (MQ=255) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg > 1:497230/1‑66 (MQ=25) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg > 1:487866/1‑66 (MQ=25) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg > 1:7816748/1‑66 (MQ=25) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg > 1:5940894/1‑66 (MQ=25) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGcgcg > 1:1229900/1‑68 (MQ=25) aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTCCTGAAGGTCAGcgcg > 1:6715600/1‑68 (MQ=25) | AACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCG > AP009048/657732‑657861 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |