Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,770,058:1	+A	100%	intergenic (‑250/+46)	ywzF ← / ← csbD	conserved protein of unknown function/stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,770,058	1	.	A	100.0%	91.6 / NA	29	intergenic (‑250/+46)	ywzF/csbD	conserved protein of unknown function/stress response protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (19/10);  total (19/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTG  >  BS168_reference_genome/3769926‑3770190                                                                                                                                      |                                                                                                                                     gTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑aaaaa                                                                                                                                 >  2:50634/1‑138 (MQ=255) gTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCATATTAATTAAAAGAATAATTAAAAAATTTAGTTACATTGACCTAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑aaaaa                                                                                                                                 >  5:178629/1‑138 (MQ=255)            ttACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGa                                                                                                                       >  1:68263/1‑138 (MQ=255)                 aCTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAtt                                                                                                                 <  6:59281/139‑1 (MQ=255)                     aTGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTTTTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGt                                                                                                             >  3:151403/1‑139 (MQ=255)                                      aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAtt                                                                                                                 <  4:192935/118‑1 (MQ=255)                                      aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAtt                                                                                                                 >  3:192935/1‑118 (MQ=255)                                      aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                            >  1:36186/1‑139 (MQ=255)                                      aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                            >  3:236404/1‑139 (MQ=255)                                      aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTACTTTTCCTGCCTCTAACGTa                                                                                            >  8:78719/1‑139 (MQ=255)                                           gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  5:142841/1‑138 (MQ=255)                                           gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  1:128862/1‑138 (MQ=255)                                           gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  7:252802/1‑138 (MQ=255)                                           gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  1:191803/1‑138 (MQ=255)                                           gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  5:295330/1‑138 (MQ=255)                                                ttAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCt                                                                                  <  8:267247/139‑1 (MQ=255)                                                     tCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATc                                                                             >  4:232700/1‑139 (MQ=255)                                                            atatTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTCTTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAAt                                                                      <  7:120175/139‑1 (MQ=255)                                                                        tttAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCt                                                          <  5:89486/139‑1 (MQ=255)                                                                                 cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                                 >  4:247631/1‑139 (MQ=255)                                                                                 cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                                 <  2:36186/139‑1 (MQ=255)                                                                                          aTCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTGTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTg                                         >  8:145883/1‑138 (MQ=255)                                                                                               cttcttTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTAGTTTCCTGCCTCTTCCGTAAGATTAATCTTTATATGAAAATTTATCCTTGGCTTTCCCTATATCATGTTTGATTTc                                   >  4:11197/1‑139 (MQ=255)                                                                                                                    ttCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCttttt               <  2:191803/138‑1 (MQ=255)                                                                                                                     tCACAACAGAAAAAGCGAAAAAAAAAGACAAGAAACAATTCGTTTCCTGCCTCTTCCGTACGATTAATATTTATCTGAAAATTTATCCTTGGCTTTCCAGATATCAATATGGATATCGCCTATTGCCTTAACttattt              >  4:86194/1‑138 (MQ=255)                                                                                                                        caTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCt          >  4:289492/1‑139 (MQ=255)                                                                                                                        caTCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCt          >  2:3529/1‑139 (MQ=255)                                                                                                                          tCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGAttt                                    <  1:260511/111‑1 (MQ=255)                                                                                                                          tCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGAttt                                    >  2:260511/1‑111 (MQ=255)                                                                                                                           cAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCa       <  6:295330/139‑1 (MQ=255)                                                                                                                                aaaGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTg  <  2:68263/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     GTTATAAAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTG  >  BS168_reference_genome/3769926‑3770190

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: ‑