Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,457,105 | 0 | A | G | 24.2% | 227.5 / 23.1 | 91 | I441V (ATC→GTC) | mpl | UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase |
Reads supporting (aligned to +/- strand): ref base A (20/49); new base G (22/0); total (42/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.31e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCA > NC_000913/4457058‑4457154 | aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTc > 1:4352674/1‑50 (MQ=38) aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGAt < 1:2285140/49‑1 (MQ=255) aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTc > 1:7725884/1‑50 (MQ=38) gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTcc > 1:5270486/1‑50 (MQ=38) gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATcc < 1:5825260/50‑1 (MQ=255) ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCa > 1:7811905/1‑50 (MQ=255) ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCa < 1:274022/50‑1 (MQ=255) cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAt < 1:8263574/50‑1 (MQ=255) tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTccct > 1:3284028/1‑47 (MQ=38) tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGATccct > 1:3274279/1‑47 (MQ=38) tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTc > 1:6035369/1‑50 (MQ=37) gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCa < 1:4645832/49‑1 (MQ=255) gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCa < 1:7316161/49‑1 (MQ=255) gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAg < 1:5388488/50‑1 (MQ=255) gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAg < 1:2539369/50‑1 (MQ=255) gCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTccct > 1:8519593/1‑45 (MQ=37) cGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAg > 1:1565679/1‑49 (MQ=38) cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGa < 1:7029607/50‑1 (MQ=255) aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGaaa > 1:6152395/1‑50 (MQ=38) aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGCGGGGTCCACCGGaa > 1:728583/1‑49 (MQ=25) aCCATATTCTGGTGATGAGCAACGGCGGTCTTGGTGGGGACACTCAGaaa > 1:5970668/1‑50 (MQ=25) ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGaaaa > 1:4267101/1‑49 (MQ=255) ccATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAc < 1:2717502/50‑1 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct < 1:7903819/50‑1 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct < 1:5703070/50‑1 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct < 1:6110674/50‑1 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAct < 1:5812339/50‑1 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGATCCATCAGAAAct > 1:5229695/1‑50 (MQ=255) cATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAGaaaat > 1:830956/1‑48 (MQ=25) atatTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAActg > 1:8424894/1‑50 (MQ=38) atatTCTGGTGATGAGCAACGGCGGTTTTGGTTGGGTTCATCCAAAAccc > 1:2990218/1‑48 (MQ=21) atatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAAatg > 1:277372/1‑50 (MQ=255) tatTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAatgc > 1:6488145/1‑50 (MQ=38) tatTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAatgc > 1:5728598/1‑50 (MQ=38) tatTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgc < 1:3757538/50‑1 (MQ=255) atTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgct > 1:2192992/1‑50 (MQ=255) ttCTGGTGATGAGCAACGGCGGTTTTTGTTGGGTGCCTCAGAAAacgctg > 1:3323618/1‑50 (MQ=16) ttCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAActgctg > 1:4396435/1‑50 (MQ=38) ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgct < 1:1861710/49‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGGTCCATTAGAAACTGCTGGGTgg > 1:6981410/1‑50 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGAt < 1:3823472/48‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATgg < 1:1721866/50‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATgg < 1:8281117/50‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg < 1:3499182/49‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg < 1:775352/49‑1 (MQ=255) tgatgaGCAACGGCGGTTTTTGTGGGGTCCATCAGAAACTGCTGGATggg > 1:518174/1‑50 (MQ=255) tgatgaGCAACGGCGGTTTTGGTGGGATCCCTCAGAAACTGCTGGATggg > 1:4120201/1‑50 (MQ=255) tgatgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATggg < 1:3211212/50‑1 (MQ=255) tgatgaGCAACGGCGGTTTTGGGGGGGACCCTCAGAAACTGCtgggtggg > 1:3459456/1‑50 (MQ=255) atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGc < 1:912846/49‑1 (MQ=255) atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGc < 1:4156113/49‑1 (MQ=255) atgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCt < 1:190354/50‑1 (MQ=255) tgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTg < 1:764526/50‑1 (MQ=255) tgaGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTg < 1:2578403/50‑1 (MQ=255) aGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTgg < 1:7167091/49‑1 (MQ=255) aGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGc < 1:686045/50‑1 (MQ=255) aGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGc < 1:7479794/50‑1 (MQ=255) aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaa < 1:6342093/50‑1 (MQ=255) aaCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaa < 1:1521963/50‑1 (MQ=255) cggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaaga < 1:3014533/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa < 1:6628527/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa < 1:6311320/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa < 1:3346391/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa < 1:1544165/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa < 1:4608020/50‑1 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaag > 1:5762715/1‑48 (MQ=255) gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag < 1:1529397/50‑1 (MQ=255) ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGc > 1:1047528/1‑50 (MQ=255) ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGc > 1:6584283/1‑50 (MQ=255) ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGc > 1:5015432/1‑50 (MQ=255) gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg < 1:3379050/50‑1 (MQ=255) gTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg < 1:2598080/50‑1 (MQ=255) ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCg > 1:3148072/1‑49 (MQ=255) tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg > 1:1208307/1‑49 (MQ=255) tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg > 1:1455354/1‑49 (MQ=255) tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg > 1:4007787/1‑49 (MQ=255) ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCgg > 1:8615945/1‑48 (MQ=255) ttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGa < 1:380298/49‑1 (MQ=255) tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa < 1:28172/49‑1 (MQ=255) tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAg > 1:6471087/1‑50 (MQ=255) tggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAg > 1:2104533/1‑50 (MQ=255) gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc > 1:1770118/1‑50 (MQ=255) gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc > 1:3167460/1‑50 (MQ=255) tggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcg > 1:6531669/1‑50 (MQ=255) tggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcg < 1:1486101/50‑1 (MQ=255) tggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcg > 1:5991351/1‑50 (MQ=255) ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg < 1:4153337/50‑1 (MQ=255) ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg < 1:8510608/50‑1 (MQ=255) ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg < 1:4852377/50‑1 (MQ=255) ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc < 1:7914437/49‑1 (MQ=255) ggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc < 1:3790495/49‑1 (MQ=255) aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg > 1:5944667/1‑48 (MQ=255) aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg > 1:4026328/1‑48 (MQ=255) aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgcg > 1:8397498/1‑48 (MQ=255) aTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCa < 1:8536382/50‑1 (MQ=255) | AGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCA > NC_000913/4457058‑4457154 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |