Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,161 | 0 | C | T | 58.2% | 0.9 / 75.1 | 67 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (7/21); new base T (39/0); total (46/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.13e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-08 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAG > NZ_CP009273/4594076‑4594247 | aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt > 2:15937/1‑90 (MQ=255) aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt > 2:1456022/1‑90 (MQ=255) cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCTCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 1:650006/90‑1 (MQ=255) cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg < 1:538107/90‑1 (MQ=255) ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttttc > 1:1708334/1‑87 (MQ=255) ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 2:1708334/87‑1 (MQ=255) gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag < 2:984853/90‑1 (MQ=255) gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag < 1:187932/90‑1 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTCTTTTTTTTCCCTgaga > 2:98723/1‑90 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGaa > 2:1199864/1‑90 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGaa > 1:469095/1‑90 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAg > 1:241832/1‑90 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTgaggag > 2:180875/1‑90 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCAt > 1:1512973/1‑84 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt < 2:1512973/84‑1 (MQ=255) cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAgg > 2:599065/1‑89 (MQ=255) aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt > 2:1688494/1‑75 (MQ=255) aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 1:1688494/76‑1 (MQ=255) ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc < 1:946676/90‑1 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGgcgc > 1:1613503/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCTTGAGGAGGGGgcga > 2:275944/1‑89 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTg > 1:967917/1‑72 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCa > 2:1580981/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc > 1:708793/1‑74 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCa > 2:719700/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg < 2:967917/72‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc < 2:708793/74‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATTGGAAGGGAGGGCa > 2:598882/1‑90 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCaga > 1:360955/1‑90 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGAGAGCGCaga < 1:1570980/90‑1 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga < 2:962996/90‑1 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga < 2:1344912/90‑1 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga < 1:98723/90‑1 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATTAGAAGGGGGGGCagaggg > 1:1096674/1‑88 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGGGAAGTGAGCGCagaggg > 2:1716399/1‑88 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAAGGAGCGGagagag > 2:253676/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt > 1:438587/1‑64 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 2:438587/65‑1 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGTATTTTTTTTTTTTTTGCCCTGAGGAGGGAGGggagggag > 2:1230383/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCCGAGAgg > 1:930813/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAgggg > 2:624214/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAgg > 1:513669/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGGGGGCAGAGgggg > 1:461942/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGGAGAGgggg > 1:973348/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAggg > 1:129396/1‑90 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGTGGTCATTTTTTCTTTTTCTGCCATg < 2:1388584/67‑1 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg > 1:1388584/1‑67 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTGTGCCATGAGGAGTGAGCGCAGAGAggggta > 1:644273/1‑87 (MQ=255) aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt < 1:719700/90‑1 (MQ=255) gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCGGAGAGGGGGTATc > 1:43126/1‑90 (MQ=255) gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATc > 1:1175561/1‑90 (MQ=255) gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGAGAAAGGGGGGCAGGGAGGGGGTATc > 1:687610/1‑90 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCtt > 2:1626683/1‑90 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGGGGGGTTATCtt > 2:290004/1‑90 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc > 1:324138/1‑59 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc < 2:324138/59‑1 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc < 2:901387/59‑1 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc > 1:901387/1‑59 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc < 1:290004/90‑1 (MQ=255) aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGGGGGGGTATCTTGCctga > 2:1665240/1‑90 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 2:1680726/45‑1 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc > 1:1680726/1‑45 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGAGGGCAGGGGGGGGTTATCTTTCcttactga > 2:1440592/1‑90 (MQ=255) gCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCATGAGGAGAGAGCGCAGAGAGGGGTTATCTTTCCTGACTGaa > 1:1720169/1‑90 (MQ=255) gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGGGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCggt > 1:97118/1‑88 (MQ=255) ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGGGGGGGGGGGGGGGTTTTTTTGCCTTACTTGACTTGCGAttt > 2:598268/1‑89 (MQ=255) ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTa > 1:1416466/1‑90 (MQ=255) tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGGAATGAGCTCTGAGGGGGGTTATTTTGCCTGACTTGACTGGCGTTTaaat > 1:531558/1‑88 (MQ=255) tcGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGGGGGGTTATCTTGGCTGACTTAACTGGCGGTTTAGttt > 1:882298/1‑90 (MQ=255) tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt < 2:469095/90‑1 (MQ=255) ttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc < 1:459370/90‑1 (MQ=255) gcggcgGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTcgcgg > 1:1595185/1‑89 (MQ=255) gcggcgGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTTTTGCCTGACTGAACTGGGGATTAAGTTTTcgcgg > 1:1645105/1‑89 (MQ=255) gcggcgGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGt > 1:32488/1‑90 (MQ=255) gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc < 1:1559829/90‑1 (MQ=255) gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc < 1:1230383/90‑1 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAAAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 2:937318/90‑1 (MQ=255) tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc > 1:978515/1‑90 (MQ=255) ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa < 1:447847/90‑1 (MQ=255) tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg > 1:655414/1‑90 (MQ=255) | AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAG > NZ_CP009273/4594076‑4594247 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |