Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,097,080:1 | +A | 100% | intergenic (‑104/+46) | azoJ ← / ← yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,097,080 | 1 | . | A | 96.3% | 171.2 / ‑1.2 | 54 | intergenic (‑104/+46) | azoJ/yocK | FMN‑dependent NADH‑azoreductase/putative general stress protein |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base A (29/23); total (29/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
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. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
ATTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCTTTTAC > BS168_reference_genome/2096952‑2097212 | aTTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑a‑aaa‑aaaagc < 6:248136/139‑4 (MQ=255) ttttACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑a‑aaa‑aaaagc < 1:290572/138‑4 (MQ=255) aTAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATc < 6:147660/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACTAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGTATATGTGTTTTTTTTTCA‑A‑AAA‑AAAGCTAGTACATCTAatat < 1:323878/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAatat < 4:313598/139‑1 (MQ=255) cTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAatat < 2:110765/139‑1 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCgg > 8:172635/1‑139 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑a‑aaa‑aaaagcaa < 7:317143/123‑6 (MQ=255) tAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑a‑aaa‑aaaagcaa > 8:317143/1‑118 (MQ=255) aTTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGt < 1:257195/120‑1 (MQ=255) aTTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGt > 2:257195/1‑120 (MQ=255) tGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAA‑AAA‑AAAGCAAGTACATCTAATATCGGTGTActt > 1:248064/1‑139 (MQ=255) ttCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTActtgcttg > 6:356511/1‑139 (MQ=255) ttCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTActtgctt < 3:135984/138‑1 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 8:47154/1‑139 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 6:254421/1‑139 (MQ=255) ctccCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCa > 7:51225/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 4:239966/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtt > 6:355416/1‑139 (MQ=255) cccGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAAGATCGGTGTACTTGCTTGGCAtt > 4:383667/1‑139 (MQ=255) cGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt < 8:356643/139‑1 (MQ=255) ggTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAtttt > 1:287823/1‑138 (MQ=255) ttttATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGaa > 7:415876/1‑139 (MQ=255) tttaATATTTACTTACAAATAGTATGATACTTTCTTTTTTACGTAAAGTCAATCAATTTATTTTGAATATGAGATTATTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATAACGGTGTACTTGATTGGCATTTTTACTTCGa > 6:331121/1‑138 (MQ=255) ttacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCTTTTTTACTTCGAATCTTGTTc > 4:76966/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 4:287224/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 8:241379/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCa > 2:181788/1‑139 (MQ=255) tacttacAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGAATTTTTACTTCGAATCTTGTTCa > 4:195662/1‑139 (MQ=255) actttcAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑A‑AAATAAAGCAAGTACATCTAATATCGGTGTTCTTGATTTTCATTTTTACTTCGAATCTTGTTCa < 4:83344/138‑1 (MQ=255) tATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTAAATCTAATATCGGTGTACTTGCTTGGCTTTTTTACTTCGAATCTTGTTCACGATCCaatcttt > 8:204802/1‑134 (MQ=255) ttACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATa < 1:176109/139‑1 (MQ=255) ctttctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGt > 8:91364/1‑139 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAttt > 8:53934/1‑97 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCAttt < 7:53934/97‑1 (MQ=255) tctttTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTg > 2:285849/1‑139 (MQ=255) ttttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTAACGATCAAAGATTTCATACGTCTGCTAg > 7:28714/1‑139 (MQ=255) tttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAg < 1:409881/138‑1 (MQ=255) tttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGc > 7:269925/1‑139 (MQ=255) tttACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑ATAAA‑AAAGCTAGTACATCTACTATCGGTGTACTTTCTTGGAATTTTTACTTCGAATCTTGTTCAGGATCCAAGATTTCATACGTCTGCGAg < 5:331121/138‑1 (MQ=255) aaGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCAAA‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTg < 2:248064/139‑1 (MQ=255) aaGTCAATCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGa < 6:164610/139‑1 (MQ=255) aTCACTTTATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAAtt > 2:147770/1‑139 (MQ=255) tATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCaa > 3:323434/1‑130 (MQ=255) tATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCaa < 4:323434/130‑1 (MQ=255) tATTTTGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCATTTGTTTCtt > 1:153011/1‑139 (MQ=255) ttttGAATATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTttg < 1:181788/139‑1 (MQ=255) aaTATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCaa > 4:183602/1‑91 (MQ=255) aaTATGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCaa < 3:183602/91‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 7:47154/139‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 7:364983/139‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 5:356511/139‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 5:208227/139‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 8:51225/139‑1 (MQ=255) tatGAGATTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATa < 1:147770/139‑1 (MQ=255) gagaTTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTAGATCAATGGTTTCTTTGTTGGATAAGt > 1:83902/1‑139 (MQ=255) gagaTTTTTTTTCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCAGCGAGCTTTTTTGATAAATTGTTTCTTTGTTGGATAAg < 8:415876/138‑1 (MQ=255) ttttCA‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCTTTTAc < 7:172635/138‑1 (MQ=255) | ATTTTACAAATAAAACTGTAGACATTTGTATTCCTCCCGGTAAACTTTTTATATTTACTTACAAATAGTATGTTACTTTCTTTTTTACGTAAAGTCAATCACTTTATTTTGAATATGAGATTTTTTTTC‑‑A‑AAA‑AAAGCAAGTACATCTAATATCGGTGTACTTGCTTGGCATTTTTACTTCGAATCTTGTTCACGATCCAAGATTTCATACGTCTGCGAGCTTTTTTGATCAATTGTTTCTTTGTTGGATAAGTCTTTTAC > BS168_reference_genome/2096952‑2097212 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |