Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	389,065	T→G	14.5%	L441R (CTT→CGT)	srfAB →	surfactin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	389,065	0	T	G	14.5%	77.3 / 3.5	45	L441R (CTT→CGT)	srfAB	surfactin synthetase
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (5/2);  total (16/29)
Fisher's exact test for biased strand distribution p-value = 3.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGGACCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCACCGCCGTCGTATATGA  >  BS168_reference_genome/388933‑389203                                                                                                                                     |                                                                                                                                           cGTGGACCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTggg                                                                                                                                      <  4:5265/139‑1 (MQ=255) cGTGGACCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTggg                                                                                                                                      <  3:405531/139‑1 (MQ=255) cGTGGACCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTggg                                                                                                                                      <  2:297427/139‑1 (MQ=255)      aCCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACg                                                                                                                                 <  6:100570/139‑1 (MQ=255)              aaTTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGagcag                                                                                                                         <  2:18619/139‑1 (MQ=255)                        aaCGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTTTGCGTGGGTTCGGGCGTGTGGGGGCCGGAGCGCAGCGGCCAATTg                                                                                                               >  1:25998/1‑139 (MQ=255)                          cGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTTTGGGTGTGGAAGGGCGTGGGGGAGACGAACAGCAGCGGCGAATTTTc                                                                                                             >  3:401581/1‑139 (MQ=255)                           ggcgccgccTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                                            <  4:349871/136‑1 (MQ=255)                           ggAGCCGCCTTTGAATCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGAGCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCg                                                                                                            <  5:134947/139‑1 (MQ=255)                             cgccgccTTTGATTCCGCTTTTATCGAACGGCCAGCGGAGCACCTTACCCGCATGATGGAAGCAGCGGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAg                                                                                                          <  4:228960/138‑1 (MQ=255)                             aGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAg                                                                                                          <  1:80406/139‑1 (MQ=255)                                            cGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTGCGTGAGTACGGGCTTGTGGGTGACGCACAGCAGCGGCAAATTGTAGAGGTGTTTAAAAGCCCg                                                                                           >  7:166115/1‑139 (MQ=255)                                             gCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGa                                                                                          <  1:280742/139‑1 (MQ=255)                                                  tATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTGTGCGTGTGTACGGGGGTGTGGGGGGCGCAGAGCAGCGGCAAATTGTCGAGGGGTTTTACCGCCCGCAAGcc                                                                                     >  6:305915/1‑139 (MQ=255)                                                     acagcGGTCAGCGGAGCACCTTACCCGCATGATGGAAGCTGCTGTCGATCAGCCGGCCGATTTTGTGCTTGAGTACGGGATTGTGGGAGACGTAGAGCAGCGGGAAATTGTCGAGGTGTTTAACAGCACGAAATCCGaa                                                                                  <  1:138672/135‑1 (MQ=255)                                                      gcaCGGTCTGCGGAGCACCTTTCCCGCATGATGGAAGCAGCTGTCGAGCAGCGGGCGGCTTTTGTGCGTTAGTCCGGGATTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGGGTTTAACAGCACGAAAGCCGaa                                                                                  <  2:117045/136‑1 (MQ=255)                                                      gAACGGACCGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCGTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGaa                                                                                  <  6:138464/138‑1 (MQ=255)                                                      acagggACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAAc                                                                                 <  7:186813/135‑1 (MQ=255)                                                           gACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCTTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCt                                                                            <  7:344430/139‑1 (MQ=255)                                                           gACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCt                                                                            <  1:7741/139‑1 (MQ=255)                                                             cAGCGGAGCTCCTTACCCGCATGTTGGTAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGGGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAATTTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGa                                                                          <  5:427505/139‑1 (MQ=255)                                                             cAGCGGAGCACCTTACCCGCATGATCGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGa                                                                          <  6:251342/139‑1 (MQ=255)                                                                     cACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCCGCAAAATTTAGAGGTGTTTTACAGCCCCAAAGCCGCACACCCAGAAggcccgg                                                                  >  8:186813/1‑134 (MQ=255)                                                                      aacTTACCGGCATGATGGAAGCAGCTGTCGATCAGCGGGCCGAGTTTGTGCGTGAGTACTGGCTTGTGTGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGc                                                                 <  6:393753/137‑1 (MQ=255)                                                                      aCCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGc                                                                 <  4:76798/139‑1 (MQ=255)                                                                      aCCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGc                                                                 <  5:309742/139‑1 (MQ=255)                                                                        cTTACCCGCATGATGGCAGCAGCTGGCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTg                                                               <  1:334273/139‑1 (MQ=255)                                                                          tACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGGGTGATTACGGGCTTGTGGGTGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGtt                                                             >  4:270607/1‑139 (MQ=255)                                                                          tACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGGGTGTGGGAGACGAAGAGCAGCGGGAAATTGGGGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTggt                                                             >  1:50520/1‑137 (MQ=255)                                                                          tACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTTTGCGGGGGTACGGGGTTTTGGGAGACGCAGAGGCGCGGCAAAATTTCGAGGTGTTTATAACAAAGAAAGCCACACTCACTGCATGAATGGCTGt                                                              >  5:393753/1‑138 (MQ=255)                                                                             ccGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTGTTGTGGGGGAGGTCGGGGATGTGGGTGACGCAGAGCAGCGGCAAATTGTCGATGGGTTTAACAGCACGAAAGCCGCACGCCCTGAAGGAATGGCTGTTCAc                                                          >  7:328127/1‑139 (MQ=255)                                                                               gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTTTTGGTGAGTACCGGCTTTTTGGGGACCCAGAGCAGCGGCAAAATTTCGAGGTGGTTAAAAGCACCAAATCCGCACTCCCAGAAGGAATAGCTGGTGAcc                                                         >  6:309742/1‑138 (MQ=255)                                                                               gCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCGTTTTTGCGAGAGTACCGGGTTGTGTGAGAGAAATAGCACCGGGAAAATGTCGCGGTGTTTAAAAACACCAAACCCCCACACCCTGAAGGACAGGGTGTTCAcc                                                         >  2:334273/1‑138 (MQ=255)                                                                               gCATGATGGAAGCAGCTGTCGATAAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGCCGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCa                                                        <  8:328127/139‑1 (MQ=255)                                                                               gCATGATGGAAGCAGCTGTCGAGCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCa                                                        <  3:295639/139‑1 (MQ=255)                                                                                             gcTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGaa                                                                                  >  7:279230/1‑99 (MQ=255)                                                                                             gcTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGaa                                                                                  <  8:279230/99‑1 (MQ=255)                                                                                              cTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCGTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGGAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCa                                         <  8:166115/139‑1 (MQ=255)                                                                                                           ggccgCGTTTTTGGGTGGGTAAGGGCGTGGGGGGGGCGAAGAGCAGCGGCAAATTTTCGGGGTGGTTTACAGCACGTAAGCCGCACTACCTGAAGGAATGGGTGTTGACCAACTGTTTGAAGAGCCACCGAAACGGAcc                            >  2:111003/1‑138 (MQ=255)                                                                                                                       gCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAcc                >  5:365636/1‑139 (MQ=255)                                                                                                                        cGTGAGTACGGGCTTGTGGGAGACGATGAGCAGCGGAAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCGGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccg               <  4:401581/139‑1 (MQ=255)                                                                                                                         gTGAGTACGGGCTTGGGGGAGACGAAGAGCAGCGGCAAATTGGCGAGGGGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccgc              <  6:16534/139‑1 (MQ=255)                                                                                                                          tGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccgcc             >  7:272023/1‑139 (MQ=255)                                                                                                                          tGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccgcc             <  3:182722/139‑1 (MQ=255)                                                                                                                           gAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCAccgccg            >  7:220435/1‑139 (MQ=255)                                                                                                                                     ttGTGGGTGACGACGTGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCTAACTCCCTTAAGTAATGGCTGTTCACCATGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCACCGCCGTCGTATATGa  <  4:71608/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                           CGTGGACCGTCAAAATTAAATATAACGGAGCCGCCTTTGATTCCGCTTTTATCGAACGGACAGCGGAGCACCTTACCCGCATGATGGAAGCAGCTGTCGATCAGCCGGCCGCTTTTGTGCGTGAGTACGGGCTTGTGGGAGACGAAGAGCAGCGGCAAATTGTCGAGGTGTTTAACAGCACGAAAGCCGAACTCCCTGAAGGAATGGCTGTTCACCAAGTGTTTGAAGAGCAAGCGAAACGGACGCCGGCGAGCACCGCCGTCGTATATGA  >  BS168_reference_genome/388933‑389203

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