Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	113,389	A→G	7.7%	E197G (GAA→GGA)	cysE →	serine O‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	113,389	0	A	G	7.7%	143.0 / 2.8	51	E197G (GAA→GGA)	cysE	serine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCATGATGTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGAAAAGTGAAA  >  BS168_reference_genome/113253‑113524                                                                                                                                         |                                                                                                                                        ctgcATGATGTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGaa                                                                                                                                        >  3:285298/1‑138 (MQ=255) ctgcATGATGTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGAGTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAAc                                                                                                                                       >  4:48739/1‑139 (MQ=255) ctgcAAGATGTTCCTGATTTTTCAACAGTTGTCGTTATCCCTGGACGGGTCGTTGTTCAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGTCCGATCCAGTATCTGACCCCGTTAAGTCATTAGaa                                                                                                                                        >  2:291572/1‑138 (MQ=255)     atgatgTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAAcagca                                                                                                                                   <  4:285298/139‑1 (MQ=255)      tgatgTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAAcag                                                                                                                                     >  3:155135/1‑136 (MQ=255)         tgTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGAtt                                                                                                                               >  5:39859/1‑139 (MQ=255)          gTTCCTGATTTTTCAACAGTGGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGAttt                                                                                                                              <  6:194442/139‑1 (MQ=255)                gATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTGGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGAGTTGCCCGAGCCAGTTGCTGACCGCGTTAAGTCTTTGGAACATCAGATTTTAGAGc                                                                                                                        >  6:57292/1‑139 (MQ=255)                                tCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGGTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGa                                                                                                        <  2:53574/139‑1 (MQ=255)                                            gACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCt                                                                                                                       >  4:9505/1‑112 (MQ=255)                                            gACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCt                                                                                                                       <  3:9505/112‑1 (MQ=255)                                                 gTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGa                                                                                       >  4:137331/1‑139 (MQ=255)                                                 gTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGa                                                                                       >  2:76820/1‑139 (MQ=255)                                                 gTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGa                                                                                       >  8:185413/1‑139 (MQ=255)                                                 gTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGa                                                                                       >  6:162122/1‑139 (MQ=255)                                                 gTCGGTGTACAAAATGGAAAGAAAGTAAGACGCGATTTAACCCATAAGGATTTGCCCGATCCAGTTGCTGAACGCTTTACGTCTTTGGGACCTCAGATTGTAGAGCTGAATGAAGAACGTGAAAAAAGAaacgaaagtc                                                                                       >  4:150460/1‑137 (MQ=255)                                                  tCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGat                                                                                      >  6:156327/1‑139 (MQ=255)                                                     ttGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGatcaa                                                                                   <  2:35046/139‑1 (MQ=255)                                                           aaaaTGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTTCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCaaa                                                                              >  3:28261/1‑138 (MQ=255)                                                             aaTGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAAt                                                                           >  1:28994/1‑139 (MQ=255)                                                                  cAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTa                                                                                                                            >  1:127788/1‑85 (MQ=255)                                                                  cAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTa                                                                                                                            <  2:127788/85‑1 (MQ=255)                                                                     gAAAGTAAGACTCGATTTAAAACATCATGATTTGCCCGCTCCAGGTGCTGCCCGCTTTAAGTCTTTGGGACCGCAGATTTTAGAGCTGAATGCAGAACTTGAAGACAGAACAGAAAGGATCAATCAAAACTGACAATca                                                                   <  3:150460/139‑1 (MQ=255)                                                                     gAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATca                                                                   <  2:28994/139‑1 (MQ=255)                                                                     gAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATca                                                                   <  7:135900/139‑1 (MQ=255)                                                                       aattAATACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATcaca                                                                 <  7:6359/136‑1 (MQ=255)                                                                       aaGTAAGACGCGATTTAAACCATCAGGTTTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATcaca                                                                 <  1:113446/139‑1 (MQ=255)                                                                       aaGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATcaca                                                                 <  7:185413/139‑1 (MQ=255)                                                                       aaGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATcaca                                                                 <  7:287604/139‑1 (MQ=255)                                                                       aaGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATcaca                                                                 <  5:156327/139‑1 (MQ=255)                                                                         gTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGTCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACt                                                               <  8:98297/139‑1 (MQ=255)                                                                              aCGCGATTTAAACGATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTtata                                                          <  6:39859/139‑1 (MQ=255)                                                                                  gATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTataata                                                       >  4:205538/1‑138 (MQ=255)                                                                                  gATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATAc                                                      <  5:162122/139‑1 (MQ=255)                                                                                  gATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATAc                                                      <  3:205538/139‑1 (MQ=255)                                                                                            aTCAGGATTTGTCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGa                                            <  1:291572/139‑1 (MQ=255)                                                                                              cAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGa                                            >  6:223994/1‑137 (MQ=255)                                                                                              cAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACa                                          >  8:155136/1‑139 (MQ=255)                                                                                                           aTCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTc                             >  1:69379/1‑139 (MQ=255)                                                                                                           aTCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTc                             >  5:101667/1‑139 (MQ=255)                                                                                                                 ttGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGtt                          <  2:262510/136‑1 (MQ=255)                                                                                                                 ttGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGtt                          >  1:262510/1‑136 (MQ=255)                                                                                                                      gACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATca                                                                   >  2:217261/1‑90 (MQ=255)                                                                                                                      gACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATca                                                                   <  1:217261/90‑1 (MQ=255)                                                                                                                      gACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTg                   >  6:254030/1‑138 (MQ=255)                                                                                                                      gACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTg                   >  8:106896/1‑138 (MQ=255)                                                                                                                        ccGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAag                <  6:101667/139‑1 (MQ=255)                                                                                                                               aaGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGaa          <  1:76820/138‑1 (MQ=255)                                                                                                                                      tGGGACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGa                                            >  7:212825/1‑97 (MQ=255)                                                                                                                                      tGGGACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGa                                            <  8:212825/97‑1 (MQ=255)                                                                                                                                       ggAACAGCAGATTTTAGAGCTGAAGTCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGAAAAGTGaaa  <  5:57292/138‑1 (MQ=255)                                                                                                                                       ggAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGAAAAGTGaaa  <  2:69379/138‑1 (MQ=255)                                                                                                                                       ggAACAGCAGATTTTAGAGCAGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGAAAAGTGaaa  <  4:155135/138‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CTGCATGATGTTCCTGATTTTTCAACAGTTGTCGGTATCCCTGGACGGGTCGTTGTACAAAATGGCAAGAAAGTAAGACGCGATTTAAACCATCAGGATTTGCCCGATCCAGTTGCTGACCGCTTTAAGTCTTTGGAACAGCAGATTTTAGAGCTGAAGGCAGAACTTGAAGACAGAAAAGAAAGGATCAATCAAAAATGACAATCACACTTTATAATACATTGACTAGACAGAAGGAAACATTCGTTCCTCTTGAAGAGGGAAAAGTGAAA  >  BS168_reference_genome/113253‑113524

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