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%	168.8 / NA	49	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 (28/21);  total (28/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTG  >  BS168_reference_genome/3769932‑3770190                                                                                                                                |                                                                                                                                     aaaGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑aaaaaaaaagac                                                                                                                          <  2:169088/139‑5 (MQ=255)    gTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAg                                                                                                                        <  8:35257/139‑1 (MQ=255)      ttACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCATCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGa                                                                                                                       >  1:94472/1‑138 (MQ=255)       tACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGaaa                                                                                                                     >  7:157232/1‑139 (MQ=255)               aTGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGt                                                                                                             >  8:184987/1‑139 (MQ=255)               aTGTCATCATTATATAGAAGAGGCCTATTAATTAAAAGAATACTAAAAAAATTAAGTTACCTTCACCAAATAAGCTACTTTATAAGTTATATATATTCACATCAGAAACAGAGAAAAAAAAAGACAAGAAACATTTCga                                                                                                             >  3:79449/1‑138 (MQ=255)                    atcatTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCt                                                                                                        >  1:143277/1‑139 (MQ=255)                      catTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGc                                                                                                      >  4:182473/1‑139 (MQ=255)                          atataGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCctc                                                                                                   <  6:29861/138‑1 (MQ=255)                               gaagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCg                                                                                              <  6:176232/138‑1 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAtt                                                                                                                 >  5:196393/1‑118 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAtt                                                                                                                 <  6:196393/118‑1 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                           >  1:12061/1‑124 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                           <  2:12061/124‑1 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                            >  7:39714/1‑139 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                            >  7:91909/1‑139 (MQ=255)                                aagaGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTa                                                                                            >  6:219768/1‑139 (MQ=255)                                 agagGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTTCGTTGACGAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCATATAAAGCGAAAAAAAAAGACTTGTTACATTTCGTTTCCTTCCTCTTCCGTa                                                                                            <  4:73972/138‑1 (MQ=255)                                   agGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACAGAAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGa                                                                                         >  6:45196/1‑139 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                           <  2:110495/119‑1 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  2:72347/1‑138 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  6:142310/1‑138 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  3:57903/1‑138 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGAt                                                                                        >  1:167157/1‑138 (MQ=255)                                     gCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACATTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGttt                                                                                                           >  1:110495/1‑119 (MQ=255)                                          ttAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAttt                                                                                                                >  5:58395/1‑109 (MQ=255)                                          ttAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACAttt                                                                                                                <  6:58395/109‑1 (MQ=255)                                               tCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATc                                                                             >  6:18228/1‑139 (MQ=255)                                                      atatTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAAt                                                                      <  3:182473/139‑1 (MQ=255)                                                      atatTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAAt                                                                      <  2:94472/139‑1 (MQ=255)                                                                 aTTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGc                                                           <  8:91909/139‑1 (MQ=255)                                                                  tttAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCGTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCATTGGCt                                                          >  3:191338/1‑139 (MQ=255)                                                                  tttAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCt                                                          <  5:18228/139‑1 (MQ=255)                                                                          ccTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGat                                                  <  4:215619/139‑1 (MQ=255)                                                                          ccTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGa                                                   <  5:219768/138‑1 (MQ=255)                                                                           cTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGata                                                 <  1:72347/139‑1 (MQ=255)                                                                               aCCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCt                                             >  1:144611/1‑139 (MQ=255)                                                                                    aTCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTg                                         >  3:212327/1‑138 (MQ=255)                                                                                         cttcttTTTTTGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGAtt                                     >  2:48115/1‑137 (MQ=255)                                                                                              ttttttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCtt                              >  7:17040/1‑139 (MQ=255)                                                                                               tttttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCtt                              <  7:184987/138‑1 (MQ=255)                                                                                               tttttGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCtt                              <  1:48115/138‑1 (MQ=255)                                                                                                   tGTTATATTTATTCACATCAGAAAAAGCGAAAAAAAAAGACAATAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGcc                         >  5:37846/1‑139 (MQ=255)                                                                                                              ttCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCttttt               <  5:44192/138‑1 (MQ=255)                                                                                                               tCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATAAATCTTTATCTGAAAATTTATCATTGGCTTTCACGAGATCATTTTGGATGTCGCCTTTTGCCTTATCTTTTTTg             >  3:177860/1‑139 (MQ=255)                                                                                                               tCACATCAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCGGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTg             >  8:16062/1‑139 (MQ=255)                                                                                                                     cAGAAAAAGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCa       <  2:61072/139‑1 (MQ=255)                                                                                                                        aaaaaGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCtt         >  5:119290/1‑134 (MQ=255)                                                                                                                        aaaaaGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCtt         <  6:119290/134‑1 (MQ=255)                                                                                                                          aaaGCGAAAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTg  <  4:57903/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                     AAAGTTTACAACTTATGTCATCATTATATAGAAGAGGCCTATTAATTCAAAGAATATTTAAAAAATTTAGTTACCTTGACCAAATCAGCTTCTTTTTTTGTTATATTTATTCACATCAGAAAAAGCG‑AAAAAAAAGACAAGAAACATTTCGTTTCCTGCCTCTTCCGTACGATTAATCTTTATCTGAAAATTTATCCTTGGCTTTCCCGATATCCTTTTGGATTTCGCCTTTTGCCTTATCTTTTTTGCCTTCAGCCTG  >  BS168_reference_genome/3769932‑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: ‑