Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,511,543	C→G	10.4%	A255G (GCT→GGT)	lutP →	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,511,543	0	C	G	10.4%	145.3 / 7.9	47	A255G (GCT→GGT)	lutP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base G (3/2);  total (25/22)
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.97e-01

ATCATCATCATGAGCGGGTTCAGGAAAGCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCTGCATACAGAG  >  BS168_reference_genome/3511410‑3511680                                                                                                                                      |                                                                                                                                          atcatcatcATGAGCGGGTTCAGGAAAGCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTg                                                                                                                                      >  8:208319/1‑139 (MQ=255)       atcatGAGCGGGTTCAGGAAAGCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATATCAGTACCAATTTTTTAGGGACTGAGCTGCCCGATGTTTTGTCAGATCTTGTTTGaa                                                                                                                                >  5:32497/1‑139 (MQ=255)            gAGCGGGTTCAGGAAAGCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCt                                                                                                                           >  6:307574/1‑139 (MQ=255)                            gCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGGTTGCTGTGtttt                                                                                                           >  5:339132/1‑139 (MQ=255)                            gCCTTGGAGATCTGGCCGGCGATTCTAGTATCCGGTGTTTCCTTTGCTGTTTTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATATTTTGTACGCGATTGTTTCAATGGCTGCGCTTGCGGTGtttt                                                                                                           >  8:3396/1‑139 (MQ=255)                             ccTTGGAGATCTGGCCGGCGATTCGCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATATCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCCGCTCTTGTTTCAATGGCTGCGCTTGCTGTGttttt                                                                                                          >  6:30023/1‑139 (MQ=255)                                tGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGa                                                                                                        >  1:11941/1‑138 (MQ=255)                                 ggAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGaa                                                                                                       <  1:93369/138‑1 (MQ=255)                                       cTGGCCGGCGATTCTCGTATCAGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTGAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAAtggtg                                                                                                 >  7:249611/1‑138 (MQ=255)                                               cGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGcc                                                                                           <  1:200679/136‑1 (MQ=255)                                               cGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGcc                                                                                           >  2:200679/1‑136 (MQ=255)                                               cGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGa                                                                                         >  7:186122/1‑138 (MQ=255)                                                              gtgtTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGtt                                                                         <  8:177543/139‑1 (MQ=255)                                                              gtgtTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGtt                                                                         <  4:269916/139‑1 (MQ=255)                                                               tgtTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGGGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGttt                                                                        >  2:49687/1‑139 (MQ=255)                                                               ggtTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAAGTTTTTAGGGCCTGGGGGGCCGGTGGTTGTTGCAGGTGTTGTTTCATGGGCTGCGCTGGCTTGGGTTGTGATTTGGTGGATGCCGAAACGACCATGCCGGttg                                                                        >  8:367407/2‑138 (MQ=255)                                                                 gttCCGTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGc                                                                      <  7:367407/138‑1 (MQ=255)                                                                    ccTTTGCTGTTGTTCACTATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCg                                                                    <  4:138427/138‑1 (MQ=255)                                                                    ccTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCg                                                                    <  5:30023/138‑1 (MQ=255)                                                                         gCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGaa                                                               <  1:22159/138‑1 (MQ=255)                                                                               gttCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGt                                                        >  2:46615/1‑139 (MQ=255)                                                                                  cAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCgg                                                     >  3:45499/1‑139 (MQ=255)                                                                                     tatCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGgccg                                                  >  4:97905/1‑139 (MQ=255)                                                                                       tctcAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGgccgc                                                 <  6:126979/138‑1 (MQ=255)                                                                                             ttCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGATATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGGCGGCTCaaat                                          >  5:33205/1‑137 (MQ=255)                                                                                               ccAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTg                                         >  6:332415/1‑138 (MQ=255)                                                                                                  aTTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGGTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAAc                                     >  6:31803/1‑139 (MQ=255)                                                                                                  aTTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGGTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAAc                                     >  6:26826/1‑139 (MQ=255)                                                                                                       ttAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTc                                                                                                                                  <  4:300176/41‑1 (MQ=255)                                                                                                       ttAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTc                                                                                                                                  >  3:300176/1‑41 (MQ=255)                                                                                                                      cccGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCcggcgg                 >  5:76349/1‑139 (MQ=255)                                                                                                                      cccGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCcggcgg                 <  5:332415/139‑1 (MQ=255)                                                                                                                      cccGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCcggcgg                 <  1:46615/139‑1 (MQ=255)                                                                                                                      cccGATGTTTTGTCAGCTCTTGTTTCAATGGATGCGCGTGCTGTGGTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGGGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCcggcgg                 <  5:70374/139‑1 (MQ=255)                                                                                                                        cGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGgcggcg               <  8:186122/139‑1 (MQ=255)                                                                                                                         gATGTTTTGTCAGGTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGTGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGgcgc              <  5:31803/139‑1 (MQ=255)                                                                                                                         gATGTTTTGTCAGGTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGgcgc              <  5:26826/139‑1 (MQ=255)                                                                                                                         gATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGgcgc              <  6:127455/139‑1 (MQ=255)                                                                                                                         gATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGgcgc              <  4:256572/139‑1 (MQ=255)                                                                                                                         gATGTTTTGGCAGCTCTTGTTTCAATGGCTGCGCTTGCTGGGGTTTTGAAATGGTGTAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATGCGGCGgcgc              <  6:351567/139‑1 (MQ=255)                                                                                                                           tGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGGGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCt            <  5:180963/139‑1 (MQ=255)                                                                                                                            gTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCt            <  3:97905/138‑1 (MQ=255)                                                                                                                                tGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCTGCATa       >  6:124059/1‑139 (MQ=255)                                                                                                                                tGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAAGTTAAACGGCGCGCCCAAAGCCGGCCGCGccgggcga       >  3:27810/1‑133 (MQ=255)                                                                                                                                tGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCACAATTGAAACGGCGCGCACCAAACCCGCGGCGCCGGCATa       >  6:191676/1‑139 (MQ=255)                                                                                                                                  tCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCTGCATACa     <  8:229002/139‑1 (MQ=255)                                                                                                                                     gCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCTGCATACagag  >  3:41162/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                          ATCATCATCATGAGCGGGTTCAGGAAAGCCTTGGAGATCTGGCCGGCGATTCTCGTATCCGGTGTTTCCTTTGCTGTTGTTCAATATCTCAGTTCCAATTTTTTAGGGCCTGAGCTGCCCGATGTTTTGTCAGCTCTTGTTTCAATGGCTGCGCTTGCTGTGTTTTTGAAATGGTGGAAGCCGAAAACAACATTCCGGTTTGCCGGTGAACAAGAGTCGGCCGCCTCAATTGAAACGGCGCGCACAAATCCGGCGGCGCCTGCATACAGAG  >  BS168_reference_genome/3511410‑3511680

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