Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 276,387 | 0 | T | C | 25.3% | 174.1 / 62.5 | 91 | G20G (GGA→GGG) | yagB | type IV toxin‑antitoxin system antitoxin YagB |
Reads supporting (aligned to +/- strand): ref base T (29/39); new base C (12/11); total (41/50) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTC > NZ_CP009273/276300‑276477 | tAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCg < 2:43362/91‑1 (MQ=255) gTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAac > 2:298425/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:723065/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 2:32252/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:22417/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 2:335767/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 1:735374/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:17850/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 2:491352/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 2:179818/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 2:193470/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:377604/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:452552/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac < 1:605006/91‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 2:434083/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAacac > 1:352742/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAac < 2:219256/89‑1 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAac > 1:219256/1‑89 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAacac > 1:183040/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAacac > 1:175692/1‑91 (MQ=255) cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAacac > 1:91083/1‑91 (MQ=255) tGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGATCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACg < 1:785643/91‑1 (MQ=255) gCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACgg > 2:508262/1‑91 (MQ=255) aCCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGgtg < 1:352301/91‑1 (MQ=255) cACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACg > 1:463306/1‑84 (MQ=255) cACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACg < 2:463306/84‑1 (MQ=255) gAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATAtct > 2:245102/1‑91 (MQ=255) gAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATAtct > 2:455601/1‑91 (MQ=255) gAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATAtct < 2:223736/91‑1 (MQ=255) gAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATAtct > 2:4895/1‑91 (MQ=255) gAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATAtct > 1:276093/1‑91 (MQ=255) aGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGttt > 2:842727/1‑91 (MQ=255) aGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGttt > 1:149851/1‑91 (MQ=255) aGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGttt > 1:78681/1‑91 (MQ=255) aGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGttt > 1:802639/1‑91 (MQ=255) gCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTa > 2:835083/1‑91 (MQ=255) gCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTa < 2:206718/91‑1 (MQ=255) gCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACgg > 2:88396/1‑75 (MQ=255) gCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACgg < 1:88396/75‑1 (MQ=255) ccGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTaa < 1:508262/91‑1 (MQ=255) cGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAg < 1:347582/91‑1 (MQ=255) gTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGatat < 2:638414/79‑1 (MQ=255) gTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGatat > 1:638414/1‑79 (MQ=255) aaTCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGccccc < 1:173318/91‑1 (MQ=255) aaTCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGccccc < 1:732247/91‑1 (MQ=255) aaTCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGccccc < 2:89237/91‑1 (MQ=255) tCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCaa < 1:853362/91‑1 (MQ=255) tCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCaa < 2:83256/91‑1 (MQ=255) aGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCGCCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGtt < 2:442776/91‑1 (MQ=255) cccGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTg < 1:434083/91‑1 (MQ=255) gTCAGCCAGAAAATGCAGACGATGCCCTTCCTGCACCAGTCGTGCTCCGAACCACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTa > 1:102830/1‑91 (MQ=255) gTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTa > 2:132354/1‑91 (MQ=255) gTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTa > 2:528438/1‑91 (MQ=255) cAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTaga < 2:232193/91‑1 (MQ=255) cAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTaga < 2:449590/91‑1 (MQ=255) cAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTaga < 1:564395/91‑1 (MQ=255) cAGCCAGAAAATGCAGACGATGACCTTCCTGAACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTaga < 1:706614/91‑1 (MQ=255) aaTGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTAtt < 2:447551/91‑1 (MQ=255) gCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCAt < 2:379822/91‑1 (MQ=255) gCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCAt < 2:6205/91‑1 (MQ=255) cAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATa < 2:245933/91‑1 (MQ=255) cAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATa < 2:78681/91‑1 (MQ=255) cAGACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATa < 1:324739/91‑1 (MQ=255) gACGATGACCTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAAt < 2:599562/91‑1 (MQ=255) aCGATGACCTTCCTGCACCAGTCGTGCTCCGAAac < 1:671451/35‑1 (MQ=255) aCGATGACCTTCCTGCACCAGTCGTGCTCCGAAac > 2:671451/1‑35 (MQ=255) tGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTAt > 1:565013/1‑91 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 1:736732/91‑1 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 1:765379/91‑1 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 1:290243/91‑1 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 2:286823/91‑1 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 1:402505/91‑1 (MQ=255) aCCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 2:413676/91‑1 (MQ=255) ccTTCCTGCACCAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTcc < 2:831783/91‑1 (MQ=255) ttCCTGCACCAGTCGTGCTCCGAAACACGGTGGGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTg > 2:425236/1‑91 (MQ=255) tCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGt < 1:239303/91‑1 (MQ=255) ccAGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAaagaa > 1:406933/1‑91 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTag < 2:576789/55‑1 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTag > 1:576789/1‑55 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc < 2:352471/78‑1 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTc > 1:352471/1‑78 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGt > 2:334793/1‑82 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGt < 1:334793/82‑1 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAaag < 2:788856/87‑1 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAaag < 2:79404/87‑1 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAaag > 1:79404/1‑87 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAaag > 1:788856/1‑87 (MQ=255) aGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAagaaga > 1:81252/1‑91 (MQ=255) aGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAagaaga > 1:81745/1‑91 (MQ=255) aGTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAagaaga > 1:290555/1‑91 (MQ=255) gTCGTGCCCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGAc > 1:51167/1‑91 (MQ=255) tCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTc > 2:862571/1‑91 (MQ=255) | TAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTTC > NZ_CP009273/276300‑276477 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |