Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,306,179	G→A	6.9%	G268S (GGC→AGC)	metAA →	homoserine O‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,306,179	0	G	A	6.9%	182.1 / 3.3	58	G268S (GGC→AGC)	metAA	homoserine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (2/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TCTAAAGACGAAAAACAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACGCTGAGCAAA  >  BS168_reference_genome/2306047‑2306314                                                                                                                                     |                                                                                                                                        tCTAAAGACGAAAAACAAGTCTTTTTGACAGGGCAGCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGgcagc                                                                                                                                    >  4:131914/1‑138 (MQ=255) tCTAAAGACGAAAAACAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGgcagc                                                                                                                                    >  4:127423/1‑138 (MQ=255) tCTAAAGACGAAAAACAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGgcagc                                                                                                                                    >  7:218056/1‑138 (MQ=255)     aaGACGAAAAACAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACg                                                                                                                                <  1:283630/138‑1 (MQ=255)           aaaaaCAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGt                                                                                                                         >  5:216055/1‑139 (MQ=255)                    tCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTg                                                                                                                >  1:213636/1‑139 (MQ=255)                             cAGCGCACCCTGAATATGAAACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCATTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCa                                                                                                       <  5:271764/139‑1 (MQ=255)                               gggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATg                                                                                                     <  1:50529/139‑1 (MQ=255)                               gggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATg                                                                                                     <  7:194472/139‑1 (MQ=255)                               gggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATg                                                                                                     <  5:245347/139‑1 (MQ=255)                               gggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACAATATTTTGCAAAAGGCAGCAACGAGCCAGGGAATCGCTGGAAAGCGCATg                                                                                                     <  1:29678/139‑1 (MQ=255)                                ggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGc                                                                                                    >  1:194220/1‑139 (MQ=255)                                ggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGc                                                                                                    >  3:123259/1‑139 (MQ=255)                                ggCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGc                                                                                                    >  6:134178/1‑139 (MQ=255)                                               aTACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTAt                                                                                     >  2:89209/1‑139 (MQ=255)                                               aTACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTAt                                                                                     >  3:135320/1‑139 (MQ=255)                                                   aaaTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGaa                                                                                  <  7:56156/138‑1 (MQ=255)                                                   aaaTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGaa                                                                                  >  8:56156/1‑138 (MQ=255)                                                   aaaTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAAt                                                                                 >  1:73879/1‑139 (MQ=255)                                                   aaaTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAAt                                                                                 >  7:136239/1‑139 (MQ=255)                                                   aaaTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAAt                                                                                 >  3:102998/1‑139 (MQ=255)                                                     aTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACatta                                                                                             >  1:71364/1‑125 (MQ=255)                                                     aTACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACatta                                                                                             <  2:71364/125‑1 (MQ=255)                                                      tACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTgg                                                                              <  3:141740/139‑1 (MQ=255)                                                        cGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCt                                                                            <  1:172807/139‑1 (MQ=255)                                                         gCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTAGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCt                                                                            >  2:7145/1‑138 (MQ=255)                                                          cttcttCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTg                                                                           >  8:42849/1‑138 (MQ=255)                                                             cttcAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAat                                                                        >  5:83505/1‑138 (MQ=255)                                                                     gTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGtt                                                               >  4:98525/1‑139 (MQ=255)                                                                     gTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAAGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGtt                                                               >  5:60686/1‑139 (MQ=255)                                                                     gTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAAGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGtt                                                               >  5:56172/1‑139 (MQ=255)                                                                        tGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTa                                                             <  4:18369/138‑1 (MQ=255)                                                                                   tGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTc                                                  <  5:97899/138‑1 (MQ=255)                                                                                          aaaCCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGaa                                          >  8:232674/1‑139 (MQ=255)                                                                                          aaaCCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGaa                                          <  5:91498/139‑1 (MQ=255)                                                                                          aaaCCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGaa                                          >  5:129557/1‑139 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     <  2:47143/114‑1 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     <  4:886/114‑1 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     >  1:47143/1‑114 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     >  3:886/1‑114 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAAttttt                            >  3:73186/1‑139 (MQ=255)                                                                                                        tCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAAttttt                            >  6:76407/1‑139 (MQ=255)                                                                                                          gAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCt                          >  6:8966/1‑139 (MQ=255)                                                                                                               cgcTAAACTTTATTTTTCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGTCTTAATTATTACGTATACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtttttt                     <  3:127423/137‑1 (MQ=255)                                                                                                               cccTTAACATTAATTTTCCAAAGGCAGCAACGAGCCAGTGCATCTCTGTAAACCGCATGCCACATTAATGTTAATGAAATTGCTGAATTATTAAGTATACCAAGAAACTACTTATGAATTTGACTAAATTTTAattttt                     <  4:120949/139‑1 (MQ=255)                                                                                                               cccTAATAATTATTTTTCAAAAGGCAGCAACGATCCAGTGAATCTCTTGAAATCGAATGCCACATTATTGTTTATTAATTTTCTTAATTATTACGTTTACCAAGAAACTCCTTATTAATGGGACTAATTTTTCtttttt                     <  3:131914/139‑1 (MQ=255)                                                                                                               cccTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGACTTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtttttt                     <  4:222708/139‑1 (MQ=255)                                                                                                               cccTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtttttt                     <  1:89209/139‑1 (MQ=255)                                                                                                               cccTAAACATTATTTTGCAAAAAGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtttttt                     <  6:60686/139‑1 (MQ=255)                                                                                                               cccTAAACATTATTTTGCAAAAAGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtttttt                     <  6:56172/139‑1 (MQ=255)                                                                                                                         tattTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     >  6:232310/1‑97 (MQ=255)                                                                                                                         tattTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGaaa                                                     <  5:232310/97‑1 (MQ=255)                                                                                                                         tattTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACg           >  7:150670/1‑139 (MQ=255)                                                                                                                           ttttGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACGCt         >  7:18848/1‑139 (MQ=255)                                                                                                                           ttttGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACGCt         >  3:148799/1‑139 (MQ=255)                                                                                                                              tGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtt                         <  6:30016/120‑1 (MQ=255)                                                                                                                              tGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCtt                         >  5:30016/1‑120 (MQ=255)                                                                                                                                  aaaGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACGCTGAGCaaa  <  7:42849/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        TCTAAAGACGAAAAACAAGTCTTTTTGACAGGGCACCCTGAATATGATACAAATACGCTTCTTCAAGAGTATGAAAGAGATTTGGAACGAAACCTTTCTACTGTCGAGGCCCCTAAACATTATTTTGCAAAAGGCAGCAACGAGCCAGTGAATCGCTGGAAAGCGCATGCCACATTATTGTTTATGAATTGGCTGAATTATTACGTTTACCAAGAAACTCCTTATGAATGGGACTAATTTTTCTTTTTTTAGTGATACGCTGAGCAAA  >  BS168_reference_genome/2306047‑2306314

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