Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,440,194 | 0 | C | A | 52.8% | 8.6 / 104.7 | 53 | G279V (GGT→GTT) | rpoA | RNA polymerase subunit alpha |
Reads supporting (aligned to +/- strand): ref base C (11/14); new base A (10/18); total (21/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAC > NC_000913/3440145‑3440240 | aaGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc < 1:1567258/50‑1 (MQ=255) gTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa < 1:287937/50‑1 (MQ=255) ttAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa > 1:3611034/1‑49 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata > 1:2823250/1‑50 (MQ=255) tAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata > 1:2608411/1‑50 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTGCCAGATCAACGatat < 1:1967220/49‑1 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGAtata < 1:871883/50‑1 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGAtata < 1:2718969/50‑1 (MQ=255) ggCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGAtata < 1:905397/50‑1 (MQ=255) gCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAg < 1:161869/50‑1 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTTCCAGATCACCGATATAGTg < 1:337524/50‑1 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg > 1:1811432/1‑50 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTg > 1:2075773/1‑50 (MQ=255) gTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAACGATATAGTg > 1:582786/1‑50 (MQ=255) ttttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTgg < 1:174231/50‑1 (MQ=255) ttAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa < 1:2087820/49‑1 (MQ=255) aaGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATa < 1:1010014/49‑1 (MQ=255) aGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGc < 1:549637/50‑1 (MQ=255) ggAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCt < 1:1622788/50‑1 (MQ=255) gAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt < 1:1237882/50‑1 (MQ=255) aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc < 1:3462377/48‑1 (MQ=255) aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt > 1:1264264/1‑50 (MQ=255) aaCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt < 1:3068865/50‑1 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGc < 1:3222561/48‑1 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt > 1:2623736/1‑50 (MQ=255) aaCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCtt > 1:319845/1‑50 (MQ=255) aCCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt > 1:3777422/1‑49 (MQ=255) aCCTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt < 1:3130946/50‑1 (MQ=255) ccTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTa < 1:3798079/50‑1 (MQ=255) cTCGGTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCttt > 1:4118840/1‑48 (MQ=255) cGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg > 1:1594272/1‑50 (MQ=255) ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg > 1:3806157/1‑49 (MQ=255) ggTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGc < 1:2599396/50‑1 (MQ=255) ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg < 1:2070844/49‑1 (MQ=255) ggTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGc < 1:1798874/50‑1 (MQ=255) gTACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAgg > 1:111819/1‑48 (MQ=255) tACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg > 1:1265523/1‑50 (MQ=255) tACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAg < 1:3656136/50‑1 (MQ=255) cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt < 1:2477123/50‑1 (MQ=255) cGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt < 1:741114/50‑1 (MQ=255) gCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa < 1:1021584/50‑1 (MQ=255) cTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg > 1:1371301/1‑50 (MQ=255) tACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa < 1:4104237/49‑1 (MQ=255) tACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa < 1:1484642/49‑1 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag < 1:4009395/50‑1 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag > 1:4188538/1‑50 (MQ=255) ccAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag < 1:2599796/50‑1 (MQ=255) cAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc > 1:1926557/1‑50 (MQ=255) aGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg > 1:4159656/1‑50 (MQ=255) gATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg < 1:2251308/50‑1 (MQ=255) aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa > 1:517691/1‑50 (MQ=255) aTCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa < 1:3446030/50‑1 (MQ=255) tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa > 1:3976257/1‑49 (MQ=255) tCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAc < 1:766457/50‑1 (MQ=255) | AAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAC > NC_000913/3440145‑3440240 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |