Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 60.0% | 15.4 / 27.5 | 25 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (9/1); new base C (0/15); total (9/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336711‑2336874 | ttATTCCAAATAACCAGCTTTTTCCCGTCGGAAATAAAAAAGGAAACTCCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGaggaa < 1:684830/89‑1 (MQ=255) gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa > 1:244524/1‑90 (MQ=255) gCCAATACCCAGCTTTTTTCCGTCGGAAATACAAAAGGTAACTTCCCTAAACCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGAcc < 2:878378/90‑1 (MQ=255) gCCAATACCCAGCTTTTTCCCGCCGGAAATATAAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGAcc < 2:808122/90‑1 (MQ=255) cagaTACCCAGCTTTTTTCCGTCGCCAATACCAAAGGTACCTTCACTTAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCa < 1:71112/87‑1 (MQ=255) aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc > 1:194106/1‑90 (MQ=255) aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc < 2:322589/90‑1 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 2:526468/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:839304/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:545532/1‑90 (MQ=255) tttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctg > 2:71112/1‑90 (MQ=255) ttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:347750/90‑1 (MQ=255) tCCGTCGGAAATATCAAAGGTACCTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc < 2:320125/90‑1 (MQ=255) cGGAAATAAAAAAGGTAACTTCCCTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt < 2:393589/90‑1 (MQ=255) ggTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGAc > 1:225097/1‑90 (MQ=255) aCTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGcc > 1:878516/1‑90 (MQ=255) ttaaaGCCGAAAAAAATAACCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCt < 1:481988/87‑1 (MQ=255) taaaGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTc < 2:545532/88‑1 (MQ=255) ccGAGATACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 2:10738/90‑1 (MQ=255) ccGAAATAAATAAGGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 2:294595/90‑1 (MQ=255) ccGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 1:70894/90‑1 (MQ=255) gagaTACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg < 2:225097/90‑1 (MQ=255) gagaTACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg < 2:495035/90‑1 (MQ=255) cATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTa > 1:64104/1‑90 (MQ=255) tAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGc < 2:653041/90‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt > 1:137621/1‑44 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt < 2:137621/44‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 2:879586/1‑90 (MQ=255) | GTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336711‑2336874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |