Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	4,068,272	T→A	6.9%	K89* (AAG→TAG)	yxeA ←	conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	4,068,272	0	T	A	6.9%	173.8 / 4.5	58	K89* (AAG→TAG)	yxeA	conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (22/32);  new base A (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

GAAAATAAAAAATAAAAAAACGGCACAGTCATGACGCTGTGCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCATCTCGATCAAT  >  BS168_reference_genome/4068134‑4068407                                                                                                                                           |                                                                                                                                        caaaaaaaaaaaaaaaaaaACGGCACAGTCATGACGCTGTGCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGcctt                                                                                                                                         <  7:97018/138‑1 (MQ=255) aaaaaaaaaaaataaAAACACGGCACAGTCATGACGCTGTGCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGcctt                                                                                                                                         <  3:185466/138‑1 (MQ=255)                             tCATGACGCTGTGCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTGGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGc                                                                                                             >  7:142712/1‑139 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGCCTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   >  3:75516/1‑139 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   >  3:279037/1‑139 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   >  8:42615/1‑139 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   >  6:63801/1‑139 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   <  5:100469/139‑1 (MQ=255)                                       gtgCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCgg                                                                                                   >  7:47394/1‑139 (MQ=255)                                             ttttttATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGaaa                                                                                              >  1:30805/1‑138 (MQ=255)                                             ttttttATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGaaa                                                                                              >  8:100952/1‑138 (MQ=255)                                                  tATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGt                                                                                         <  1:286116/138‑1 (MQ=255)                                                                     acaggTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTACACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCgg                                                                     <  6:21462/135‑1 (MQ=255)                                                                     acaggTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTACACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCgg                                                                     <  6:21463/135‑1 (MQ=255)                                                                      aCTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCg                                                                      >  8:120437/1‑137 (MQ=255)                                                                      aCTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  7:42615/139‑1 (MQ=255)                                                                      aCTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  1:63146/139‑1 (MQ=255)                                                                      aCTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  1:276769/139‑1 (MQ=255)                                                                      aCTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  7:100952/139‑1 (MQ=255)                                                                               aGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATtata                                                           <  4:83899/139‑1 (MQ=255)                                                                               aGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATtata                                                           <  1:272737/139‑1 (MQ=255)                                                                               aGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATtata                                                           >  3:216531/1‑139 (MQ=255)                                                                               aGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATtata                                                           >  2:16822/1‑139 (MQ=255)                                                                               aGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATtata                                                           <  1:65309/139‑1 (MQ=255)                                                                                gTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          <  2:299494/139‑1 (MQ=255)                                                                                gTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          <  6:206141/139‑1 (MQ=255)                                                                                    ggCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          >  6:8218/1‑135 (MQ=255)                                                                                    ggCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          >  6:274401/1‑135 (MQ=255)                                                                                    ggCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          <  5:274401/135‑1 (MQ=255)                                                                                    ggCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTatat                                                          <  5:8218/135‑1 (MQ=255)                                                                                    ggCATATCTTCAAACTTGACTTCTTAAAAAGTCTCGACATACTTTCCTTTCGCCTACACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCTGt                                                      >  5:21462/1‑139 (MQ=255)                                                                                    ggCATATCTTCAAACTTGAATTCTTCAAAAGTCTCGACATACTTTCCTTTCGCCTACACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCTGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCCTTATATCTGt                                                      >  5:21463/1‑139 (MQ=255)                                                                                      cATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATATGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCt                                                    <  5:98228/139‑1 (MQ=255)                                                                                      cATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCt                                                    <  1:16822/139‑1 (MQ=255)                                                                                      cATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCt                                                    <  8:41349/139‑1 (MQ=255)                                                                                      cATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCt                                                    <  8:47394/139‑1 (MQ=255)                                                                                          tCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAAt                                                 <  4:216531/138‑1 (MQ=255)                                                                                              cAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGt                                               <  7:121928/136‑1 (MQ=255)                                                                                               aaaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGCAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTatat                                           >  5:132179/1‑139 (MQ=255)                                                                                               aaaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAAt                                                 >  5:233210/1‑133 (MQ=255)                                                                                               aaaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAAt                                                 <  6:233210/133‑1 (MQ=255)                                                                                               aaaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAAt                                                 <  4:121529/133‑1 (MQ=255)                                                                                               aaaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAAt                                                 >  3:121529/1‑133 (MQ=255)                                                                                                aaCTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATAtt                                          <  8:142712/139‑1 (MQ=255)                                                                                                      aCTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCg                                                                                                                 <  4:154485/62‑1 (MQ=255)                                                                                                      aCTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCg                                                                                                                 >  3:154485/1‑62 (MQ=255)                                                                                                           ttcCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGttcttc                                                                                                                   <  8:94167/55‑1 (MQ=255)                                                                                                           ttcCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGttcttc                                                                                                                   >  7:94167/1‑55 (MQ=255)                                                                                                             cccAAGTCTCGACATACTTTCCTTTAGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGcc                                <  8:151137/136‑1 (MQ=255)                                                                                                              ccAAGTCTAGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCgg                             <  7:310358/138‑1 (MQ=255)                                                                                                                             ctttcctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGcc              <  4:279037/138‑1 (MQ=255)                                                                                                                               ttcctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCAtc           >  7:136194/1‑139 (MQ=255)                                                                                                                                 cctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    >  7:238280/1‑80 (MQ=255)                                                                                                                                 cctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  8:238280/80‑1 (MQ=255)                                                                                                                                 cctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGt                                               >  2:98680/1‑101 (MQ=255)                                                                                                                                 cctttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGt                                               <  1:98680/101‑1 (MQ=255)                                                                                                                                    ttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    >  2:294204/1‑77 (MQ=255)                                                                                                                                    ttcGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGa                                                                    <  1:294204/77‑1 (MQ=255)                                                                                                                                     tcGCCATCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCATCTCGa       <  5:75038/137‑1 (MQ=255)                                                                                                                                      cGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCATCTCGATc     <  3:143402/138‑1 (MQ=255)                                                                                                                                        ccttcaccttcaAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCATCTCGATCAAt  <  5:63801/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        GAAAATAAAAAATAAAAAAACGGCACAGTCATGACGCTGTGCCGTTTTTTATGATTCACTTTAATTTTGACTGTACGGAGTCCGGCATATCTTCAAACTTGACTTCTTCCCAAGTCTCGACATACTTTCCTTTCGCCTTCACCTTCAAATACGCGTTCTTCCGAAGCTCTTTGCCGGCGAAAAAGGTCACTTCTTCTTTTTTCCCGGATGCATTATATCCGTCTAATGTATATTCAGTACCGCCCGGACTAAGATGTCTGCCATCTCGATCAAT  >  BS168_reference_genome/4068134‑4068407

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑