Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,919,423	C→A	9.8%	L404L (CTG→CTT)	galT ←	galactose‑1‑phosphate uridyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,919,423	0	C	A	9.8%	146.8 / 4.2	61	L404L (CTG→CTT)	galT	galactose‑1‑phosphate uridyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base A (5/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCCACTCACTATGCTTTGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGTCTCCCCTGCGCCGCGCAAT  >  BS168_reference_genome/3919287‑3919557                                                                                                                                         |                                                                                                                                        cccACTCACTATGCTTTGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAg                                                                                                                                       >  4:233987/1‑139 (MQ=255)           aTGCTTTGCAGTGAGCGGATTTTGTTCTATGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGGAAGCCGCCAGGCCATATCGCCTATCCCATGACCTGAATAAGGCCGA‑AAttt                                                                                                                            >  1:33572/1‑139 (MQ=255)           aTGCTTTGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGGTTCGTTGCTTTCCTCCTGTAATCCGCCCGGCAATATCCCCACACCCCAGCCCTCAAGCATGCCCA‑TCttt                                                                                                                            >  3:109877/1‑139 (MQ=255)                ttGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTGCTCCTGTAAGCGGCCCGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCttt                                                                                                                       >  5:94152/1‑139 (MQ=255)                ttGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCGGCCGTTTCGTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCCTGGCCGCAATCAGGCCGA‑TAttttagtt                                                                                                                       >  1:21421/1‑135 (MQ=255)                ttGCAGTGAGCGGA‑TTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCATTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCtttt                                                                                                                      >  1:287529/1‑139 (MQ=255)                 tGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCtttt                                                                                                                      >  8:312356/1‑139 (MQ=255)                  gCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTGTTTGATTTTCTCCGGTAATCGGGCCGGCCATATAGCCCAAACCCTGACCGCAAACACTGCGA‑GATTTTAttttt                                                                                                                     >  4:123511/1‑139 (MQ=255)                    aGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTGTTTCATTTCCTCCTGTAAGCGGCCAGGCAAAAACGCCAATACCATGTACTCAATCAGGCCGA‑TATTTTGtttttt                                                                                                                    >  1:250310/1‑138 (MQ=255)                    aGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCCA‑TATTTTCtttttt                                                                                                                    >  3:293995/1‑138 (MQ=255)                        aGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTGTTTAATTTACACCGGTAAGCCGCCCCGGAATATCGCCAAAACCCAGGCCCCCATAATGCCGA‑AATTTTTtttattggtca                                                                                                               >  8:354823/1‑134 (MQ=255)                        aGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTGTTTAAATTCCGCCTGTAAGCGGCCCTGCAATATCGCCAATACCATCACCTCAATCAGGCGGA‑TTTTTTCtgtttaata                                                                                                                 >  6:401352/1‑133 (MQ=255)                        aGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTGATTTGCTGCGTTAAGCGGCCCGGCAAAATAGCCAATACCAATACCTCAAACAGGCCCATTATTTTTTTTTTTaaaa                                                                                                                >  1:276812/1‑137 (MQ=255)                        aGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAatat                                                                                                               >  6:289731/1‑139 (MQ=255)                            gATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATAtggtg                                                                                                           <  3:123511/139‑1 (MQ=255)                                       aaGGCTGTTTTCGGATCAGCAGAACATAGCGCTGCTGCCGTGTCTTTCATTTCCTCCGGTAAGCGGCCCGGCAATATCGCCAAACCCATGGCCGCCATCAGGCCGA‑TATTTTATTTTTTAATTTTGTGGAGTTCTtgat                                                                                                >  6:124399/1‑139 (MQ=255)                                       aaGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTTAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATAAGGCCGA‑TTTTTTATTTTTTTATATTGTGGTATTGTtgag                                                                                                >  6:77004/1‑138 (MQ=255)                                        aGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTtgat                                                                                                <  7:312356/138‑1 (MQ=255)                                             gTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCAGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTGTTTTTTAATTTGGTTGACTTCTTGATGAGGAt                                                                                          >  8:77255/1‑139 (MQ=255)                                                 gcGGATCAGCATAACATAGCGCCGCTGCCGTTGCTTTCAATTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGaa                                                                                       <  4:100713/137‑1 (MQ=255)                                                      tCAGCAGAAAATAGCGGCGATGCCGTTTATTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTc                                                                                 <  2:285141/139‑1 (MQ=255)                                                              caattgcgcGGTTGCGGTTTCTTAAATTACATCGTTTACGCGGCCAGGCAATATCGCCAATCCCTGGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  2:396680/134‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGTCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCAGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  5:401352/139‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  8:236863/139‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  3:233987/139‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGGGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  1:226436/139‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTGCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  5:168553/139‑1 (MQ=255)                                                              aCATAGCGCCGCTGCCGTTTCTTTCATTTCCTCATGTAAGCGGCCAGGCAATATCGCCAATCCCTTGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGa                                                                         <  3:97444/139‑1 (MQ=255)                                                               cATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGAt                                                                        >  1:5011/1‑139 (MQ=255)                                                                  aGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGct                                                                     >  7:48050/1‑139 (MQ=255)                                                                   gcgcCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGctc                                                                    <  4:293995/139‑1 (MQ=255)                                                                       cgcTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATAGCCAATCCCATGACCTCAATCACGCCGA‑TATTTTCTTTTTTAAAATGGTGGACTTCATGATGAGGATGAAATATTCCTAATTGATGctccta                                                                 >  3:206195/1‑137 (MQ=255)                                                                         cTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGa                                                                                                         >  2:30975/1‑96 (MQ=255)                                                                         cTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGa                                                                                                         <  1:30975/96‑1 (MQ=255)                                                                         cTGCCGTTTCTTTCAATTCCTCCTGTTAGCGGCCAGGCAATATCGCCAATCCCATGCCCTCAATCAGGCCGA‑TATTTTGTTTTTTAATATTGTGGACTTCTTTATGAGGATGAAATATTGCTAATTGATGGTCctcctc                                                              >  8:226330/1‑139 (MQ=255)                                                                            ccGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCAtccgt                                                           <  2:250310/139‑1 (MQ=255)                                                                              atttctttcATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCAtccgtcc                                                         <  2:287529/138‑1 (MQ=255)                                                                                       aTTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTc                                                 <  2:5011/138‑1 (MQ=255)                                                                                        tttCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCt                                               <  1:31043/139‑1 (MQ=255)                                                                                        tttCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCt                                               <  5:124399/139‑1 (MQ=255)                                                                                        tttCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCt                                               <  5:108386/139‑1 (MQ=255)                                                                                            gtcctGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCAGA‑TATTTTCTTTTTTATTATGGTTGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGt                                           <  4:15425/138‑1 (MQ=255)                                                                                             tcctGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGAATTCTTGATGAGGATGAAATATTCATAATGGATGCTCCTCATCCGTCAGGTTGTTACTCAGt                                           >  2:25796/1‑138 (MQ=255)                                                                                               ctGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACg                                        <  1:64483/139‑1 (MQ=255)                                                                                                    agaggacAGTCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATc                                   <  1:325695/133‑1 (MQ=255)                                                                                                            ggCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGctcct                                                                  <  1:260280/100‑1 (MQ=255)                                                                                                            ggCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGctcct                                                                  >  2:260280/1‑100 (MQ=255)                                                                                                             gCAATATCGCCAATCCCATGACCCCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTa                          <  4:245766/139‑1 (MQ=255)                                                                                                                  atCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGt                     <  8:48050/139‑1 (MQ=255)                                                                                                                   tctccAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGt                     <  4:206195/135‑1 (MQ=255)                                                                                                                   tccccAATCCCATTACCACAATCAGGCCGA‑TTTTTTCTTTTTTAATATTGTGGACTTCATGAAGATGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGt                     <  1:116891/135‑1 (MQ=255)                                                                                                                   tCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGt                     <  6:94152/138‑1 (MQ=255)                                                                                                                   tCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGt                     <  7:77255/138‑1 (MQ=255)                                                                                                                      ccAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGTCTcc                 >  7:282168/1‑139 (MQ=255)                                                                                                                         aTCCCATGACCTCAATAAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTc                                                                                 >  6:277120/1‑72 (MQ=255)                                                                                                                         aTCCCATGACCTCAATAAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTc                                                                                 >  6:282808/1‑72 (MQ=255)                                                                                                                         aTCCCATGACCTCAATAAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTc                                                                                 <  5:282808/72‑1 (MQ=255)                                                                                                                         aTCCCATGACCTCAATAAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTc                                                                                 <  5:277120/72‑1 (MQ=255)                                                                                                                              aTGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGa                                 <  4:254008/115‑1 (MQ=255)                                                                                                                              aTGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGa                                 >  3:254008/1‑115 (MQ=255)                                                                                                                                      aaTCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGTCTCCCCTGCGCCGCGCAAt  <  2:276812/138‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CCCACTCACTATGCTTTGCAGTGAGCGGATTTTGTTCTAAGGCTGTTTTCGGATCAGCAGAACATAGCGCCGCTGCCGTTTCTTTCATTTCCTCCTGTAAGCGGCCAGGCAATATCGCCAATCCCATGACCTCAATCAGGCCGA‑TATTTTCTTTTTTAATATGGTGGACTTCTTGATGAGGATGAAATATTCCTAATGGATGCTCCTCATCCGTCCGGTTGTTCCTCAGTACGATATCGAGTTCGTACAAGTCTCCCCTGCGCCGCGCAAT  >  BS168_reference_genome/3919287‑3919557

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