| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 310,934 | 0 | A | C | 42.9% | 28.6 / 22.3 | 28 | F135V (TTC→GTC) | ecpR | DNA‑binding transcriptional dual regulator MatA |
| Reads supporting (aligned to +/- strand): ref base A (13/3); new base C (0/12); total (13/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACG > NC_000913/310871‑310998 | tCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAAc < 1:1575288/65‑1 (MQ=255) cATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAAcc > 1:2197760/1‑65 (MQ=255) aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCAtt > 1:1279461/1‑65 (MQ=255) tAATTTCCATTCCCCGGTCCGTTATTTTATCTTTCTTAATATTTTTTCTCAGGACCCTGCCATTa < 1:2894726/65‑1 (MQ=255) tttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATc > 1:3355930/1‑65 (MQ=255) tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATcac < 1:3099610/65‑1 (MQ=255) tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATcac > 1:2445340/1‑65 (MQ=255) tccattcccCGTTCCGTGATTTTTTCTTTTTTAATATCTTTTATCAGGACCCTGCCATTAATcac < 1:3916914/65‑1 (MQ=255) ttccCGGTCCGTTATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa < 1:2446151/65‑1 (MQ=255) ttccCGGTCCGTTATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa < 1:2421820/65‑1 (MQ=255) ttccCGGTCCGTGATTTTATTTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa < 1:1627300/65‑1 (MQ=255) tccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGAc > 1:2851013/1‑65 (MQ=255) ggTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt < 1:3552066/65‑1 (MQ=255) ggTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGt < 1:2541520/65‑1 (MQ=255) tttATCTTTCTTAATATTTTTTCTCATGCACCTGCCATTAATCACATAGACCAGTTCTTTTCTaa < 1:2352688/65‑1 (MQ=255) tATCTTTTTTAATATCTTTCCCCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAata < 1:3307650/65‑1 (MQ=255) tATCTTTCTTAATATCTTTTCTCAGGACTCTGCCATTAATCACATAGACCAGTTCTTTTCTAata < 1:695345/65‑1 (MQ=255) aTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAatat > 1:1011524/1‑65 (MQ=255) tctttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc < 1:2667159/65‑1 (MQ=255) ttcttAATATCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg < 1:1994416/65‑1 (MQ=255) ttcttAATATCTTTTCACAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg < 1:1960825/65‑1 (MQ=255) cttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc < 1:856251/65‑1 (MQ=255) atatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAg < 1:184354/65‑1 (MQ=255) atatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAg < 1:379247/65‑1 (MQ=255) tttCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:3018417/63‑1 (MQ=255) tatTTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:1952006/65‑1 (MQ=255) tatCTTTTCTCGGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:728465/65‑1 (MQ=255) tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:1920425/65‑1 (MQ=255) tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:702997/65‑1 (MQ=255) tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:3707477/65‑1 (MQ=255) tatCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:438081/65‑1 (MQ=255) tatCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:3404060/65‑1 (MQ=255) tatCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:2240252/65‑1 (MQ=255) taccTTTCCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt < 1:2647100/62‑1 (MQ=255) atCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTc > 1:3649818/1‑65 (MQ=255) tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc > 1:3021409/1‑65 (MQ=255) tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc > 1:3131376/1‑65 (MQ=255) tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc > 1:3595917/1‑65 (MQ=255) tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc > 1:41902/1‑65 (MQ=255) tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc > 1:2383907/1‑65 (MQ=255) ttttCTCGGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg < 1:1069267/65‑1 (MQ=255) ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg < 1:657735/65‑1 (MQ=255) ttttCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg > 1:2791519/1‑65 (MQ=255) ttttCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg < 1:175388/65‑1 (MQ=255) tttcCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg < 1:3049536/65‑1 (MQ=255) tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc < 1:3150662/65‑1 (MQ=255) tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc < 1:177419/65‑1 (MQ=255) tttCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc > 1:1116793/1‑65 (MQ=255) tttCTCAGCACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc < 1:3395033/65‑1 (MQ=255) tttCCCAGGCCCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc < 1:1956938/65‑1 (MQ=255) tttCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc < 1:3256960/65‑1 (MQ=255) ttCACAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCg < 1:1507396/65‑1 (MQ=255) tctcAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt < 1:1624451/65‑1 (MQ=255) tctaAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt < 1:3330057/65‑1 (MQ=255) aaCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCAc > 1:970686/1‑65 (MQ=255) aCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATAccctg > 1:3074148/1‑62 (MQ=255) aCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACg > 1:695919/1‑65 (MQ=255) | TCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACG > NC_000913/310871‑310998 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |