Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,579,597	A→G	H35R (CAT→CGT)	yhhZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,579,597	0	A	G	100.0%	256.4 / NA	80	H35R (CAT→CGT)	yhhZ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/44);  total (36/44)

TTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAA  >  NC_000913/3579508‑3579685                                                                                          |                                                                                         ttACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGa                                                                                       <  1:6228186/93‑1 (MQ=255) ttACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGa                                                                                       <  2:203838/93‑1 (MQ=255)    ccTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGa                                                                                    <  2:6720550/93‑1 (MQ=255)    ccTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGa                                                                                    <  2:5616685/93‑1 (MQ=255)    ccTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGa                                                                                    <  2:5029469/93‑1 (MQ=255)    ccTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGa                                                                                    >  2:548376/1‑93 (MQ=255)     cTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGAt                                                                                   <  2:735038/93‑1 (MQ=255)     cTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGAt                                                                                   <  1:6208498/93‑1 (MQ=255)      tGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATg                                                                                  >  1:6114550/1‑93 (MQ=255)      tGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATg                                                                                  >  1:5922682/1‑93 (MQ=255)       gACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGa                                                                                 >  2:2463623/1‑93 (MQ=255)       gACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGa                                                                                 <  1:2675031/93‑1 (MQ=255)        aCAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGAGGGCGTGAGGATGaa                                                                                <  2:3963562/93‑1 (MQ=255)         cAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGaaa                                                                               <  2:2327863/93‑1 (MQ=255)          aGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAat                                                                              >  2:128906/1‑93 (MQ=255)           gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                             <  1:6352748/93‑1 (MQ=255)           gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                             >  1:3957338/1‑93 (MQ=255)           gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                             >  1:3697565/1‑93 (MQ=255)                 ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                       <  1:1300068/93‑1 (MQ=255)                 ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                       >  2:3617392/1‑93 (MQ=255)                 ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                       <  2:6625457/93‑1 (MQ=255)                   agaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACAtt                                                                     <  2:4519475/93‑1 (MQ=255)                   agaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACAtt                                                                     >  2:1487005/1‑93 (MQ=255)                     aacaacAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATtct                                                                   <  1:1925953/93‑1 (MQ=255)                     aacaacAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATtct                                                                   >  1:6453959/1‑93 (MQ=255)                        aacaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCAc                                                                >  2:4120937/1‑93 (MQ=255)                         acaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCAct                                                               >  1:1489364/1‑93 (MQ=255)                           aaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCActct                                                             >  2:2309881/1‑93 (MQ=255)                           aaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCActct                                                             <  2:6148076/93‑1 (MQ=255)                            aGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCtt                                                            >  2:5072494/1‑93 (MQ=255)                            aGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCtt                                                            <  1:5706201/93‑1 (MQ=255)                            aGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCtt                                                            >  1:1384166/1‑93 (MQ=255)                               aaGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCCCTCTTaaa                                                         <  1:4164849/93‑1 (MQ=255)                                aGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaat                                                        >  1:3503131/1‑93 (MQ=255)                                   aTCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataat                                                     <  2:1221819/93‑1 (MQ=255)                                   aTCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataat                                                     >  2:4129657/1‑93 (MQ=255)                                   aTCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataat                                                     >  2:1011568/1‑93 (MQ=255)                                   aTCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataat                                                     <  1:641704/93‑1 (MQ=255)                                    tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                    <  2:896678/93‑1 (MQ=255)                                    tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                    <  1:3516426/93‑1 (MQ=255)                                    tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                    <  2:2764753/93‑1 (MQ=255)                                      tcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATAtt                                                  >  2:535750/1‑93 (MQ=255)                                      tcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATAtt                                                  <  2:2121173/93‑1 (MQ=255)                                        cGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaa                                                <  1:4813936/93‑1 (MQ=255)                                         gCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaat                                               >  1:3624245/1‑93 (MQ=255)                                           aGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaataa                                             <  2:3123811/93‑1 (MQ=255)                                              ttGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATAc                                          >  1:4412352/1‑93 (MQ=255)                                                gtgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgg                                        >  1:1493057/1‑93 (MQ=255)                                                gtgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgg                                        <  2:3982290/93‑1 (MQ=255)                                                gtgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgg                                        <  2:3080549/93‑1 (MQ=255)                                                 tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                       <  1:5467299/93‑1 (MQ=255)                                                 tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                       <  2:2988405/93‑1 (MQ=255)                                                  gggACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgggg                                      >  1:6240416/1‑93 (MQ=255)                                                    gACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGTATGAAATATTTACATTCTCTCTCTTAAATAATATTAATAATACGGGGCt                                    <  2:6325141/93‑1 (MQ=255)                                                      cTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTg                                  >  2:255302/1‑93 (MQ=255)                                                         cTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGtt                               >  1:3871663/1‑93 (MQ=255)                                                          tGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTc                              >  2:2617307/1‑93 (MQ=255)                                                          tGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTc                              <  1:6014610/93‑1 (MQ=255)                                                            aGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcac                            >  2:6058271/1‑93 (MQ=255)                                                             gTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcaca                           <  1:5102359/93‑1 (MQ=255)                                                              tCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAg                          >  2:1597404/1‑93 (MQ=255)                                                                tACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                        <  1:4674519/93‑1 (MQ=255)                                                                tACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                        <  1:2444810/93‑1 (MQ=255)                                                                    ggTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCAt                    >  1:3824818/1‑93 (MQ=255)                                                                       aaTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGt                 <  1:3245115/93‑1 (MQ=255)                                                                        aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta                >  1:6590922/1‑93 (MQ=255)                                                                        aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta                >  2:1780761/1‑93 (MQ=255)                                                                        aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta                <  1:5238939/93‑1 (MQ=255)                                                                         tCGTGGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGtat               <  1:3721545/93‑1 (MQ=255)                                                                          cGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGtata              <  2:1384877/93‑1 (MQ=255)                                                                           gTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATaa             <  1:4289020/93‑1 (MQ=255)                                                                           gTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATaa             <  1:5167613/93‑1 (MQ=255)                                                                             tGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACa           <  2:4507250/93‑1 (MQ=255)                                                                              ggCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAt          <  1:2300616/93‑1 (MQ=255)                                                                              gccAGAGCGGGCGTGACGATGAAATATTTACATTCTCGCTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAt          <  2:1034859/91‑1 (MQ=255)                                                                                 agagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtttt       >  2:5693524/1‑93 (MQ=255)                                                                                  gagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTg      >  2:471022/1‑93 (MQ=255)                                                                                  gagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTg      >  1:6197577/1‑93 (MQ=255)                                                                                  gagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTg      >  1:3298443/1‑93 (MQ=255)                                                                                   agCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGt     >  1:3544620/1‑93 (MQ=255)                                                                                      gggCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTaaa  <  2:1193700/93‑1 (MQ=255)                                                                                          |                                                                                         TTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAA  >  NC_000913/3579508‑3579685

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