Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,012,951 | 0 | A | T | 39.3% | 36.1 / 18.6 | 28 | L144F (TTA→TTT) | ssnA | putative aminohydrolase SsnA |
Reads supporting (aligned to +/- strand): ref base A (7/10); new base T (11/0); total (18/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGAC > NZ_CP009273/3012819‑3013079 | agagTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGtt < 2:523410/139‑1 (MQ=255) gTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGTCCTCTCCCCATTTCGC‑GACC‑CATCTTTATAAGTTggg > 2:85231/1‑138 (MQ=255) cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGAcc < 1:404989/139‑1 (MQ=255) gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGTCCCCCCCATGCGC‑GGCG‑CCTTTTTTTTAAATGGCGCTCGCGCGGTGACCTg > 2:404989/1‑139 (MQ=255) gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACACTTCGCGGCCC‑CTTTTTTTAAAGTT‑GCCTTCGCGCGCTGGCCTg > 2:659485/1‑139 (MQ=255) gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCACCTTTCGC‑GCCG‑CATTTTTTTAAATTGGGCGGCGCGCGCGGACCtc > 1:431542/1‑138 (MQ=255) tttGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGa < 1:659485/139‑1 (MQ=255) ttCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCCCCCTGCGC‑GGCG‑CCTTTTTTTAAATTTGCCTGCGCGCGCTGCCCTGtttttttac > 2:557743/1‑133 (MQ=255) gAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc < 1:274687/139‑1 (MQ=255) agCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCCATTCGC‑GGCGCCCTTTTTTAAAGTGGGCC‑GCGCGCGGGGGCCCGTTTTTGAACTACTTGCCCTaacaaaagc > 2:10511/1‑134 (MQ=255) gATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGTCCTCTCCCTTCGC‑GCCG‑CCTTTTTTTAAATTTGCCTGCCCCCGCTGTACTCTTTTGTAACGACTTAACCTACCAACGCCagcc > 2:176594/1‑135 (MQ=255) aTGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCaaa < 2:613380/139‑1 (MQ=255) gTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAgg < 1:557743/139‑1 (MQ=255) ttATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCCCCTTTCGC‑GACG‑CCTTTTTTTAAATTTGCCCTCGCGCGCTGGCCCCCTTTTTAACAAATGTCCGCAACAACCGCAACAAAAAGTTGGAgg > 2:364682/1‑138 (MQ=255) tCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCGCTTCACATTGACC‑GCCG‑CCCTTTTTTAAGATGGCCTGCGCGCGGTGCCCTGTTTTTTAAATACTATTCGTAAAAAACGCATCACAGAGGTGTCGGAAgg > 2:258878/1‑139 (MQ=255) ccATGCCTCTCCGGCGTATATCGGCGGGGCGGTCTCCCCATTCCGC‑GCCG‑CCTTTTTTAAAATTGGCCTGCGCGCGGTGAACTGGTTTTAAACTACTTTCCGGAACAACAGCATCAAAGAGATGCAGGAAgggggag > 2:189695/1‑132 (MQ=255) tGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa < 1:210450/139‑1 (MQ=255) tGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa < 1:75553/139‑1 (MQ=255) tGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACAATGCGC‑GAAG‑CACTTTTTTAAGTTGGGCTGCGCGCGGGGAGCTGTTTTTAAAATAATGCACGGAAAAACCGCAACCAAAAGGGGGAGGAAGGGGGggaaga > 1:517354/1‑138 (MQ=255) ctcCGGCGTATATCGGCGGGGAGCTCTCCACACTGCGC‑GGCG‑CAATTTTTAAAAATGGGCTGGGCGCGGGGAACGGGTTTGAAACTACTGACCGTGACAACGGCATCAAAAAGTTGCAGGAAGGAGTGGAAGAAAACaa > 1:573418/1‑138 (MQ=255) gCTCACCACATTGCGA‑TACG‑CATTTTTAAAAACGGACATGCGCGCGATGACCTGTTATGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAATACATCCGTTTCGCCCGTTtgattgat < 2:89121/139‑1 (MQ=255) tGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCgaaga > 1:632322/1‑138 (MQ=255) tGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCgaaga > 1:43152/1‑138 (MQ=255) tGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCgaaga > 1:225593/1‑138 (MQ=255) gc‑gACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGc > 1:413661/1‑139 (MQ=255) gACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGa < 2:431542/139‑1 (MQ=255) gACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGAAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGa > 1:64562/1‑139 (MQ=255) cG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGAc > 1:408797/1‑138 (MQ=255) | AGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGAC > NZ_CP009273/3012819‑3013079 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |