Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 228,181 | 0 | C | T | 60.6% | 7.6 / 36.4 | 33 | F191F (TTC→TTT) | yafD | endonuclease/exonuclease/phosphatase family protein |
Reads supporting (aligned to +/- strand): ref base C (1/12); new base T (20/0); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.66e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-06 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA > NZ_CP009273/228094‑228261 | gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa < 2:401758/90‑1 (MQ=255) ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAAt < 2:204184/90‑1 (MQ=255) ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTTATTCCTg > 1:386489/1‑90 (MQ=255) aTGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGa < 2:964544/90‑1 (MQ=255) aTGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGa < 1:609761/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc < 2:1049163/90‑1 (MQ=255) tAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAATCCTGGAGGCGGAGAAg > 2:562949/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATTCCCGGGGCCGGAGAAgg > 1:1096117/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATTCCTGGAGGCGGAGAAgg > 1:131177/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTTAATGCCCTGGGGCGGAGAAgg > 2:882921/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTaggagg > 2:1080174/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCCTGGGCCGCGGAAGgaagaa > 1:900080/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAAATCCCGGGGCCGGAGAAGGATGaa > 2:1222388/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGaa < 2:821083/90‑1 (MQ=255) cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTTAATGCCTGGGGCCGGAGAAGGGTGAAAGCGGtat > 2:106869/1‑90 (MQ=255) tGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGGAGAAGGAAGAAAGCGGTATATCg > 2:6071/1‑90 (MQ=255) tGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGCCGTAGAAGGATGAAAGCGGTATATCg > 2:461810/1‑90 (MQ=255) gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCttt < 1:489850/90‑1 (MQ=255) aTCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTTATGCCCGGGGCCCCAGAAGGGTGAAAGCGGTATATCGCcttt > 2:1082633/1‑89 (MQ=255) tCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCTGGGGGCGGGGAAGGGGGGAAGCGGTATATTGCttttt > 1:659580/1‑88 (MQ=255) tAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGGAGAAGGGTGAAAGCGGTATATCGCTTTTCGCggg > 2:28667/1‑90 (MQ=255) gCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGATTTAATTCCTGGGGGCGGAGAAGGATGAAAGCGGTATATTGCTTTTCGCgggga > 2:613054/1‑88 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCCGGGGGCGGGGGAGGGGGGAAGCGGTTTATTGCTTTTCGCGGGgaaag > 1:1114583/1‑88 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTTGGGCCGGAGAAGGATGAAAGCGTTATATTGCTTTTGGCGGGgaaag > 2:302187/1‑88 (MQ=255) ccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGTAGAAGGAAGAAAGCGGTATATTGCTTTTCGCGGGGAAAGtt > 2:915456/1‑89 (MQ=255) caGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCt < 2:744629/90‑1 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGGTTTTAATGCCTGGGGCCGGAGAAGGGGGAAAGCGGTTTATTGCTTTTTGCGGGGAAAGTCGCTgc > 1:481942/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGGATTTAATGCCTGGGGCCGGAGAAGGGGGAAAGCGGTATATTGCTTTTTGCGGGGAAAGTCGCTgc > 1:279974/1‑90 (MQ=255) ccGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGGAGAAGGGGGAAAGCGTTATATTGCTTTTCGCGGGGAAAGTCGCTGCGcccg > 1:64515/1‑88 (MQ=255) aTTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGt < 2:475226/90‑1 (MQ=255) aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt < 1:461810/90‑1 (MQ=255) gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 2:1073190/90‑1 (MQ=255) gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 2:917051/90‑1 (MQ=255) | GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA > NZ_CP009273/228094‑228261 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |