Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,405,231	A→C	18.1%	S133R (AGC→CGC)	yvrD →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,405,231	0	A	C	18.1%	94.2 / 2.9	52	S133R (AGC→CGC)	yvrD	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (27/14);  new base C (3/7);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 7.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTGGTTTCTTCGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACAAAAGGCACAAATGTG  >  BS168_reference_genome/3405101‑3405365                                                                                                                                   |                                                                                                                                       ccTTGGTTTCTTCGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCg                                                                                                                                 >  7:159057/1‑138 (MQ=255)     ggTTTCTTCGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCAttt                                                                                                                            >  7:6433/1‑139 (MQ=255)      gTTTCTTCGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCAttt                                                                                                                            >  8:336622/1‑138 (MQ=255)          cttcGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTaaa                                                                                                                        >  5:233899/1‑138 (MQ=255)           ttcGAGGTAGAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAata                                                                                                                      >  6:342629/1‑139 (MQ=255)           ttcGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAata                                                                                                                      >  2:36272/1‑139 (MQ=255)           ttcGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAata                                                                                                                      >  6:231621/1‑139 (MQ=255)            tcGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAatat                                                                                                                     >  1:199894/1‑139 (MQ=255)                         cTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCa                                                                                                        <  7:336622/139‑1 (MQ=255)                            cGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGc                                                                                                     >  6:368687/1‑139 (MQ=255)                                      aCAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCg                                                                                            <  4:405050/138‑1 (MQ=255)                                          aTGAATTGTGTAATAAATATTTCGATGTGAACGTGATGAGCGCCGTTCGGACAAGCGGTCATTTCCTTCCAAAATTGCTGGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGcc                                                                                       <  8:6433/139‑1 (MQ=255)                                            gAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCaa                                                                                     >  3:58281/1‑139 (MQ=255)                                               gattgtaaTCAATATTTCGAAGTGAACCTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTGGCCAAAAACCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACga                                                                                  <  7:164180/133‑1 (MQ=255)                                               gattGGAATCAATATTTCTAAGTTAAGGTGATCAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACga                                                                                  <  5:206777/136‑1 (MQ=255)                                               gagTGGAATCAATATTTCGAAGTGAACGTGATTAGCGCCGTTCGGACAAGCCGGCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACTGCTGCCAACga                                                                                  <  7:64771/139‑1 (MQ=255)                                               gagTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACga                                                                                  <  7:406767/139‑1 (MQ=255)                                                agTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAACCCGCTGCCAACgat                                                                                 >  6:309846/1‑139 (MQ=255)                                                 gTGGAATCAATATTGCGAAGTGAACGGGATGCGCGCCGTTCGGACAAGCCGTCACTTCCTTCAAAAAATGCTTGCCAAAAACCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACgatg                                                                                <  3:192222/139‑1 (MQ=255)                                                  tGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACgatg                                                                                <  1:36272/138‑1 (MQ=255)                                                    gAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGAtt                                                                             >  1:204599/1‑139 (MQ=255)                                                       tCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCa                                                                          <  8:159057/139‑1 (MQ=255)                                                         tatttttCGAAGTGAACGTGATGAGCGCGGTTCGTAAAAGCCGTCACTTCCTTCAAAAAATGCTTGCAAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATa                                                                        <  6:142282/135‑1 (MQ=255)                                                         ccgatTTCGACGTGAAGGTGCTGCGCGCCGTTGGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAAGCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATa                                                                        <  6:35168/136‑1 (MQ=255)                                                         aaTATTTGCATGTGAACGTGATGAGCGCCGTGGGGACACGCGGGAAATTCATTCCAAAAATGCTTGCCAAAAAACGCGGCCTCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATa                                                                        <  3:163751/139‑1 (MQ=255)                                                         aaTATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATa                                                                        <  6:285637/139‑1 (MQ=255)                                                          atatTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATAc                                                                       >  6:289335/1‑139 (MQ=255)                                                          atatTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATACCATa                                                                        >  1:309021/1‑138 (MQ=255)                                                             tttCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCa                                                                    <  7:249269/139‑1 (MQ=255)                                                                cGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCTACGATGATTCCATACTCAAt                                                                  >  5:376464/1‑138 (MQ=255)                                                                cGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAAt                                                                  >  1:352526/1‑138 (MQ=255)                                                                     tGAACGTGATGAGCGCCGTTCGGACAAGCCGTCCCTTCCTTCCAAAAATGCTTGCCAAAAACCGCGGCCGCATTTTTAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCCACGATGATTCCATGCTCAATGACaa                                                             >  4:29307/1‑138 (MQ=255)                                                                     tGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACaac                                                            >  5:142282/1‑138 (MQ=255)                                                                       aaCGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACCGCGGCCGCAGTATAAATATTGCAAGTGAAGCAGGCGTCACACCGCTGCCCACGATGATTCCATACTCAAGGACaactg                                                          >  5:35168/1‑136 (MQ=255)                                                                       aaCGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAATCCGCTGCCAACGATGATTCCATACTCAATGCCaaaaa                                                          >  1:413609/1‑139 (MQ=255)                                                                                  cgcgGTTCGGCCAAGCCGTCACTTCCATCAAAAAATGCTTCCAAAAAACCGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATTATTCCATACTCAATGACAAAAACGGCACTGATc                                               <  3:51435/139‑1 (MQ=255)                                                                                  cgcCGTGCGGACAACCCGGCACTTCCTGCCAAAAAGGCTTCCCAAAACCCGCGGCCGCATTTTAAATATTGCAAGTGACGCAGGCGTCAACCCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATc                                               <  8:91170/139‑1 (MQ=255)                                                                                    ccGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAg                                             >  7:147761/1‑139 (MQ=255)                                                                                    ccGTCCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAg                                             <  3:181272/139‑1 (MQ=255)                                                                                       ttCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCt                                          >  3:184284/1‑139 (MQ=255)                                                                                       ttCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCt                                          >  7:77818/1‑139 (MQ=255)                                                                                           caCAAGCCGCAACTACCTACAAAAAATGCTTCCAAAAAACAGCGGCCGAATTTTAAATATTGCAAGTTAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCATTTACAAAAACGGCACTGATCAGCCTTTCa                                      <  3:183184/138‑1 (MQ=255)                                                                                               aGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAACCCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAagag                                  >  3:315618/1‑139 (MQ=255)                                                                                               aGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAagag                                  >  7:204978/1‑139 (MQ=255)                                                                                               aGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAagag                                  >  3:93869/1‑139 (MQ=255)                                                                                                   gTCAATACCTACAAAAACTGCTTCAAAAAACCTGCGGCCGAATTTTAAATATTGCAATTGAAGCCGGCGTCAACCCGCGGCCAACGATGATACCATACTCAATTAAAAAAACGGCACTGATCAGCCTTTCAAGAGGAAt                              <  2:352526/139‑1 (MQ=255)                                                                                                              ccAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGa                   >  7:126982/1‑139 (MQ=255)                                                                                                              ccAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGa                   >  1:323391/1‑139 (MQ=255)                                                                                                                  aaaGGCTTGCCAAACCCAGCGGCCGAATTTAAAATATTGCAAGTGACGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACaaa               <  6:103266/139‑1 (MQ=255)                                                                                                                      gCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACAAAAgg            >  3:50593/1‑138 (MQ=255)                                                                                                                         tGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACAAAAGGcac         >  5:140989/1‑138 (MQ=255)                                                                                                                               aaaCAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACAAAAGGCACAAAtgtg  <  1:41221/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       CCTTGGTTTCTTCGAGGTAAAAGACTTCGCCGATGTGACAGATGAAGAGTGGAATCAATATTTCGAAGTGAACGTGATGAGCGCCGTTCGGACAAGCCGTCACTTCCTTCCAAAAATGCTTGCCAAAAACAGCGGCCGCATTTTAAATATTGCAAGTGAAGCAGGCGTCAAACCGCTGCCAACGATGATTCCATACTCAATGACAAAAACGGCACTGATCAGCCTTTCAAGAGGAATGGCTGAAATGACAAAAGGCACAAATGTG  >  BS168_reference_genome/3405101‑3405365

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