| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | Exported | 75,335 | 2 | . | A | 33.3% | 19.8 / 15.6 | 12 | intergenic (‑76/‑1018) | BCHEFKBB_00065/BCHEFKBB_00066 | Opacity protein opA54 [NGO_11100]/peptidase A8 [NGO_0071] |
| Reads supporting (aligned to +/- strand): ref base . (5/3); new base A (2/2); total (7/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTCCGGATAATCGTGGGTAATGCGTTCGGCGGCGTAGGCTAAATCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACAAAAACC > Exported/75205‑75484 | gttCCGGATAATCGTGGGTAATGCGTTCGGCGGCGTAGGCTAAATCCGCCTGCACATACGGGCCGCGGCCATTTCCTTCAATTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTg > 2:1908818/1‑150 (MQ=255) aTCGTGGGTAATGCGTTCGGCGGCGTAGGCTAAATCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGAtt > 1:976942/1‑151 (MQ=255) aTGCGTTCGGCGGCGTAGGCTAAATCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGACGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGc < 2:976627/151‑1 (MQ=255) ggCTAAATCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGACGAGAAGAGAAGGTTTGTTGCGGGCTGGATTCCTTTTCGGCTCCTTATTCGGt > 2:1132037/1‑151 (MQ=14) ctatctCCGCCTGCACCTCCGGGCCGCGGCCAGTGACTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCg < 2:1106659/146‑1 (MQ=18) tCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGACGAGAAGAGAAGGTTTTTTGCGGGCGGGATTCATTTTCGGCTCCTTATTCGGTTTAACCg < 2:1762300/151‑1 (MQ=17) cGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTaaaa > 2:1667230/1‑151 (MQ=255) ggCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGACGAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTCTTCGGTTTATCCGGTTAAAAAAAGAtttt > 2:1121634/1‑151 (MQ=14) ggCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCa > 1:1266614/1‑148 (MQ=255) ggCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCa > 1:1950521/1‑148 (MQ=255) cTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGACGAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAgg < 2:661835/151‑1 (MQ=11) aagcGAAGAGAAGAGAAGAGAAGAGCAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCggg < 2:162946/146‑1 (MQ=11) aagagaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAac < 2:1004922/151‑1 (MQ=11) aagagaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAac < 1:861495/151‑1 (MQ=11) aagagaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAac < 1:861478/151‑1 (MQ=11) aagagaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTGTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAac < 1:560987/151‑1 (MQ=11) agagaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAaca > 2:453760/1‑151 (MQ=31) agaagagaagagaagagaagagaa‑‑‑‑‑gagacgagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAacac < 1:371272/150‑1 (MQ=17) agaagagaagagaagagaagagaa‑‑‑‑‑gagaagagatgGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAacac < 1:1827557/150‑1 (MQ=17) agaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAacac < 1:371128/150‑1 (MQ=32) agaacagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAacac < 1:1859940/150‑1 (MQ=17) gaagagaagagaagagaatagaa‑‑‑‑‑tagaagagaagGTTTTTTGCGGGCTGGATTCTTTATCGGCTCCTTATTCGGTGTAACCGGTTAAAAAATGATTTTCAATGATGTTGAAGGGCGGAGTAGATCGGGTTCCGCGCGGTGTGTCAACACaa > 2:872494/1‑151 (MQ=11) gaagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACaa > 2:836844/1‑151 (MQ=32) aagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAat < 1:1522764/151‑1 (MQ=32) aagagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTAGTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAat < 1:1522759/151‑1 (MQ=17) gagaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat > 2:597491/1‑151 (MQ=32) agcagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat < 2:1778887/147‑1 (MQ=17) agaagagacgagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat < 1:474152/150‑1 (MQ=17) agaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat < 1:948994/150‑1 (MQ=32) agaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat < 1:658864/150‑1 (MQ=32) agaagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAatat < 1:678457/150‑1 (MQ=32) aagggaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 1:1099001/147‑1 (MQ=17) aagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 2:1703478/151‑1 (MQ=32) aagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 2:1698862/151‑1 (MQ=32) aagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 2:1235447/151‑1 (MQ=32) aagagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 1:74574/151‑1 (MQ=32) gagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAGTATGGCgg > 2:1957600/1‑151 (MQ=17) gagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg > 1:510940/1‑151 (MQ=32) gagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg > 1:602391/1‑151 (MQ=32) gagaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg > 1:1043196/1‑151 (MQ=32) agaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGc < 2:1235455/148‑1 (MQ=32) gaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg > 2:1543418/1‑149 (MQ=32) gaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg > 2:1543547/1‑149 (MQ=32) gaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg < 1:1544081/149‑1 (MQ=32) gaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCgg < 1:1543952/149‑1 (MQ=32) gaagagaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGAt > 1:628084/1‑151 (MQ=32) gaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACa > 2:1020041/1‑151 (MQ=33) gaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACa > 2:1020039/1‑151 (MQ=33) gaagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAAAATGGCGGATGAACa > 2:1582147/1‑151 (MQ=18) aagagaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACaa < 1:120416/151‑1 (MQ=33) gaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCTGATGAACaaaaa > 2:1539139/1‑150 (MQ=18) gaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACaaaaa > 1:910567/1‑150 (MQ=33) gaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACaaaaa > 1:1130272/1‑150 (MQ=33) gaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACaaaaa > 1:1130263/1‑150 (MQ=33) gaa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACaaaaa > 1:928907/1‑150 (MQ=33) aa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGTATTCATTTTCGGATCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACAAAAAcc < 2:760981/151‑1 (MQ=16) aa‑‑‑‑‑gagaagagaagGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGATAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACAAAAAcc < 1:445456/151‑1 (MQ=18) | GTTCCGGATAATCGTGGGTAATGCGTTCGGCGGCGTAGGCTAAATCCGCCTGCACATACGGGCCGCGGCCATTGCCTTCACTTGCCGCCTGCGCTGCCGAAGAGAAGAGAAGAGAAGAGAAGAGAAGAGAA‑‑‑‑‑GAGAAGAGAAGGTTTTTTGCGGGCTGGATTCATTTTCGGCTCCTTATTCGGTTTAACCGGTTAAAAAAAGATTTTCACTGATGTTGAAGGGCGGATTATATCGGGTTCCGGGCGGTGTTTCAACACAATATGGCGGATGAACAAAAACC > Exported/75205‑75484 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |