Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,041,654	G→A	intergenic (+29/‑125)	dsbA → / → b3861	protein disulfide isomerase I/putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,041,654	0	G	A	100.0%	342.9 / NA	87	intergenic (+29/‑125)	dsbA/b3861	protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (3/84);  total (3/84)

ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCA  >  NC_000913/4041564‑4041730                                                                                           |                                                                              atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  2:895997/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  1:5303682/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  2:859692/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  1:5466874/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  1:5850148/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  2:5014666/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  2:4697888/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                             <  1:882915/93‑1 (MQ=255)   atGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATAtt                                                                           <  2:702919/93‑1 (MQ=255)    tGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTc                                                                          <  2:3643507/93‑1 (MQ=255)    tGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTc                                                                          <  1:6725630/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                         <  1:833742/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                         <  2:5925322/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                         <  2:7480201/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                         <  1:1612112/93‑1 (MQ=255)      gATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAc                                                                        <  2:4416400/93‑1 (MQ=255)      gATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAc                                                                        <  1:6736110/93‑1 (MQ=255)      gATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAc                                                                        <  1:2949006/93‑1 (MQ=255)      gATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAc                                                                        <  1:5092826/93‑1 (MQ=255)       aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACg                                                                       <  1:7195102/93‑1 (MQ=255)       aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACg                                                                       <  1:3836689/93‑1 (MQ=255)       aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACg                                                                       <  1:4605431/93‑1 (MQ=255)        tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACgg                                                                      <  1:1180044/93‑1 (MQ=255)        tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACgg                                                                      <  1:2761485/93‑1 (MQ=255)        tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACgg                                                                      <  1:3797471/93‑1 (MQ=255)         gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACggg                                                                     <  2:1294874/93‑1 (MQ=255)         gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACggg                                                                     <  2:7160258/93‑1 (MQ=255)         gTTTTTGTTCAGCAGTAAGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACggg                                                                     <  2:6132992/93‑1 (MQ=255)          tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                    <  2:3909230/93‑1 (MQ=255)           ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGC‑t                                                                  <  2:4723719/93‑1 (MQ=255)           ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                                                  <  2:6651758/93‑1 (MQ=255)            tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                                                 <  1:4764421/93‑1 (MQ=255)             ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                <  2:3988259/93‑1 (MQ=255)             ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                <  2:1026264/93‑1 (MQ=255)              tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                               <  1:2074507/93‑1 (MQ=255)              tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                               <  1:7430343/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  1:2547441/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:4090169/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  1:3322499/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:2731165/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:6396754/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:384467/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  2:3688686/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:2542141/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:4423446/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  2:7087095/93‑1 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:7220882/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:599164/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:5450894/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:5107459/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:7206717/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:4690199/93‑3 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAATAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  1:3395823/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:6834227/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:4846640/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:7500769/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:6272823/93‑4 (MQ=255)                    gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                         <  2:1286216/93‑5 (MQ=255)                      agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                        <  1:3283474/93‑1 (MQ=255)                      agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                        <  2:3678096/93‑1 (MQ=255)                       gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:972897/93‑1 (MQ=255)                       gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:1086984/93‑1 (MQ=255)                       gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:608877/93‑1 (MQ=255)                       gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:3269379/93‑1 (MQ=255)                        tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  1:4558408/93‑1 (MQ=255)                        tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  2:1632553/93‑1 (MQ=255)                        tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  2:1307840/93‑1 (MQ=255)                         aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                                     <  2:4862367/93‑1 (MQ=255)                          tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  2:2574690/93‑1 (MQ=255)                          tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  2:1820286/93‑1 (MQ=255)                          tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  1:4461691/93‑1 (MQ=255)                          tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  1:6527821/93‑1 (MQ=255)                           gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                   <  2:1797766/93‑1 (MQ=255)                            cTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                  <  2:2191865/93‑1 (MQ=255)                             tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  2:3065511/93‑1 (MQ=255)                             tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  1:4241208/93‑1 (MQ=255)                             tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  1:3538297/93‑1 (MQ=255)                             tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  2:3834225/93‑1 (MQ=255)                              gATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaa                                                <  2:5122377/93‑1 (MQ=255)                                  cAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatata                                            <  1:3879526/93‑1 (MQ=255)                                     tGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataaat                                         >  1:2332700/1‑93 (MQ=255)                                      gAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                        <  2:3374021/93‑1 (MQ=255)                                          tatCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACAtt                                    <  1:3988715/93‑1 (MQ=255)                                           atCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTc                                   <  1:7029741/93‑1 (MQ=255)                                                      aaaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCAt                        >  1:560180/1‑93 (MQ=255)                                                                      tAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACatt        <  1:6773556/93‑1 (MQ=255)                                                                            aCATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCa  >  1:6590349/1‑93 (MQ=255)                                                                                           |                                                                              ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCA  >  NC_000913/4041564‑4041730

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