Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 1,132,851 | 0 | A | G | 20.0% | 359.3 / 16.9 | 151 | N65N (AAT→AAC) | yceI | YceI family protein |
Reads supporting (aligned to +/- strand): ref base A (27/93); new base G (30/0); total (58/93) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.36e-16 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 190 total aligned reads displayed.
TTATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTG > NC_012967/1132780‑1132922 | ttATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGa < 1:4535456/72‑1 (MQ=255) tATCGNGTTCGGCGTGATTAGTATCGNCGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAt < 2:842815/72‑1 (MQ=255) tATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAt > 2:5237041/1‑72 (MQ=255) aTCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGagt > 1:5618960/1‑70 (MQ=255) aTCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGtt > 2:5737966/1‑72 (MQ=255) aTCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAtt < 2:2114886/72‑1 (MQ=255) tCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGttt > 2:7280951/1‑72 (MQ=255) cgcgTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGtttt > 1:6902179/1‑72 (MQ=255) cgcgTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAtttt < 1:5923394/72‑1 (MQ=255) cgcgTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAtttt < 2:7326021/72‑1 (MQ=255) gcgTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAttttt < 2:6811203/72‑1 (MQ=255) gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGgtttttt > 2:6642658/1‑72 (MQ=255) gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAttttttt < 2:5021796/72‑1 (MQ=255) gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAttttttt < 2:4550940/72‑1 (MQ=255) gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAttttttt < 2:1606017/72‑1 (MQ=255) gTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGAttttttt < 2:305035/72‑1 (MQ=255) ttCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGgttttttt > 1:1543388/1‑71 (MQ=255) ttCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTATTTTTc > 1:3178694/1‑72 (MQ=255) ttCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTc < 2:7033761/72‑1 (MQ=255) tCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTGACTTTATCGGCC‑GCCGGGgtttttttg > 2:4471540/1‑70 (MQ=255) cGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCAGGCCGGTattttttgt > 2:6309252/1‑70 (MQ=255) cGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGt < 1:4722367/72‑1 (MQ=255) ggCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTc < 2:5958668/72‑1 (MQ=255) ggCGTGATTAGTATCGACGCTGGGGGTGTTAATTGTCACATTCACTTTATCGGCA‑GGCGGGTTTTTTTTGtg > 1:2559381/1‑71 (MQ=255) agCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTc > 1:127257/2‑72 (MQ=255) tGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTT‑aaaa > 2:3015850/1‑71 (MQ=255) tGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAg < 2:4571473/72‑1 (MQ=255) atgaGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGt < 1:4447929/69‑1 (MQ=255) aTTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGGTTTTTTTGTC‑AAAGGt > 1:4228282/1‑72 (MQ=255) aTTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGCTTTTTTTCGTC‑TAAggg > 1:3920224/1‑71 (MQ=255) gTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCNGATTTTTTTCGTC‑AAAGGTaaaa < 1:4832490/72‑1 (MQ=255) tATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTCGTCAAAAGGTaaaa > 2:5559162/1‑72 (MQ=255) aTCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTCTTC‑AAAGGGAAAAga > 2:1306519/1‑71 (MQ=255) aTCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGt < 2:5786775/72‑1 (MQ=255) aTCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCCGAATTTTTTTTTC‑CAAGGTnaaata > 2:2883330/1‑70 (MQ=255) cGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGATTTTTTTTTC‑CAAGGTAAAAGTaa > 2:3927661/1‑71 (MQ=255) cGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTAc < 2:2600507/72‑1 (MQ=255) cGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCAGGGTTTTTTTCGGC‑AAAGTTGAAAGTAc > 2:32286/1‑72 (MQ=255) gACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTAcc < 1:1028533/72‑1 (MQ=255) gACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTAcc < 1:3682171/72‑1 (MQ=255) gACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTAcc < 1:354393/72‑1 (MQ=255) cggtggtggGGTTAATTGGCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGt < 1:5967886/69‑1 (MQ=255) gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGTA‑NCCGGATTTTTTTCGTC‑ANNNGTAAAAGTACCGTc > 1:5328904/1‑72 (MQ=255) gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTT‑CAAGGTAAAAGTAACGTc > 2:4122551/1‑72 (MQ=255) gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTCGTC‑AAAGGTAAAAGTACCGTc > 1:7130940/1‑72 (MQ=255) gCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTc < 2:6786990/72‑1 (MQ=255) cTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTTTC‑ACAAGAAAAAATACCGTCg > 2:4952450/1‑72 (MQ=255) cTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCg < 1:5858195/72‑1 (MQ=255) ggtggtGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTTTC‑AAAGGTAAAAATACCCTTGaa > 1:4253059/1‑72 (MQ=255) ggtggtGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGGTTTTTTCGTC‑AAAGGGAAAAGTAACGTCgga > 2:6849828/1‑70 (MQ=255) ggtggtGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGaa < 2:7230021/72‑1 (MQ=255) ggtggtGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGaa < 1:5847337/72‑1 (MQ=255) gtggtgTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAGGTAAAAGTACCGTCGaaa > 2:6251988/1‑72 (MQ=255) tggtgTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTT‑TAAGGTAAAAATAACGTCGAAAt > 2:2396885/1‑72 (MQ=255) tggtgTTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAGGGAAAAGTACCGTCGaaaa > 2:1855051/1‑71 (MQ=255) tggtgTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAAt < 1:3242679/72‑1 (MQ=255) tggtgTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAAt < 1:6357587/72‑1 (MQ=255) gtgtTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑t < 1:5414981/72‑1 (MQ=255) tgtTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑tt < 1:1747978/72‑1 (MQ=255) tgtTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑tt > 1:1246833/1‑72 (MQ=255) gtTAATTGTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAGGTAAAAGTACCGTCGAAATC‑ttt > 1:41728/1‑72 (MQ=255) ttAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTa < 2:3647344/72‑1 (MQ=255) aaTTGTCACATTCACTTTNTCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTaaa < 1:2892055/72‑1 (MQ=255) aaTTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTaaa > 1:492069/1‑72 (MQ=255) aTTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAg < 1:1857281/72‑1 (MQ=255) aTTCTCACATCCACTTTATCGGCA‑GCCGAATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAg < 1:4674793/72‑1 (MQ=255) tGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGt < 1:7355169/72‑1 (MQ=255) tGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGt < 1:6207161/72‑1 (MQ=255) gTCACATTCACTTTATCGGCA‑GCCGGGTTTTTTTCGTC‑AAAAGTAAAAGTACCGTCGAAATC‑TTTAAAGGTg > 1:5994346/1‑72 (MQ=255) gTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTg < 1:2993368/72‑1 (MQ=255) cacaTTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAAGGAAAAATACCGTCGAAATC‑TTTAAAAGGGcc > 2:7127195/1‑72 (MQ=255) cacaTTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGcc < 2:4911631/72‑1 (MQ=255) caTTCACTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAGGGAAAAGTACCGTCGAAATCTTTTAAAGGGGcgt > 2:5412222/1‑70 (MQ=255) caTTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGt < 1:2201209/72‑1 (MQ=255) caTTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGt < 2:7378612/72‑1 (MQ=255) tCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGtata < 1:5692894/72‑1 (MQ=255) cACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATaa > 1:7523608/1‑72 (MQ=255) aCTTTATCGGCA‑GCCGGGTTTTTTTTGTC‑AAAGGGAAAAGTACCGTCGAAATC‑TTTAAAAGTGCCGTATAAc > 2:1523630/1‑72 (MQ=255) aCTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAAc < 1:5840307/72‑1 (MQ=255) ttATCGGCA‑GCCGGGTTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAg > 1:3609514/1‑72 (MQ=255) tATCGGCA‑GCCGGGTTTTTTTTGTC‑CAACGGAAAAATACCCTCGAAATC‑TTTTAAAGTTCCGTAATACCAGc > 2:993772/1‑72 (MQ=255) tATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGc < 2:3880056/72‑1 (MQ=255) tATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGc < 2:1272082/72‑1 (MQ=255) aTCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCt < 2:6527590/72‑1 (MQ=255) ggCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCtata < 1:1346924/72‑1 (MQ=255) cA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGc < 2:4610568/72‑1 (MQ=255) cA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGc < 2:3780602/72‑1 (MQ=255) cA‑GCCGGA‑TTTTTTCGTC‑AAAGGTAAAAGTNCCGTCGAAATC‑TTNAAAGGTGCCGTATAACCAGCTATAGcc < 1:373720/72‑1 (MQ=255) a‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGcc > 2:5597623/1‑72 (MQ=255) a‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGcc > 1:2779081/1‑72 (MQ=255) gCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCa < 2:3211128/72‑1 (MQ=255) cGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAg > 2:1254952/1‑72 (MQ=255) cGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAg < 1:1504463/72‑1 (MQ=255) ggATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAgg < 1:101808/72‑1 (MQ=255) ggATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAgg > 2:5293228/1‑72 (MQ=255) gATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGt > 1:279649/1‑72 (MQ=255) aTTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTg < 1:6271044/72‑1 (MQ=255) aTTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTg < 1:238743/72‑1 (MQ=255) aTTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTg < 1:7191379/72‑1 (MQ=255) aTTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTg < 2:1579507/72‑1 (MQ=255) | TTATCGCGTTCGGCGTGATTAGTATCGACGCTGGTGGTGTTAATTGTCACATTCACTTTATCGGCA‑GCCGGATTTTTTTCGTC‑AAAGGTAAAAGTACCGTCGAAATC‑TTTAAAGGTGCCGTATAACCAGCTATAGCCAAGGTG > NC_012967/1132780‑1132922 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |