| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 197,945 | T→C | 11.8% | S324P (TCA→CCA) | cmpC → | Bicarbonate transport ATP‑binding protein CmpC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 197,945 | 0 | T | C | 11.8% | 97.0 / 3.0 | 42 | S324P (TCA→CCA) | cmpC | Bicarbonate transport ATP‑binding protein CmpC |
| Reads supporting (aligned to +/- strand): ref base T (15/21); new base C (3/2); total (18/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GACATCCTATGGTGATTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAA > USA300TCH1516_ALE/197809‑198082 | gataTCGTATGGTGATTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATg < 2:370995/138‑1 (MQ=255) ggTGATTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCca < 1:100690/141‑1 (MQ=255) tGATTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATTGCGTTCATGACACGTCCcaca < 1:20839/141‑1 (MQ=255) tGATCTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCcaca < 1:10776/141‑1 (MQ=255) aTTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCACTTGTTTAATCCACCTGCATCTGATGCCTTTGTTGAACCGTCATTGTATAAGCATGCATCGCGTGGATTACGCGTCCCACaaa < 2:377594/141‑1 (MQ=255) ttAACAATTCAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGAACGCATCGCGTCCCTACCACGCCCCACAAATa > 2:189238/1‑141 (MQ=255) aaCAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAAc > 1:214486/1‑141 (MQ=255) aaCAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAAc > 1:61627/1‑141 (MQ=255) aCAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACa < 1:274695/141‑1 (MQ=255) aCAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCAGTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACa > 2:141586/1‑141 (MQ=255) ttAAGCCATCCGGCTATCAAGAAATTACGACATTGGCAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAAtt < 1:61380/141‑1 (MQ=255) atAAGCCATCCGTCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAAtt < 2:351472/140‑1 (MQ=255) taggCCAGCCGGCTATCAAGAAATTACGACATGGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAAttt < 1:305303/138‑1 (MQ=255) tAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAATGAGGCATCGCGTCCATGACACGTCCCACAAATAACAAAttt > 1:27637/1‑141 (MQ=255) tAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTCCATGACACGTCCCACAAATAACAAAttt > 1:27625/1‑141 (MQ=255) tAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGACATCGCGTTCATGACACGTCCCACAAATCACAAAttt > 2:119210/1‑141 (MQ=255) aGAACTTACGACATTGGGAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATcaca < 1:137570/141‑1 (MQ=255) accaCATTGGTAAAACAACATCATTTGTTTAATCCACCCGCATATGATGACTTTCTTGAACCGTCATTGTATAAGGAGGCAGCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATttatta < 2:8135/138‑1 (MQ=255) cATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATTACACGTCCCACAAATAACAATTTTATATGACCTATTATCACATTTATTAtttt < 1:256132/141‑1 (MQ=255) cATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTAtttt < 1:72947/141‑1 (MQ=255) cATTGGTAAAACAACATCATTTGTTTAATCCAACTGACTACGTTGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTAtttt < 1:206799/141‑1 (MQ=255) acaTCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTgg < 1:210731/141‑1 (MQ=255) acaTCATTTGTTTAATCCACCTGCATATGATGACTTTGGTAAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACCAATAACAAATTATAATTCCCTATTATCCCGCTTAATATTTTCTTGGGCATTTgg > 1:279580/1‑141 (MQ=255) tcatTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGaa < 1:165529/141‑1 (MQ=255) catTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCCCAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGaaa > 2:369469/1‑141 (MQ=255) gcgggcTTTGTTGAACCGTCATTGTATAAGGAGCCATCGTGTTCAAGACACGTCACACAGATAACAAGTTCATCTTACCTATTATCACAGTTATTATTTTCTTAGGCATGTGGGAAATGGTCAATATTATTGGGCATGAcc < 2:58895/136‑1 (MQ=255) gatgaCTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTCCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTAcc < 2:27625/141‑1 (MQ=255) gatgaCTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTCCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTAcc < 2:27637/141‑1 (MQ=255) atgaCTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCa < 1:178761/141‑1 (MQ=255) tttGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAAGAAATTTATATTACCGATTATCACATTTATTATTTTCTTAGGCATCTGGGAAATGGTCATTATTATTGGGCATTACCAACCTg < 2:262271/141‑1 (MQ=255) gAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGt < 2:154641/141‑1 (MQ=255) ttGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGTGCATGACCAACCTGTATTGTTACCGGGTc > 2:317963/1‑141 (MQ=255) ttGTATAAGGAGGCATCGCGTTCAGGACACGTCCCACTAATAACAAACTTATATTACGTATTATCACATTTATTATTGTCTTAGGCATTTCGGAAATGGTCATTATTATTGGGCATTACAAACCTGTATTGTTACCGGGTc > 2:95193/1‑141 (MQ=255) gtctaAGGAGGCATCGCGTTCATGACCCGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCt < 2:15586/138‑1 (MQ=255) ataAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGc > 1:291026/1‑141 (MQ=255) ggCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTa < 1:141586/141‑1 (MQ=255) gcgTTCATGACACGTCCCCCAAATAACAAATTTATATTACCGATTATAACGTTTATTTTTTTCTTAGGCATTTGGGAAATGGTCATTATTCTTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGGTATTGTAGGaaaa > 1:79295/1‑141 (MQ=255) gcgTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGGGCCGGTGAGTAGGATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGcaaa > 1:91842/1‑141 (MQ=255) gcgTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTTGGCATTACCAACCTGTATTGTTACCGGGTCGTGCTCTTGTAGGgaaa > 1:86136/1‑141 (MQ=255) gcgTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGaaaa > 2:234836/1‑141 (MQ=255) gcgTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTGATTATTTTCTTAGGCATTTGGGAAATGGGCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGaaaa > 1:67154/1‑141 (MQ=255) cgTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATTACATTTATTATTTTCTTCGGCATTTGGGATATGGTCATTATTATTGGGCATTACCACACTGTCTTGTTACCGGGTCCTGCTGTTGTAGgaagaa > 2:295448/1‑141 (MQ=255) | GACATCCTATGGTGATTTAACAATTAAGCCATCCGGCTATCAAGAAATTACGACATTGGTAAAACAACATCATTTGTTTAATCCACCTGCATATGATGACTTTGTTGAACCGTCATTGTATAAGGAGGCATCGCGTTCATGACACGTCCCACAAATAACAAATTTATATTACCTATTATCACATTTATTATTTTCTTAGGCATTTGGGAAATGGTCATTATTATTGGGCATTACCAACCTGTATTGTTACCGGGTCCTGCTCTTGTAGGAAAAA > USA300TCH1516_ALE/197809‑198082 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |