Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,028,269 | G→C | 100% | P126A (CCA→GCA) | USA300HOU_RS10230 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,028,269 | 0 | G | C | 96.7% | 95.2 / ‑2.8 | 30 | P126A (CCA→GCA) | USA300HOU_RS10230 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base C (11/18); total (11/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAC > USA300TCH1516_ALE/2028134‑2028402 | ccTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATATACTGCTGTAc < 1:554687/141‑1 (MQ=255) cACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTg < 1:389576/141‑1 (MQ=255) aaaTGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAAcc < 1:655435/141‑1 (MQ=255) aaTGTGTCTAACAACAGTGGTATTATTTTCCTTAATATCTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTAATATACATAAACTGCTGTACGTTGCTCACCAAACCa > 1:424274/1‑141 (MQ=255) caacaaTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAACTAACTATCTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTtaa < 1:750029/141‑1 (MQ=255) ggTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGACAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGa < 1:64715/141‑1 (MQ=255) gTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGAt < 1:490574/141‑1 (MQ=255) aGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGt > 1:297829/1‑141 (MQ=255) gCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTc > 1:57908/1‑141 (MQ=255) ttttGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTGTATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAGCCTCAGCcgccgc > 1:254560/1‑141 (MQ=255) aaGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactact < 1:102675/141‑1 (MQ=255) gTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACttt > 1:667623/1‑141 (MQ=255) tCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtttt < 1:764992/141‑1 (MQ=255) tatGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTg > 1:308281/1‑141 (MQ=255) aGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGCCGCAACTACTTTTTTTGGTTTa > 1:595595/1‑141 (MQ=255) cAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCACCTACTTTTTTTGGTTTAAACTtgt > 1:44771/1‑141 (MQ=255) gTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAAt < 1:687711/141‑1 (MQ=255) aTTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATTGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGc < 1:479999/141‑1 (MQ=255) ttGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGcac < 1:767290/141‑1 (MQ=255) tACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAAt > 1:61038/1‑141 (MQ=255) tACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAAt < 1:50184/141‑1 (MQ=255) tCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGt < 1:661357/141‑1 (MQ=255) cctCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGta < 1:254838/139‑1 (MQ=255) ttCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaa < 1:360696/141‑1 (MQ=255) ttCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaa < 1:751062/141‑1 (MQ=255) cACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTTTTTTTTTTTTGTGTATTACATAACGCACAATACGtaatact > 1:27583/1‑139 (MQ=255) ataCAGAAACTGCTGGACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACTTCATCCGTCGCAACTTCTTTTTTTGGTTTAAACTTGTGTCTTCGGTCTCTCACAATAAGTAATAATAAAAAGa < 1:182407/141‑1 (MQ=255) taCATAAACTGCTGTACGTTGCTCACCAAACCACTTATGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATACGAAATAATAAAAAGAc < 1:423664/141‑1 (MQ=255) cATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGa > 1:746786/1‑141 (MQ=255) tacaCTGCTGTACGATGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAc < 1:72338/138‑1 (MQ=255) | CCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAAC > USA300TCH1516_ALE/2028134‑2028402 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |