Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,755,184	Δ1 bp	6.7%	coding (108/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,184	0	C	.	6.7%	211.3 / 7.5	60	coding (108/1359 nt)	ywoD	putative efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (17/39);  new base . (2/2);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AATAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGAAA  >  BS168_reference_genome/3755046‑3755311                                                                                                                                           |                                                                                                                                aaTAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGc                                                                                                                                 <  1:233069/139‑1 (MQ=255)  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)                                                                                                                                           |                                                                                                                                AATAAGCCCGAATAAAAACAAGCGTTTCCGCCCAATCAAATCACCGAGCTTTCCCGTTACGGGCTGCGCTACGGCACTCGCCAAGTAAAACGAGGAAATCAGCCAGGAAACCGTTGTAAACGAGAGGTGAAACTCATGCTGGATGCTATGGAGCGCAAGAGAGATCATAGACGAGTTCAGCGGATTCAGCATTGTTCCCGAAGCCACAGCTGTTAAAAAAACAGCTCGGCTTTTTTTCGATAACATCATTTCACTTGTTCCAGAAA  >  BS168_reference_genome/3755046‑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: ‑