Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	378,514	T→G	16.0%	V516G (GTT→GGT)	srfAA →	surfactin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	378,514	0	T	G	16.0%	111.9 / 6.8	58	V516G (GTT→GGT)	srfAA	surfactin synthetase
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base G (7/3);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCTGAACTGCTTTCTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGATTTCA  >  BS168_reference_genome/378376‑378651                                                                                                                                           |                                                                                                                                          aGGCTGAACTGCTTTCTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGt                                                                                                                                           >  3:297190/1‑139 (MQ=255) aGGCTGAACTGCTTTCTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGCGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGt                                                                                                                                           >  6:276469/1‑139 (MQ=255)   gCTGAACTGCTTTCTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGttt                                                                                                                                         <  4:28715/139‑1 (MQ=255)            cTTTCTGAGCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGGGCTGAAACAGACAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTttgattga                                                                                                                                  <  1:121839/137‑1 (MQ=255)                cTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGct                                                                                                                             <  7:284150/138‑1 (MQ=255)                cTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGct                                                                                                                             <  2:178715/138‑1 (MQ=255)                      atcCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGAATGAACGCTCTGaaac                                                                                                                      >  1:75699/1‑137 (MQ=255)                                  tATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAg                                                                                                           <  2:300398/138‑1 (MQ=255)                                  tATTTGAAGAGCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGGTTTGATTGAACGCTCTGAACGGATGATTACAg                                                                                                           <  8:192419/138‑1 (MQ=255)                                                      tCCTATCTAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCACGGACAGTCGTGTCGCGGGGTTGATTGAACGCTATGCACGGATGATTTCAGCTATCATTGGTATTTTAAAAg                                                                                      >  2:69193/1‑139 (MQ=255)                                                        cTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACg                                                                                                                      <  7:16672/105‑1 (MQ=255)                                                        cTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACg                                                                                                                      >  8:16672/1‑105 (MQ=255)                                                          aTCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTttgattga                                                                                                                                  <  5:309145/91‑1 (MQ=255)                                                          aTCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTttgattga                                                                                                                                  >  6:309145/1‑91 (MQ=255)                                                           tCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCcggc                                                                                 >  3:269763/1‑139 (MQ=255)                                                           tCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCcgg                                                                                  <  4:297190/138‑1 (MQ=255)                                                              aaCGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCcggcg                                                                                >  4:217821/1‑137 (MQ=255)                                                               aCGTTACATGAGCAATCTGCACGCATCGCCATTGTGCTGAAACAGAAAGGGGGTGGCCCGGACAGTCCTGTCGCGGGTTTGATTGAACGCTCTTAACGGATGATAACAGCGATCATGGGAATCTTAAAAGCcggcgg                                                                               <  1:232958/137‑1 (MQ=255)                                                                cGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGCTTGAACGCTCTGAACGGTTGATTACATCTTTCATTGGAATTTTAAAAGCCGGCGGAGc                                                                            >  8:24095/1‑139 (MQ=255)                                                                cGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGCACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGc                                                                            >  6:114954/1‑139 (MQ=255)                                                                cGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGc                                                                            >  6:244412/1‑139 (MQ=255)                                                                cGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGGTTTGGTTGAGCGCGGGGCACGGCTGATTTCCGGGATGATGGGTATTTTTAAAGCCGGGGGGGc                                                                            >  4:24612/1‑139 (MQ=255)                                                                  ttACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGGTTTGATTGAACGCGCTGAACGGATGATTACAGCGATCATGGGGATTTTAAAATCCGGCGGCGcc                                                                           >  2:235165/1‑138 (MQ=255)                                                                     cgtGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTAtg                                                                       >  4:311773/3‑139 (MQ=255)                                                                     cATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGGTTTGGTTTGAAGCGGTGAACCGCTGGTTAACGCGATCGATGGAATTTTAAAAGCCGGCGGAGCCAAtg                                                                       >  4:28204/1‑139 (MQ=255)                                                                        gAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGCTTTTTATTGAACGCACTGAACGGAAGAATACAGCTAACATGGGAATTTTAAAAGCCGGCGCAGCCTAtgatc                                                                    >  8:254801/1‑136 (MQ=255)                                                                         aGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGGTTGAAAGCTCGGGACGGGTGATTTCAGCGTTGAAGGGTTTTTTATAAGCCGGCGGCGCCGATGGGcc                                                                   >  2:256757/1‑139 (MQ=255)                                                                                tGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGGTCTGAACGGATGATTTCAGCGATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGAt                                                            >  2:192071/1‑139 (MQ=255)                                                                                  cACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGGCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGAt                                                            <  3:263889/137‑1 (MQ=255)                                                                                     gCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCggg                                                       <  4:269763/139‑1 (MQ=255)                                                                                     gCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCggg                                                       <  1:162712/139‑1 (MQ=255)                                                                                          gCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTc                                                  <  3:145441/139‑1 (MQ=255)                                                                                                 tgCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCctgct                                             >  2:253781/1‑137 (MQ=255)                                                                                                   ccgAAACCGAATGGTGTTGTCCCGGACAGTCCTGTCGCGGTTTTGCTTGAACGCTCTGAACGCATGCTTACAGCTATCATGGGTATTTTAAAAGCCGGCGGAGCCTATGTGACGCTTGATCCGGGGTTTCCTGATGAgc                                         <  6:50319/137‑1 (MQ=255)                                                                                                   cTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGGACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGCTTGATCCGGGTTTTCCTGCTGAgc                                         >  1:40425/1‑139 (MQ=255)                                                                                                   cTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAgc                                         <  1:235165/139‑1 (MQ=255)                                                                                                   cTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCAAGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAgc                                         >  7:258558/1‑139 (MQ=255)                                                                                                      aaaCAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCa                                      >  6:216293/1‑139 (MQ=255)                                                                                                             acggggTTGGCCCGGACAGTCCTGGCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAata                               <  3:24612/137‑1 (MQ=255)                                                                                                             aaGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAata                               <  3:311773/139‑1 (MQ=255)                                                                                                             aaGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAata                               <  5:114954/139‑1 (MQ=255)                                                                                                              aGGGGTTGGCCCGGACAGTCCTGTCGCGGGTTTGGTTGCACGCTCGGTACGGATGATTACAGCTATAATTGGAATTTTTAAAGCCGGCGGAGCCTATGTGGGGATTGATCCGGGTTTTGCTGGGGCGCGCCTTCAagat                              >  3:230200/1‑139 (MQ=255)                                                                                                                    gggCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGTTGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGATCCTATGTGCCGATTTATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTg                        <  4:230200/138‑1 (MQ=255)                                                                                                                     ggCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTg                        <  5:244412/138‑1 (MQ=255)                                                                                                                      gCCCGGACAGTCCTGTCGCGGTTTTGATTTAACGCTCTGAACGGATGATTAAAGCTATAATGGGAATTTTAAAAGCCCGCTGAGCCTATGTGCCGATTTATCCGGGTTTTACTGCTGAGCGCATTCAATATAatttt                        >  2:198162/1‑136 (MQ=255)                                                                                                                         cGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGTTTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTggagga                   <  3:261886/139‑1 (MQ=255)                                                                                                                           gACAGTCCTGTCGCGGTTTTGATTTAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAATCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACt                 >  5:272747/1‑139 (MQ=255)                                                                                                                                tCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTa                                                                                                              >  6:166674/1‑41 (MQ=255)                                                                                                                                tCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTa                                                                                                              <  5:166674/41‑1 (MQ=255)                                                                                                                                tCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCggg            >  2:195889/1‑139 (MQ=255)                                                                                                                                tCCTGTCGCGGGTTTTAATTAACGCGCTGAACCGATTATTACAGCTATAAAGGGAATTTTAAAAGCCCGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAAAAAATTTTTGAGGACTGCggg            >  4:261886/1‑139 (MQ=255)                                                                                                                                tCCTGTCGCGGGTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTAGGGGCCGAGTGATGCGGGTTTTCCTGCTGAGCGCATTCAATATTATTTTGAGGACTGCggg            >  3:169090/1‑139 (MQ=255)                                                                                                                                  cTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGc          <  1:6549/139‑1 (MQ=255)                                                                                                                                  cTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGc          <  8:258558/139‑1 (MQ=255)                                                                                                                                     tCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAata                               >  5:137535/1‑115 (MQ=255)                                                                                                                                     tCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAata                               <  6:137535/115‑1 (MQ=255)                                                                                                                                      cgcgGGTTTGATTGAACGCTCTGTACGGATGATTACAGCTATCTTGGGAATTTTAAAAGCCGGCGGAGCGTATGTGCCGTTTGATCGGGGTTTTCGTTCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGAt      <  4:21795/139‑1 (MQ=255)                                                                                                                                        cgGTTTTGATTGAACGCTCTGGACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGAttt    <  2:40425/139‑1 (MQ=255)                                                                                                                                          gTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGATTTCa  >  8:199101/1‑139 (MQ=255)                                                                                                                                          gTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGATTTCa  >  1:24857/1‑139 (MQ=255)                                                                                                                                          gTTTTGATTGAACGCTCTGAACCGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGATTTCa  >  3:30398/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                          AGGCTGAACTGCTTTCTGATCATCCAGCGGTTGTATTTGAAGATCGCACATTGTCCTATCGAACGTTACATGAGCAATCTGCACGCATCGCCAATGTGCTGAAACAGAAAGGGGTTGGCCCGGACAGTCCTGTCGCGGTTTTGATTGAACGCTCTGAACGGATGATTACAGCTATCATGGGAATTTTAAAAGCCGGCGGAGCCTATGTGCCGATTGATCCGGGTTTTCCTGCTGAGCGCATTCAATATATTTTGGAGGACTGCGGGGCGGATTTCA  >  BS168_reference_genome/378376‑378651

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