| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,335,773 | 0 | G | A | 56.2% | 25.1 / 62.3 | 48 | A241V (GCG→GTG) | murA | UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (11/10); new base A (20/7); total (31/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGATAGAC > NC_000913/3335725‑3335822 | tCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATcgccgccgc > 1:1587012/1‑50 (MQ=255) ggCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCcaccac > 1:3264798/1‑50 (MQ=255) gCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCcaccac < 1:759210/49‑1 (MQ=255) gCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCcaccac < 1:2068909/49‑1 (MQ=255) cTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCcaccac > 1:1070739/1‑48 (MQ=255) cTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGccacca < 1:2521464/50‑1 (MQ=255) tGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAg > 1:290759/1‑50 (MQ=255) cgcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAgg > 1:2265410/1‑49 (MQ=255) cgcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAgg > 1:2526754/1‑49 (MQ=255) cgcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAgg > 1:2946353/1‑49 (MQ=255) cgcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGa < 1:3921080/50‑1 (MQ=255) gcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGaa < 1:2717554/50‑1 (MQ=255) gcgTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGa > 1:153438/1‑49 (MQ=255) gTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAg > 1:2536089/1‑50 (MQ=255) ttACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAg > 1:658918/1‑49 (MQ=255) ttACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGt < 1:1195553/50‑1 (MQ=255) aCGGCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTAc > 1:2822568/1‑50 (MQ=255) ggCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCg < 1:1773542/50‑1 (MQ=255) ggCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCg > 1:4316339/1‑50 (MQ=255) ggCAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCg > 1:482920/1‑50 (MQ=255) cAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGt > 1:744491/1‑50 (MQ=255) cAGATAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGt > 1:3561681/1‑50 (MQ=255) gATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGttt > 1:752980/1‑50 (MQ=255) gATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGttt > 1:654397/1‑50 (MQ=255) tAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTc > 1:2737501/1‑49 (MQ=255) tAATTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTc > 1:2596462/1‑49 (MQ=255) aaTTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGa > 1:3614821/1‑50 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGa < 1:2811576/49‑1 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGa > 1:1605482/1‑49 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGAt < 1:901545/50‑1 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGa > 1:2607102/1‑49 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGa < 1:1829249/49‑1 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGa < 1:1913798/49‑1 (MQ=255) aTTTTGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGAt < 1:4254949/50‑1 (MQ=255) tttGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATAc > 1:771805/1‑50 (MQ=255) tGCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGa > 1:4307998/1‑50 (MQ=255) gCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGa < 1:393263/49‑1 (MQ=255) gCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGAt > 1:2922754/1‑50 (MQ=255) gCCGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGAt < 1:4279013/50‑1 (MQ=255) ccGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATc > 1:1361115/1‑50 (MQ=255) ccGCGAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATc > 1:2292494/1‑50 (MQ=255) cgAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCgg < 1:678540/50‑1 (MQ=255) cgAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCgg > 1:3482444/1‑50 (MQ=255) cgAGAAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCgg > 1:4139904/1‑50 (MQ=255) gagaAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCgg > 1:3554641/1‑49 (MQ=255) agaAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCgg > 1:242112/1‑48 (MQ=255) agaAATCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCgg > 1:324399/1‑48 (MQ=255) aaaTCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGa < 1:1112495/50‑1 (MQ=255) aaTCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGaa < 1:1997660/50‑1 (MQ=255) aaTCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGaa < 1:3692373/50‑1 (MQ=255) aaTCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGaa > 1:3580167/1‑50 (MQ=255) aaTCGCCGCCACCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGaa > 1:2656555/1‑50 (MQ=255) aTCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAAc < 1:3870440/50‑1 (MQ=255) cgccgccacCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAAcgc > 1:658792/1‑50 (MQ=255) ccgccACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGa > 1:714774/1‑47 (MQ=255) gccACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGATAGac < 1:314904/50‑1 (MQ=255) | TCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGATAGAC > NC_000913/3335725‑3335822 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |