Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,022	(C)5→4	coding (1383/1494 nt)	yidJ ←	sulfatase/phosphatase superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,018	0	C	.	100.0%	240.0 / NA	51	coding (1387/1494 nt)	yidJ	sulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/33);  total (18/33)

GAATAATGCCCGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATC  >  NC_000913/3856872‑3857163                                                                                                                                                   |                                                                                                                                                  gAATAATGCCCGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGcccc                                                                                                                                                <  2:350575/150‑1 (MQ=255)  aataatGCCCGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGccccat                                                                                                                                              >  2:106662/1‑149 (MQ=255)      atGCCCGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAg                                                                                                                                         <  1:522683/151‑1 (MQ=255)         cccGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAgcgc                                                                                                                                      <  2:330788/151‑1 (MQ=255)              cgGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCg                                                                                                                                 <  2:721825/151‑1 (MQ=255)                          aGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACg                                                                                                                     <  2:735911/151‑1 (MQ=255)                             gggCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCa                                                                                                                  <  2:828782/151‑1 (MQ=255)                               gCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCAca                                                                                                                <  2:750572/151‑2 (MQ=255)                                 aTTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACgg                                                                                                               >  2:254254/1‑150 (MQ=255)                                  ttATCGTCAGAGCTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGAc                                                                                                             >  2:415182/1‑151 (MQ=255)                                    aTCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGAc                                                                                                             >  2:589506/1‑149 (MQ=255)                                          aGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTc                                                                                                     >  2:609230/1‑151 (MQ=255)                                              cTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCAtt                                                                                                 >  1:741528/1‑151 (MQ=255)                                                  tGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGt                                                                                              >  1:785341/1‑150 (MQ=255)                                                   gTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTaa                                                                                            <  2:560154/151‑1 (MQ=255)                                                    ttttttCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTgg                                                                                               <  1:199884/147‑1 (MQ=255)                                                     tttttCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGCCGCAGACTCCATTGGTaa                                                                                            >  2:437127/1‑149 (MQ=255)                                                           tcctccACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGa                                                                                     >  2:769156/1‑150 (MQ=255)                                                            cctccACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCgg                                                                                      <  1:491293/148‑1 (MQ=255)                                                             ctccACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACg                                                                                  >  1:507284/1‑151 (MQ=255)                                                                 aCCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATc                                                                              >  1:250394/1‑151 (MQ=255)                                                                 aCCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATc                                                                              <  2:300577/151‑1 (MQ=255)                                                                     tCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATcgcg                                                                           >  2:334647/1‑150 (MQ=255)                                                                      cACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATcgcg                                                                           <  2:531713/149‑1 (MQ=255)                                                                      cACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGa                                                                          >  2:778541/1‑150 (MQ=255)                                                                         cccTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAAttt                                                                      >  1:666748/1‑151 (MQ=255)                                                                          ccTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAAtttt                                                                     <  1:12314/151‑1 (MQ=255)                                                                          ccTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAAtttt                                                                     >  2:320876/1‑151 (MQ=255)                                                                          ccTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAAtttt                                                                     <  1:637056/151‑1 (MQ=255)                                                                             tgtgTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGt                                                                   <  2:377059/150‑1 (MQ=255)                                                                                gtCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCACCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCAt                                                               <  2:284010/151‑1 (MQ=255)                                                                                  cGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATg                                                              <  2:603617/150‑1 (MQ=255)                                                                                  cGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGt                                                             <  2:400449/151‑1 (MQ=255)                                                                                      aGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATc                                                         <  2:463302/151‑1 (MQ=255)                                                                                       ggCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATcc                                                        <  1:820906/151‑1 (MQ=255)                                                                                           ggTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATa                                                    <  2:138093/151‑1 (MQ=255)                                                                                           ggTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATa                                                    <  2:749629/151‑1 (MQ=255)                                                                                           ggTGTCATAGTCGCGTACAACTGGCGAATAGCCATATTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAAt                                                     >  2:267918/1‑150 (MQ=255)                                                                                              gtCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAg                                                  <  2:497806/150‑1 (MQ=255)                                                                                                         gTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCAtt                                      <  2:103739/151‑1 (MQ=255)                                                                                                              aCTGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgct                                 <  1:426566/151‑1 (MQ=255)                                                                                                                tGGCGAATAGCCATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTtgctg                                >  1:612351/1‑150 (MQ=255)                                                                                                                          ccATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGc                     <  2:704628/151‑1 (MQ=255)                                                                                                                           cATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCa                    <  2:285726/151‑1 (MQ=255)                                                                                                                           cATCTTGTGGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCa                    <  2:507320/151‑1 (MQ=255)                                                                                                                                tgtgGACGTGGACGAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAAc                <  1:630651/150‑1 (MQ=255)                                                                                                                                      cgtggacgAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGatat          <  1:423017/150‑1 (MQ=255)                                                                                                                                       gtggacgAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatc         >  2:125350/1‑150 (MQ=255)                                                                                                                                          gacgAATCG‑CCTCATCCAGCGCGGCCGTGCATCTTTACGCCACTGACGCAGTCTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATatcatc      <  2:198638/150‑1 (MQ=255)                                                                                                                                          gacgAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCg     <  1:326275/151‑1 (MQ=255)                                                                                                                                          gacgAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCg     <  2:135279/151‑1 (MQ=255)                                                                                                                                            cgAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGAt   <  1:522600/151‑1 (MQ=255)                                                                                                                                             gAAACG‑CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATc  >  1:528876/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GAATAATGCCCGCCGGGTGATCATCAGCGGGCATTATCGTCAGAACTTCTGTTTTTTCTCCTCCACCTTCACCCCTTGTGTCGGTAGGCCGGTGTCATAGTCGCGTACAACTGGCGAATAGCCATCTTGTGGACGTGGACGAAACGCCCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATC  >  NC_000913/3856872‑3857163

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑