Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,755,185	Δ1 bp	7.5%	coding (107/1359 nt)	ywoD ←	putative efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,755,185	0	T	.	7.5%	150.4 / 7.1	53	coding (107/1359 nt)	ywoD	putative efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (15/31);  new base . (2/2);  total (17/36)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

ATAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGAAA  >  BS168_reference_genome/3755047‑3755311                                                                                                                                           |                                                                                                                               aaaaccccGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                <  2:67360/134‑1 (MQ=255) aTAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                >  1:132099/1‑139 (MQ=255) aTAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                <  5:314931/139‑1 (MQ=255) aTAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                <  3:172874/139‑1 (MQ=255) aTAAGCCCGAATAAAAACAAGCGGTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                <  4:37360/139‑1 (MQ=255)          aaTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTAt                                                                                                                       >  5:365780/1‑139 (MQ=255)          aaTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTAt                                                                                                                       >  5:120710/1‑139 (MQ=255)          aaTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATG‑‑GGATGCTATgg                                                                                                                     >  3:287789/1‑139 (MQ=255)          aaTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATG‑‑GGATGCTATgg                                                                                                                     >  3:282699/1‑139 (MQ=255)           aTAAACACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCGGCGCTACGGCACTCGACAAGTAAACCGAGGAAATCAGCCAGGAACCCGTTGTAAACGAGAGGTGAAACTCATGCAGGATGCTATg                                                                                                                      <  6:200205/139‑1 (MQ=255)           aTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAATCTCATGCTGGATGCTATg                                                                                                                      <  7:280217/139‑1 (MQ=255)           aTAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATg                                                                                                                      <  3:225221/139‑1 (MQ=255)                caCAAGCGTTTCCGCCCAATCACATCACCGAGCTTTCCCGGTACGGGCTGCGCTACGGCACTGGCCAAGTAAAACGAGGAACTCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAgcg                                                                                                                 <  1:104413/138‑1 (MQ=255)                aaCAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAgcg                                                                                                                 <  4:36965/139‑1 (MQ=255)                  cAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCa                                                                                                               <  8:99092/139‑1 (MQ=255)                     gCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACTGGCTGCGCTACGGCACTCTCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAaga                                                                                                            <  7:220565/139‑1 (MQ=255)                     gCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAaga                                                                                                            <  2:132099/139‑1 (MQ=255)                     gCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAaga                                                                                                            <  8:115936/139‑1 (MQ=255)                       gtatCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  5:118887/136‑1 (MQ=255)                       gTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  1:320714/139‑1 (MQ=255)                       gTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  3:332348/139‑1 (MQ=255)                       gTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCACAagaga                                                                                                          <  6:80625/139‑1 (MQ=255)                       gTTTACGACCAATCAAATCACCGAGCTTTCCCGTTACGGTCTTCGCTAAGGCACTCGCCTAGTAAAACTAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  7:169410/139‑1 (MQ=255)                       gTTGCCGCCCAATCAAATCACCGATCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  3:36837/139‑1 (MQ=255)                          tCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAaga                                                                                                            >  3:63417/1‑134 (MQ=255)                               ccAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATG‑‑GGATGCTATGGAGCGCAAGAGAGATCATAGAc                                                                                                <  4:282699/139‑1 (MQ=255)                               ccAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATG‑‑GGATGCTATGGAGCGCAAGAGAGATCATAGAc                                                                                                <  4:287789/139‑1 (MQ=255)                                 aatcaactcaCCGCGCTGTCCCGTGACGGGCGGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGAc                                                                                                <  6:45283/139‑1 (MQ=255)                                 aatcaaatcaCCGAGCTTTCCCGTTACGGGCTGCGCTACCGCCCTCGCCAAGTAAAACGACGAAATCAGCCAGGAAACCTTTGTAAACGAGATGTGAAACTCATGCTGGATGCTATGGATCGCAAGAGAGATCATAGAc                                                                                                <  8:10928/139‑1 (MQ=255)                                     aaatcaCCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGt                                                                                            <  8:114910/139‑1 (MQ=255)                                       atcaCCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTc                                                                                          <  5:295987/139‑1 (MQ=255)                                        tcaCCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCa                                                                                         <  6:365780/139‑1 (MQ=255)                                          aCCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGc                                                                                       <  1:133581/139‑1 (MQ=255)                                          aCCGAGCTTTCCCGTTACGGGCGGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTTAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGc                                                                                       <  2:223853/139‑1 (MQ=255)                                                   tCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAg                                                                               >  1:335024/1‑138 (MQ=255)                                                       gtgtcGGGGTGCGCTAAGGCACTCGCCAATTAATTCGAGGAAATGAGCCATGATTCAGTTGTAAAGGAGATGTGTAAGTCATTCAGGATGCTATGGTGCGCAAGAGAGATCATAGACGAGTTCAGCTGTTTCAGCATTg                                                                          <  8:194894/135‑1 (MQ=255)                                                       gTTACGGGCTTCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTg                                                                          <  1:99538/139‑1 (MQ=255)                                                        ttACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGt                                                                         >  3:167228/1‑139 (MQ=255)                                                            gggCTGCGCTACGGCACTCGACAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTccc                                                                     <  4:63417/139‑1 (MQ=255)                                                               cTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATg                                                                                                                      >  4:312451/1‑87 (MQ=255)                                                               cTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATg                                                                                                                      <  3:312451/87‑1 (MQ=255)                                                                  cgcTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTAGGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGAATTGTTCCCGAAGcc                                                               >  6:81218/1‑139 (MQ=255)                                                                      aCGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAg                                                           >  1:15047/1‑139 (MQ=255)                                                                            cTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGtt                                                      >  5:181965/1‑138 (MQ=255)                                                                               gCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                >  4:144361/1‑61 (MQ=255)                                                                               gCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCt                                                                                                                                <  3:144361/61‑1 (MQ=255)                                                                                 cAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTaaaaa                                                 >  2:324479/1‑138 (MQ=255)                                                                                    gTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTaaaaaaa                                               >  1:100182/1‑137 (MQ=255)                                                                                             ggAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGa                                                                   <  6:160687/108‑1 (MQ=255)                                                                                             ggAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGa                                                                   >  5:160687/1‑108 (MQ=255)                                                                                              caaaTCAGCCAGGAAACCGTTGTACACGAGAGGTGAAACTCATGCGGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCttt                                   <  6:214211/138‑1 (MQ=255)                                                                                                      ccAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGAt                           <  6:120710/139‑1 (MQ=255)                                                                                                      ccAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGAt                           <  2:335024/139‑1 (MQ=255)                                                                                                        aGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATaa                         <  5:149845/139‑1 (MQ=255)                                                                                                                         gagaGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          >  1:247705/1‑41 (MQ=255)                                                                                                                         gagaGGTGAAACTCATGCTGGATGCTATGGAGCGCAagaga                                                                                                          <  2:247705/41‑1 (MQ=255)                                                                                                                            agGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTc        >  1:216025/1‑136 (MQ=255)                                                                                                                              gTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGa    <  2:100182/138‑1 (MQ=255)                                                                                                                              gTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTGTTTTTCGATAACATCATTTCACTTGTTCCAGa    >  6:200486/1‑138 (MQ=255)                                                                                                                               ggAAACTCATGCTGGCTGCTATGGAGCGCAAGAGAGATCATCGACGAGGTCAGCGGATTCAGCATTGGTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGaca  <  8:308688/138‑3 (MQ=255)                                                                                                                                           |                                                                                                                               ATAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGAAA  >  BS168_reference_genome/3755047‑3755311

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