Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	272,231	T→G	25.9%	A144A (GCT→GCG)	gudD →	glucarate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	272,231	0	T	G	25.9%	79.1 / 4.9	51	A144A (GCT→GCG)	gudD	glucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (13/24);  new base G (11/3);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTTGGCGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGTCTGCGCC  >  BS168_reference_genome/272097‑272364                                                                                                                                       |                                                                                                                                      cAGTTTGGCGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGc                                                                                                                                   <  8:351016/139‑1 (MQ=255) cAGTTTGGCGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGGCGCGGCTTTGc                                                                                                                                   <  4:123284/139‑1 (MQ=255)   gTTTGGCGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTg                                                                                                                                 >  5:278693/1‑139 (MQ=255)         cGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGaa                                                                                                                           <  7:269730/139‑1 (MQ=255)         cGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGaa                                                                                                                           <  5:86641/139‑1 (MQ=255)                     gggAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGGTGGGCGCAGGGCAGCAGCGt                                                                                                               >  7:20590/1‑139 (MQ=255)                       gAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGtga                                                                                                             <  7:33000/139‑1 (MQ=255)                       gAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGtga                                                                                                             <  5:342210/139‑1 (MQ=255)                          ggCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGCACCGGTCGCGGCGTTGCTGGGGGGAGGGGGGCAGCgggttg                                                                                                           >  7:67329/1‑133 (MQ=255)                          ggCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGtgatg                                                                                                           >  6:88240/1‑138 (MQ=255)                          ggCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTTCGTGGCGAACGGACTCAGCGggatt                                                                                                           >  6:342210/1‑137 (MQ=255)                          ggCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGCTGGGCGAAGGGAAGCAGCGtgatg                                                                                                           >  8:67609/1‑138 (MQ=255)                               ggAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGGGTTGGTGGGGGAAGGGATGCAGCGGGTTgtagta                                                                                                      >  4:205812/1‑138 (MQ=255)                                    ccAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTaaaa                                                                                                   >  8:162744/1‑136 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  7:307217/138‑1 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  7:162744/138‑1 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  3:249356/138‑1 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  5:88240/138‑1 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  6:253267/138‑1 (MQ=255)                                              gATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCGGGGCGAGGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAc                                                                                       <  7:121265/138‑1 (MQ=255)                                                    cgcACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGGGTTGGTGGGCGCAGGGTAGCAGCGTGAGGAAGTAAAAATGCTAGGGTACCGCgttt                                                                                >  6:204961/1‑139 (MQ=255)                                                         cacTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGGGTTGGTTGGGGGAGGGGAGCAGCGGTGTGAAGTAAAAATGGTGGGGTACCTCTTTTTCATc                                                                           >  6:280672/1‑139 (MQ=255)                                                              tgCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGGTGGGGGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTTCCTCTGTTTCATCGGCGa                                                                      >  8:22714/1‑139 (MQ=255)                                                              tgCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTTGGGGGGGCAGGGTAGCAGCGTGATGAAGTAAAAATTCTCGGGTTCCGCGTTTTCATAGGCgg                                                                      >  3:303717/1‑138 (MQ=255)                                                              tgCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGCTGGGCGAAGGGAAGCAGCGTGATGAATTAAAAATGCTCGGTTACCTCTTTTAAATCGGCgc                                                                      >  7:123736/1‑138 (MQ=255)                                                              tgCACGCCGTCACCGCACTGGAAGCGGCTTTGGTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGGTGGGGGGAGGGAAGCCGCGTGATGAAGGAAAAATGCTCGGTTTCCTCTTTTTCCTCGGCGa                                                                      >  5:160215/1‑139 (MQ=255)                                                              tgCACGCCGTCACCGCACTGGAAGCGGCTTTGATGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTTCAGGGGGAATCGAAACAACGGGAAGAATTTAAAAATCTAGGTTAACAAATTTACATACGCGa                                                                      >  8:146797/1‑139 (MQ=255)                                                                   gCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCgg                                                                                           >  7:105763/1‑113 (MQ=255)                                                                   gCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCgg                                                                                           <  8:105763/113‑1 (MQ=255)                                                                   gCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGa                                                                  >  8:179968/1‑138 (MQ=255)                                                                   gCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGGTGGGCGCAGGGAAGCAGCGGGATGAATTAAAAATGCTCGGGTTCCGCTTTTTCATAGGCGCCCgc                                                                  >  2:116969/1‑137 (MQ=255)                                                                        cACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGGAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTTCATCGGCGACCGaaagag                                                             >  3:98330/1‑135 (MQ=255)                                                                         aCCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGga                                                           <  3:6418/139‑1 (MQ=255)                                                                         aCCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGGCGCGGCGTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGga                                                           <  3:212914/139‑1 (MQ=255)                                                                             cACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTTCTGGGCGAAGGGAACCAGAGTGATGAATGAAATATGGACGGTTAACTCGTTTACATCGGCGAACGCAATAGgacggc                                                       >  7:68792/1‑137 (MQ=255)                                                                             cACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGCTGGGCGCAGGGAAGCAGCGTGATGAAGTAAAAATTCTCGGTTACCTCTTTTTCATCGGCGCCCGAAATAGgacgcc                                                       >  2:220205/1‑137 (MQ=255)                                                                                  gAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGc                                                  <  6:78048/139‑1 (MQ=255)                                                                                               tgGATTTGCTTGGACAGTTTCTTCAGGACCCGGTCGCGGCGTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGa                                     <  7:140877/139‑1 (MQ=255)                                                                                                 gATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATg                                   <  6:278693/139‑1 (MQ=255)                                                                                                 gATTTGCTTGGAAAGTTTCTTCAGGAACCGGACGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCTAAATCGGACGACCTTGCCGTATCAAAGCGATg                                   <  2:115767/139‑1 (MQ=255)                                                                                                   tttGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCCGCGTGATGAAGTAAAAATGCTAGGTTACCTCGTTTACATCGGCGACCGAAATCGGACGCCCTTGCCGTATCAAAGCGATGAg                                 >  6:279414/1‑139 (MQ=255)                                                                                                   tttGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTAACTCTTTTACAAAGGCGACCGAAATCGGACGCCCTTGCCGTATCAAAGCGATGAg                                 >  3:106660/1‑139 (MQ=255)                                                                                                    ttGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCGTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAATTCGGACGACCTTGCCGTATCAAAGCGATGAg                                 <  7:277563/138‑1 (MQ=255)                                                                                                       cTTGTAAAGTTTCTTCAGGAACCGGTCGAGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGt                             <  3:205812/139‑1 (MQ=255)                                                                                                       cTTGGTTAGTTTATTCAGGAACCGGTCGCGGCTTTGCGGGGCGATGGGAAGCAGGGTGATGAAGTAAAAATGCTCGGTTTCCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGt                             <  7:146797/139‑1 (MQ=255)                                                                                                           gAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGGAGGGCAGGAGCGGGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGCCCGAAATAGGACGACCATGCCGTATCAAAGCGCTGCGCAGt                             >  8:188452/1‑135 (MQ=255)                                                                                                                           aCCGGGCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGt         <  5:197726/139‑1 (MQ=255)                                                                                                                            ccGGTCGCGGCTTGGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGTc        <  8:67329/139‑1 (MQ=255)                                                                                                                              ggTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGTCTg      <  7:179968/139‑1 (MQ=255)                                                                                                                                tCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGAGGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGTCTgcg    <  3:63479/139‑1 (MQ=255)                                                                                                                                  gcgGGTTGGCGGGGCGAAGTGAAGCAGCGGGATTAAGTAAAAATGCGGGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGGCAGATTGTGCGTGGTTCCGGCTGCGcc  <  5:48650/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      CAGTTTGGCGATCAGGACAGGGGAGGGCGAGGAAACCAGACGTTTGATTTGCGCACCACTGTGCACGCCGTCACCGCACTGGAAGCGGCTTTGTTGGATTTGCTTGGAAAGTTTCTTCAGGAACCGGTCGCGGCTTTGCTGGGCGAAGGGAAGCAGCGTGATGAAGTAAAAATGCTCGGTTACCTCTTTTACATCGGCGACCGAAATCGGACGACCTTGCCGTATCAAAGCGATGAGCAGTCAGATTGTGCGTGGTTCCGTCTGCGCC  >  BS168_reference_genome/272097‑272364

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