Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,368,569 | 0 | A | C | 21.9% | 199.1 / 13.9 | 66 | L177F (TTA→TTC) | arnA | fused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (45/5); new base C (0/14); total (45/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.43e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCATCATAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGAAAA > NC_000913/2368520‑2368616 | tGCATCATAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTa > 1:4138812/1‑50 (MQ=255) catcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTa < 1:4846237/48‑1 (MQ=255) catcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTa < 1:6158029/48‑1 (MQ=255) catcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTAc < 1:4967305/49‑1 (MQ=255) catcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTAc < 1:6389085/49‑1 (MQ=255) atcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTGccc < 1:2113584/50‑1 (MQ=255) atcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTAccc > 1:2177330/1‑50 (MQ=255) atcatAAATTGTGCCATGCCGCGCGCCAGCCACTGGAACAGACATTAccc > 1:694388/1‑50 (MQ=255) tcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCg > 1:3309820/1‑50 (MQ=255) tcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCg > 1:5391548/1‑50 (MQ=255) tcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCg < 1:4420915/50‑1 (MQ=255) tcatAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCg < 1:3848371/50‑1 (MQ=255) catAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGc > 1:2221744/1‑50 (MQ=255) tAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCa > 1:2291952/1‑50 (MQ=255) tAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCa > 1:7066730/1‑50 (MQ=255) tAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCa < 1:4305709/50‑1 (MQ=255) aaaTTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCAt > 1:6443909/1‑50 (MQ=255) aaTTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCAt > 1:1787551/1‑49 (MQ=255) aaTTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCAt > 1:6305982/1‑49 (MQ=255) aTTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCAtt > 1:4067234/1‑49 (MQ=255) aTTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCAtt > 1:7178669/1‑49 (MQ=255) ttGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTa > 1:31834/1‑49 (MQ=255) ttGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTa > 1:1231475/1‑49 (MQ=255) tgtgCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTaaa > 1:5307317/1‑50 (MQ=255) tgtgCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTaaa > 1:246873/1‑50 (MQ=255) gtgCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTaaa > 1:7272969/1‑49 (MQ=255) tgCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAaca > 1:2940116/1‑50 (MQ=255) aTGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAACACgg > 1:1025107/1‑49 (MQ=255) tGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGt > 1:2034045/1‑49 (MQ=255) gCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTaa > 1:4409862/1‑50 (MQ=255) cgcgcgCCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAata > 1:4786199/1‑50 (MQ=255) gcgcgcCGGCTCGTGAAACAGACATTCCCCGCCATTAAACACGGTAatat < 1:4934194/50‑1 (MQ=25) gcgcgcCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAatat > 1:2253430/1‑50 (MQ=255) gcgcccCACCTACTGAAACAGACATTCCCCGCCATTAAACACGGTAatat < 1:1466487/50‑1 (MQ=255) cgcgcCGGCTACGGAAACAGACTTTCCCCGCCATTAAACACGGTAATAtt < 1:5040256/50‑1 (MQ=255) gcgcCAGTTACTGAAACAGAAATTCCCCGCCATTAAACACGGTAATATTc < 1:5294880/50‑1 (MQ=37) gcCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTg > 1:7052457/1‑50 (MQ=255) cAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTgg > 1:2422493/1‑49 (MQ=255) cAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTgg > 1:252625/1‑49 (MQ=255) gCTACTGAAACAGCCATTCCCCGCCATTAAACACGGTAATATTCTGGaaa < 1:4765840/50‑1 (MQ=255) gCTACGGAAACAGACTTTCCCCGCCATTAAACACGGTAATATTCTGGaaa < 1:510973/50‑1 (MQ=255) gCTACGGAAACAAACATTCCCCGCCATTAAACACGGTAATATTCTGGaaa < 1:1675399/50‑1 (MQ=255) cTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAAt > 1:3135041/1‑50 (MQ=255) cTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAAt > 1:6456723/1‑50 (MQ=255) cTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAAt > 1:4278223/1‑50 (MQ=255) tACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATc > 1:93392/1‑50 (MQ=255) aCTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCg > 1:1024118/1‑50 (MQ=255) aCTGAAAAAAACATTCCCCGCCATTAAACACGGTAATATTCTGGAAAt < 1:4461195/48‑1 (MQ=255) cTGGAACAGACATTACCCGCCATTAAACACGGTACTATTCTGGAAATCGc > 1:3930312/1‑50 (MQ=255) tGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGcc > 1:3842671/1‑50 (MQ=255) tGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGcc > 1:5524779/1‑50 (MQ=255) ggAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGcc > 1:2817918/1‑49 (MQ=255) aaaaCAGACATACCCCGCCATTAAACACGGTAATATTCTGGAAATCGccc < 1:1492900/48‑1 (MQ=255) gAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGcc > 1:6633255/1‑48 (MQ=255) aaCAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAg > 1:3425635/1‑50 (MQ=255) aaaGACATGGCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgc < 1:7393961/48‑1 (MQ=255) aCAAACATTCCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgc < 1:5243551/50‑1 (MQ=255) cAGACATTCCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgcg < 1:6300523/50‑1 (MQ=255) cAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgcg > 1:501908/1‑50 (MQ=255) aGACATTCCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgcgc < 1:265769/50‑1 (MQ=255) gACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAgcgc > 1:904201/1‑49 (MQ=255) aCATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGcgcg > 1:4714844/1‑49 (MQ=255) aCATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGcgcg > 1:4435131/1‑49 (MQ=255) cATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGa > 1:735676/1‑49 (MQ=255) aTTCCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGaaa < 1:6954076/50‑1 (MQ=255) aTTCCCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGaaa < 1:5018226/50‑1 (MQ=255) aTTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGaa > 1:4269565/1‑49 (MQ=255) ttACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGaaaa > 1:2518691/1‑50 (MQ=255) ttACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGaaaa > 1:1157334/1‑50 (MQ=255) | TGCATCATAAATTGTGCCATGCCGCGCGCCAGCTACTGGAACAGACATTACCCGCCATTAAACACGGTAATATTCTGGAAATCGCCCAGCGCGAAAA > NC_000913/2368520‑2368616 |
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: ‑ |