Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,003,480	T→C	F60S (TTC→TCC)	yedK →	SOS response‑associated peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,003,480	0	T	C	100.0%	301.2 / NA	84	F60S (TTC→TCC)	yedK	SOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/43);  total (41/43)

AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566                                                                                        |                                                                                       aaCCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      <  1:47689/90‑1 (MQ=255) aaCCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      <  2:1577880/90‑1 (MQ=255)   cccATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTgg                                                                                    <  2:276676/90‑1 (MQ=255)     cATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      >  1:1804515/1‑86 (MQ=255)     cATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCt                                                                                      <  2:1804515/86‑1 (MQ=255)     cATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTgggg                                                                                  >  1:1993412/1‑90 (MQ=255)      aTTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGa                                                                                 >  1:1406610/1‑90 (MQ=255)      aTTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTACTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGa                                                                                 >  2:2882038/1‑90 (MQ=255)       ttGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGat                                                                                >  1:2153264/1‑90 (MQ=255)        tGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGata                                                                               <  1:476622/90‑1 (MQ=255)         ggCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:511675/90‑1 (MQ=255)         ggCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  2:1993412/90‑1 (MQ=255)             gATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTc                                                                          >  1:698482/1‑90 (MQ=255)             gATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTc                                                                          >  1:508763/1‑90 (MQ=255)                 cAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCgg                                                                      <  1:1210863/90‑1 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  1:2074565/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  1:2221976/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  1:772608/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  1:2292449/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  1:864016/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  2:1129352/1‑90 (MQ=255)                   aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                                    >  2:582999/1‑90 (MQ=255)                     gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  2:398922/90‑1 (MQ=255)                     gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  2:1651783/90‑1 (MQ=255)                     gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                                  <  1:337919/90‑1 (MQ=255)                       cgcgCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtg                                                                <  2:2074565/90‑1 (MQ=255)                        gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  1:2070167/1‑90 (MQ=255)                        gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  2:794819/1‑90 (MQ=255)                        gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  2:721214/1‑90 (MQ=255)                        gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                               >  2:209939/1‑90 (MQ=255)                           ccGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGAt                                                            <  2:508763/90‑1 (MQ=255)                                  ccAAAGTCCTACTGCTCAGTGAACGTGATGAACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAAccgc                                                     <  1:1990149/90‑1 (MQ=255)                                     aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                                  <  2:382058/90‑1 (MQ=255)                                     aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                                  <  1:1129352/90‑1 (MQ=255)                                      aGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCt                                                 <  2:2949440/90‑1 (MQ=255)                                           tACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTa                                            >  1:1025149/1‑90 (MQ=255)                                               gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        <  2:1240687/90‑1 (MQ=255)                                               gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        <  2:731498/90‑1 (MQ=255)                                               gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                        >  2:1661335/1‑90 (MQ=255)                                                 tCAGTGAACGTGATGAACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  2:1985816/90‑1 (MQ=255)                                                 tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  2:2952446/1‑90 (MQ=255)                                                    gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:2690171/1‑47 (MQ=255)                                                    gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:2690171/47‑1 (MQ=255)                                                    gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCT‑‑GGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTa                                 >  2:561506/1‑90 (MQ=255)                                                      gAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:954460/85‑1 (MQ=255)                                                      gAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  2:954460/1‑85 (MQ=255)                                                         cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              <  1:534415/42‑1 (MQ=255)                                                         cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                              >  2:534415/1‑42 (MQ=255)                                                         cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                              <  2:1674615/90‑1 (MQ=255)                                                         cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGGGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                              >  2:2630161/1‑90 (MQ=255)                                                           tgatgaACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAAc                            <  1:300992/90‑1 (MQ=255)                                                           tgatgaACCCCTTCATCTGGATCCGGTTTCCGGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAAc                            <  1:2026093/90‑1 (MQ=255)                                                                aaCACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCgg                       <  2:2070167/90‑1 (MQ=255)                                                                   cccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCaa                    <  2:276536/89‑1 (MQ=255)                                                                    ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAc                   >  2:1655507/1‑90 (MQ=255)                                                                    ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAc                   >  2:2900749/1‑90 (MQ=255)                                                                     cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc                  <  1:1848344/90‑1 (MQ=255)                                                                      ttCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCa                 >  1:2869060/1‑90 (MQ=255)                                                                       tCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCgg                       >  2:1742321/1‑83 (MQ=255)                                                                       tCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCgg                       <  1:1742321/83‑1 (MQ=255)                                                                         aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc              >  1:666461/1‑90 (MQ=255)                                                                         aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc              >  2:1922760/1‑90 (MQ=255)                                                                         aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc              >  2:1651296/1‑90 (MQ=255)                                                                           ttGGATCGGGTTTCCGGGGAATATGCTCCCGGATGGGGGGAAAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:2562734/89‑1 (MQ=255)                                                                           cgggATCCGGTTTCTGGGGGATATGCCCCCGGATGGGGGGACAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:1260078/88‑1 (MQ=255)                                                                           cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  2:2153264/90‑1 (MQ=255)                                                                           cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  2:2800928/90‑1 (MQ=255)                                                                           cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  1:368200/90‑1 (MQ=255)                                                                           cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  2:1394970/90‑1 (MQ=255)                                                                           cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt            <  2:864016/90‑1 (MQ=255)                                                                            tGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  2:750337/90‑1 (MQ=255)                                                                            tGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  2:772608/90‑1 (MQ=255)                                                                            tGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTa           <  1:1399489/90‑1 (MQ=255)                                                                             ggATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCa                 >  2:601536/1‑83 (MQ=255)                                                                             ggATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCa                 <  1:601536/83‑1 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAAcc                                                       >  1:2244538/1‑44 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAAcc                                                       <  2:2244538/44‑1 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      >  2:2491963/1‑61 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                                      <  1:2491963/61‑1 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg         >  2:1918550/1‑90 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg         >  2:849966/1‑90 (MQ=255)                                                                              gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg         >  1:436354/1‑90 (MQ=255)                                                                                tCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGtt       >  2:755799/1‑90 (MQ=255)                                                                                    gTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTaaa   >  2:1232581/1‑90 (MQ=255)                                                                                     tttcCTGGGGATATGCTCCCGGATGGGGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAc  <  1:1918550/90‑1 (MQ=255)                                                                                        |                                                                                       AACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003393‑2003566

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