Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,358,986	T→C	7.0%	E113E (GAA→GAG)	mgsA ←	methylglyoxal synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,358,986	0	T	C	7.0%	174.7 / 3.9	57	E113E (GAA→GAG)	mgsA	methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base C (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTGCATTGG  >  BS168_reference_genome/2358852‑2359119                                                                                                                                       |                                                                                                                                      tgCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCg                                                                                                                                   <  4:191046/139‑1 (MQ=255) tgCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCg                                                                                                                                   <  2:346370/139‑1 (MQ=255)   ccgccCATTCCTATCTCGACATCATCACTGTGGGCGACAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGGCACTTAATTCTAAAACACCTTAATCAAGTGTGAGCACAATAATTTCCGCa                                                                                                                                 <  5:188058/139‑1 (MQ=255)   ccgccCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                                                                 <  2:247030/139‑1 (MQ=255)   ccgccCATTCCGATCTCGACATCATCAATGTGGGCGACAAAATCAAGAACGTAATCATTATACATTCGGCTATGCTCCCCGAAGAATGTCACGGAATTAGAAAACACCTTACTCAAGTGTGAGCCCAAGAATTTcctcc                                                                                                                                 >  6:188058/1‑136 (MQ=255)                   cGACATCATCACTATGGGCGCCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTtgtag                                                                                                                 >  3:366494/1‑137 (MQ=255)                         catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                                                                 >  1:861/1‑117 (MQ=255)                         catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                                                                 <  2:861/117‑1 (MQ=255)                         catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                                                           >  7:192415/1‑139 (MQ=255)                          atcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCgg                                                                                                          <  3:183801/139‑1 (MQ=255)                          atcaCTATGGGCGCCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTTGTagcgagcg                                                                                                           >  1:262381/1‑138 (MQ=255)                           tcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCgga                                                                                                         <  7:205643/139‑1 (MQ=255)                            caCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                                                        >  1:127960/1‑139 (MQ=255)                                 tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaatggaa                                                                                                   >  6:424544/1‑139 (MQ=255)                                 tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                                                           >  5:129853/1‑131 (MQ=255)                                 tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                                                           <  6:129853/131‑1 (MQ=255)                                     cgcCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTTGTAGCGAGCAGAATGGAATaca                                                                                               <  4:366494/139‑1 (MQ=255)                                     cgcCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATaca                                                                                               <  1:51056/139‑1 (MQ=255)                                                aGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                                                                    <  3:210037/139‑1 (MQ=255)                                                       cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTAcc                                                                                                                            <  8:206598/92‑1 (MQ=255)                                                       cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTAcc                                                                                                                            >  7:206598/1‑92 (MQ=255)                                                       cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                                                                    <  1:91116/132‑1 (MQ=255)                                                       cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                                                                    >  2:91116/1‑132 (MQ=255)                                                               taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                                                        >  4:387335/1‑104 (MQ=255)                                                               taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                                                        <  3:387335/104‑1 (MQ=255)                                                               taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGcc                                                                        >  5:136025/1‑136 (MQ=255)                                                               taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGcc                                                                        >  6:392216/1‑136 (MQ=255)                                                                  aTTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACa                                                                  >  5:199240/1‑139 (MQ=255)                                                                    tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                                                        >  1:227416/1‑99 (MQ=255)                                                                    tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                                                        <  2:227416/99‑1 (MQ=255)                                                                    tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                                                                >  8:246852/1‑139 (MQ=255)                                                                    tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                                                                >  4:81931/1‑139 (MQ=255)                                                                    tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                                                                >  3:384594/1‑139 (MQ=255)                                                                    tCGGCTCTTCTCACCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCTCACAAGAATTTCCGCATTACCCATATTTGTTGCGAGCGGAATGGATTACACATCACATAAACTGATTAATGCCTAGACagc                                                                >  5:394246/1‑137 (MQ=255)                                                                      ggCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCg                                                              <  6:199240/139‑1 (MQ=255)                                                                       gCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCgg                                                             >  8:324862/1‑139 (MQ=255)                                                                                     agaagGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTgcg                                               <  5:430077/139‑1 (MQ=255)                                                                                          gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAACCGGATTAATGCCGAGACATCCGGTGCCTGCGGCTGCGCGGTc                                          >  2:303401/1‑139 (MQ=255)                                                                                          gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCgg                                            >  1:407054/1‑137 (MQ=255)                                                                                          gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCgg                                            <  2:407054/137‑1 (MQ=255)                                                                                          gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTc                                          >  1:180469/1‑139 (MQ=255)                                                                                                           cacCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCaaaaa                         <  8:173991/139‑1 (MQ=255)                                                                                                           cacCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCaaaaa                         <  5:392216/139‑1 (MQ=255)                                                                                                                        tgtgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCAGATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATc             >  8:3045/1‑138 (MQ=255)                                                                                                                        tgtgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATc             >  5:36971/1‑138 (MQ=255)                                                                                                                         gtgCGGACAAGAATCTCCGCAGTACCCATATTTGTAGCGAGCGGAATGGAATACACGTCACATAACCGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAGATGACAAGATCGa           <  2:262381/139‑1 (MQ=255)                                                                                                                          tgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTcgcg                               >  7:255159/1‑118 (MQ=255)                                                                                                                          tgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTcgcg                               <  8:255159/118‑1 (MQ=255)                                                                                                                          tgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAg          <  7:324862/139‑1 (MQ=255)                                                                                                                          tgCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAg          >  5:81019/1‑139 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                                                                    >  2:362476/1‑63 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                                                                    <  1:362476/63‑1 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTg        >  3:34990/1‑139 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATAGAGTg        >  4:296943/1‑139 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTAACCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCaaaaaaa                       <  6:311748/124‑1 (MQ=255)                                                                                                                            cgcACAAGAATTTCCGCAGTAACCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCaaaaaaa                       >  5:311748/1‑124 (MQ=255)                                                                                                                                  aGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTGCATTgg  <  6:81019/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTCAGCGGGTCGCGCAAAAAAATGACAAGATCGAGTGCATTGG  >  BS168_reference_genome/2358852‑2359119

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