Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,754,231 | 0 | G | T | 28.6% | 78.5 / 26.5 | 35 | R127S (CGC→AGC) | selA | L‑seryl‑tRNA(Sec) selenium transferase |
Reads supporting (aligned to +/- strand): ref base G (12/13); new base T (5/5); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCT > NZ_CP009273/3754154‑3754319 | gTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCgtcgtc < 2:200700/90‑1 (MQ=255) tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg > 2:355991/1‑88 (MQ=255) tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg < 1:355991/88‑1 (MQ=255) aTACAGGCATCTTCCGCCCCCGTAATACGACACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc < 2:232645/81‑1 (MQ=255) aTACAGGCATCTTCCGCCCCCGTAATACGACACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc > 1:232645/1‑81 (MQ=255) gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcg < 2:283533/60‑1 (MQ=255) gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcg > 1:283533/1‑60 (MQ=255) ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 1:38707/90‑1 (MQ=255) ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 1:155904/90‑1 (MQ=255) ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 2:563134/90‑1 (MQ=255) ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa < 2:546226/68‑1 (MQ=255) ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa > 1:546226/1‑68 (MQ=255) ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 2:490250/88‑1 (MQ=255) ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa > 1:490250/1‑88 (MQ=255) cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc < 1:504768/61‑1 (MQ=255) cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc > 2:504768/1‑61 (MQ=255) ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg < 1:590530/55‑1 (MQ=255) ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg > 2:590530/1‑55 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgt > 1:573853/1‑47 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgt < 2:573853/47‑1 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgt < 2:125465/47‑1 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgt > 1:125465/1‑47 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg > 1:377055/1‑90 (MQ=255) gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg > 2:140794/1‑90 (MQ=255) tGCACCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt > 2:326866/1‑90 (MQ=255) gATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt < 1:221429/90‑1 (MQ=255) aTCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTc > 2:310977/1‑90 (MQ=255) cgcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCg < 1:310977/90‑1 (MQ=255) gcgATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa > 2:118718/1‑51 (MQ=255) gcgATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 1:118718/51‑1 (MQ=255) ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa < 1:77830/44‑1 (MQ=255) ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa > 2:77830/1‑44 (MQ=255) ccGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc > 2:46388/1‑81 (MQ=255) ccGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc < 1:46388/81‑1 (MQ=255) ccGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa > 1:226571/1‑90 (MQ=255) cgTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCt < 2:226571/90‑1 (MQ=255) | GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCT > NZ_CP009273/3754154‑3754319 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |