Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 987,014 | 0 | A | C | 29.2% | 45.2 / 14.0 | 24 | intergenic (‑32/+571) | ompF/asnS | outer membrane porin F/asparagine‑‑tRNA ligase |
Reads supporting (aligned to +/- strand): ref base A (17/0); new base C (0/7); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.89e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGT > NC_000913/986950‑987078 | gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:3035639/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:996640/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:3979444/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:3347543/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:215795/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:2100459/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:1970384/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:1801520/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:1610156/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:1170050/1‑65 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:1324085/1‑65 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 1:3217454/65‑1 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 1:1893149/65‑1 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 1:337775/65‑1 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 1:1288771/1‑65 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 1:344334/1‑65 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 1:484222/1‑65 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCa > 1:2292161/1‑65 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 1:2474441/1‑65 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 1:3411354/1‑65 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGc > 1:1372917/1‑65 (MQ=255) aaaaTTTTGCCCTTCATCATTTTTTTTTCCCCTCTTTTTTTTTTTTATCCCCCCTGCCACTGCCg < 1:2214021/63‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 1:3524726/1‑65 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 1:865334/1‑65 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 1:349850/1‑65 (MQ=255) cgcTTCATTTTTATTTATTACCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGt < 1:34891/65‑1 (MQ=255) cgcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATCCCCCCTGCCACTGCCGTCAATAAGt < 1:1676552/65‑1 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa > 1:836740/1‑65 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:3374614/1‑65 (MQ=255) cccTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg < 1:2784433/65‑1 (MQ=255) cccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 1:3737350/1‑65 (MQ=255) cccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 1:2277365/1‑65 (MQ=255) ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:1108057/1‑65 (MQ=255) ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:2216484/1‑65 (MQ=255) ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:3614178/1‑65 (MQ=255) tCATGGTTTTTTTTAGGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 1:3850968/65‑1 (MQ=255) ttttttttttATCACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:2716303/63‑1 (MQ=255) tgTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:3999771/64‑1 (MQ=255) gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:3712255/63‑1 (MQ=255) gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:1891469/63‑1 (MQ=255) gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:111989/63‑1 (MQ=255) gtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:3234100/63‑1 (MQ=255) gtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:146488/63‑1 (MQ=255) gtttttttttAGGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:3941153/63‑1 (MQ=255) ggTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTGACGGAACTAt < 1:3027592/65‑1 (MQ=255) ttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:393357/65‑1 (MQ=255) ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat > 1:875881/1‑65 (MQ=255) ttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 1:2533364/65‑1 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 1:995196/1‑65 (MQ=255) tATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag < 1:3496678/65‑1 (MQ=255) acCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt > 1:2070541/1‑65 (MQ=255) | GACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGT > NC_000913/986950‑987078 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |