Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,355,999 | 0 | G | A | 33.8% | 170.6 / 144.9 | 136 | W425* (TGG→TGA) | gltB | glutamate synthase subunit GltB |
Reads supporting (aligned to +/- strand): ref base G (40/50); new base A (22/24); total (62/74) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.20e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 158 total aligned reads displayed.
GCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATG > NC_000913/3355952‑3356048 | gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggat > 1:2680228/1‑50 (MQ=255) gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggat < 1:1968069/50‑1 (MQ=255) gCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaat > 2:2360260/1‑50 (MQ=255) cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAACGAGtgaatg > 2:4851207/1‑50 (MQ=38) cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg < 2:4363964/50‑1 (MQ=255) cAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg > 2:2721418/1‑50 (MQ=255) aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg > 1:2912551/1‑49 (MQ=255) aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaatgg < 1:1255553/50‑1 (MQ=255) aGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtgaatgg < 1:757022/50‑1 (MQ=255) gAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg > 1:4004712/1‑48 (MQ=255) gAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg > 1:3361319/1‑48 (MQ=255) aaaCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGtggatg > 1:3511147/1‑47 (MQ=255) aaaCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGgag < 2:1971678/50‑1 (MQ=255) aCCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGgaga > 2:2645847/1‑49 (MQ=255) ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa < 2:546922/50‑1 (MQ=255) ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa < 1:2862259/50‑1 (MQ=255) ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaa > 2:690102/1‑49 (MQ=255) ccGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaa > 2:3976835/1‑49 (MQ=255) cGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaa < 2:3100200/50‑1 (MQ=255) cGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaa > 2:5048490/1‑49 (MQ=255) gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaaa > 1:3794605/1‑50 (MQ=255) gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGaaaaa > 2:2200456/1‑50 (MQ=255) gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaaaaa > 1:2778516/1‑50 (MQ=255) gATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGaaaaa < 1:5084613/50‑1 (MQ=255) aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc < 1:1559175/50‑1 (MQ=255) aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc < 1:4696304/50‑1 (MQ=255) aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc < 2:2910054/50‑1 (MQ=255) aTGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAAc < 2:4457782/50‑1 (MQ=255) aCGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGt < 1:921381/49‑1 (MQ=255) cGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTcc < 1:650945/50‑1 (MQ=255) cGATCTGAAAAACCGCCATGCATATAAAGAGAGAATGTAGAAAAACGTcc < 2:3674926/50‑1 (MQ=255) gATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccg > 2:789761/1‑50 (MQ=255) aTCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgc > 2:3696564/1‑50 (MQ=255) tCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc > 2:2337978/1‑50 (MQ=255) cTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc > 1:2536629/1‑49 (MQ=255) cTGAAAAGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTccgcc > 1:3730993/1‑49 (MQ=255) tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa > 2:4993290/1‑50 (MQ=255) tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa > 1:3139916/1‑50 (MQ=255) tGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa > 1:5195561/1‑50 (MQ=255) gAAAAGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGAc > 1:3644231/1‑50 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgccg < 2:1195790/47‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgcc < 2:2173833/46‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccgc < 2:3924431/45‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTccg < 1:4159062/44‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa < 1:3373958/48‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGa < 1:1583924/48‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACt < 1:316710/50‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACt < 2:1929681/50‑1 (MQ=255) aaaaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACt < 1:2320082/50‑1 (MQ=255) aaaGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTg < 1:1889455/50‑1 (MQ=255) aaaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTg > 1:393840/1‑50 (MQ=255) aaGCCGCCATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTgg < 1:3201071/50‑1 (MQ=255) aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTgg > 1:5060138/1‑49 (MQ=255) aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTgg > 1:661274/1‑49 (MQ=255) aGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGt < 1:1895701/50‑1 (MQ=255) gccgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTa < 1:2445497/50‑1 (MQ=255) ccgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAc < 2:3076305/50‑1 (MQ=255) cgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAcc > 1:1069122/1‑50 (MQ=255) cgccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTAcc > 2:4237562/1‑50 (MQ=255) gccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCg > 2:1425974/1‑50 (MQ=255) gccATCCATATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCg > 2:3461282/1‑50 (MQ=255) ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt > 1:5050838/1‑50 (MQ=255) ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt > 1:4281028/1‑50 (MQ=255) ccatccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGt > 2:2989734/1‑50 (MQ=255) catccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGtt < 1:1800883/50‑1 (MQ=255) catccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGtt > 2:717923/1‑50 (MQ=255) atccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGttt > 2:5171302/1‑50 (MQ=255) atccatATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGttt < 1:426343/50‑1 (MQ=255) tccatATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGttt > 2:5139081/1‑49 (MQ=255) ccatATAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 1:4136067/1‑49 (MQ=255) catATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 2:5198442/1‑48 (MQ=255) catATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 2:3096389/1‑48 (MQ=255) atatAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 2:4497427/1‑47 (MQ=255) atatAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 2:4586511/1‑47 (MQ=255) atatAAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAg < 1:251957/50‑1 (MQ=255) atatAAAGAGGGGATGGAGAAAAACGTTCGCCGACTGGTTCCGTTTGAAg < 2:3543301/50‑1 (MQ=255) ataAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTg > 1:3573645/1‑45 (MQ=255) taAAGTGTGGCTGGAGCAAAACGTCCGCCGACTAGTACCGTTTGAAGATc > 2:4020800/1‑50 (MQ=255) taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt < 2:1370190/49‑1 (MQ=255) taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt < 1:4985756/49‑1 (MQ=255) taAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATc < 2:1626354/50‑1 (MQ=255) taAAGAGTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGAt < 2:1857331/49‑1 (MQ=255) aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCt < 1:291480/49‑1 (MQ=255) aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg < 1:702219/50‑1 (MQ=255) aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg < 2:5344458/50‑1 (MQ=255) aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg > 2:4016607/1‑50 (MQ=255) aaGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTg < 1:1701149/50‑1 (MQ=255) agagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGc < 2:1189914/50‑1 (MQ=255) agagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGc > 1:521686/1‑50 (MQ=255) gagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGcc > 2:3190460/1‑50 (MQ=255) gagTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGcc > 2:2768335/1‑50 (MQ=255) agTGGATGGAGAAAAACGTCCGCCGTCTGGTACCGTTTGAAGATCTGccc > 2:2835974/1‑50 (MQ=255) agTGAATGGAGAAAAACGTCCGCCGACTGGCACCGTTTGAAGATCTGccc > 2:4721675/1‑50 (MQ=39) gTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGccc > 1:2212020/1‑49 (MQ=255) gTGAATGGAGAAAAACGTCCGTCGACTGGTACCGTTTGAAGATCTGCCCg < 2:3302219/50‑1 (MQ=39) gTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGccc > 1:4338106/1‑49 (MQ=255) gTGAATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCg < 1:1098603/50‑1 (MQ=255) ggatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGAt > 1:1207053/1‑50 (MQ=255) gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGAt > 2:4348835/1‑49 (MQ=255) gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATg < 2:1758408/50‑1 (MQ=255) | GCAGAAACCGATGACGATCTGAAAAGCCGCCATCCATATAAAGAGTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATG > NC_000913/3355952‑3356048 |
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: ‑ |