Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,617,317	T→G	D408E (GAT→GAG)	bepA →	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,617,317	0	T	G	100.0%	196.5 / NA	52	D408E (GAT→GAG)	bepA	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/27);  total (25/27)

TGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTG  >  NC_000913/2617173‑2617464                                                                                                                                                 |                                                                                                                                                    tGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGc                                                                                                                                               >  2:559541/1‑151 (MQ=255)       gATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGc                                                                                                                                         <  1:675516/151‑1 (MQ=255)           aaTCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTggg                                                                                                                                     >  2:843234/1‑151 (MQ=255)            aTCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTggg                                                                                                                                     >  2:359021/1‑150 (MQ=255)               gTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGAttt                                                                                                                                 >  1:884086/1‑151 (MQ=255)               gTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGAttt                                                                                                                                 <  2:522568/151‑1 (MQ=255)                   gAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCt                                                                                                                              <  2:140007/150‑1 (MQ=255)                     aaaaTGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGc                                                                                                                           >  2:503931/1‑151 (MQ=255)                          gCCCGCGATTTACGCACCAATCCGGTTTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAg                                                                                                                       <  1:186243/150‑1 (MQ=255)                           cccGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGc                                                                                                                     <  1:844316/151‑1 (MQ=255)                             cgcgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCgg                                                                                                                   <  2:693344/151‑1 (MQ=255)                               cgATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGaa                                                                                                                 <  2:71198/151‑1 (MQ=255)                                   ttACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcg                                                                                                             >  1:462092/1‑151 (MQ=255)                                    tACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgc                                                                                                            >  1:16784/1‑151 (MQ=255)                                     aCGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCcgcg                                                                                                           <  1:462863/151‑1 (MQ=255)                                        cACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCt                                                                                                         >  1:892808/1‑150 (MQ=255)                                                cGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACcg                                                                                                <  1:409917/151‑1 (MQ=255)                                                             cAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgct                                                                                    <  1:788404/150‑1 (MQ=255)                                                               aCCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggc                                                                                 <  1:227220/151‑1 (MQ=255)                                                                ccTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggct                                                                                <  2:487927/151‑1 (MQ=255)                                                                 cTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAgctggctg                                                                               >  2:709854/1‑151 (MQ=255)                                                                  tGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTgc                                                                              >  2:822175/1‑151 (MQ=255)                                                                               tttGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGt                                                                 >  2:457712/1‑151 (MQ=255)                                                                                ttGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCTGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGt                                                                 <  1:732860/150‑1 (MQ=255)                                                                                   cAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATg                                                             >  1:416648/1‑151 (MQ=255)                                                                                   cAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATg                                                             <  2:396292/151‑1 (MQ=255)                                                                                    aatgCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATgc                                                            >  2:465862/4‑151 (MQ=255)                                                                                        cGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTc                                                        >  1:176531/1‑151 (MQ=255)                                                                                        cGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTc                                                        <  1:361004/151‑1 (MQ=255)                                                                                               ccACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACACCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGgac                                                 <  2:315590/151‑1 (MQ=255)                                                                                                cacaAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGgacg                                                >  2:191131/1‑151 (MQ=255)                                                                                                   aagaagCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACt                                             <  2:486448/151‑1 (MQ=255)                                                                                                              aaTATTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGcc                                   >  2:371814/1‑150 (MQ=255)                                                                                                                tatTCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCAtt                                >  2:609269/1‑151 (MQ=255)                                                                                                                   tCTGAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTc                              <  2:763735/150‑1 (MQ=255)                                                                                                                      gAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGc                            <  2:46455/149‑1 (MQ=255)                                                                                                                      gAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCt                           <  1:12831/150‑1 (MQ=255)                                                                                                                      gAATCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTg                          <  1:229426/151‑1 (MQ=255)                                                                                                                        aTCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGtt                        >  1:572363/1‑151 (MQ=255)                                                                                                                        aTCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGtt                        <  2:248671/151‑1 (MQ=255)                                                                                                                         tCGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTg                       >  1:430061/1‑151 (MQ=255)                                                                                                                          cGCTACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGa                      <  1:307191/151‑1 (MQ=255)                                                                                                                             tACACCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGt                    >  2:373194/1‑150 (MQ=255)                                                                                                                               cacCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAgc                 <  1:603337/151‑1 (MQ=255)                                                                                                                                acCTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAgcg                <  2:303997/151‑1 (MQ=255)                                                                                                                                  cTTTAATAATAAAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGcc              >  2:960/1‑151 (MQ=255)                                                                                                                                      aataataaAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGt           <  1:53237/150‑1 (MQ=255)                                                                                                                                       ataataaAGAGGACAACAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTc         >  2:626168/1‑151 (MQ=255)                                                                                                                                        taataaAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCg        >  2:635272/1‑151 (MQ=255)                                                                                                                                         aataaAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCg        <  2:349528/150‑1 (MQ=255)                                                                                                                                          ataaAGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGc       >  1:99411/1‑150 (MQ=255)                                                                                                                                             aaGAGGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAgg    >  1:94329/1‑150 (MQ=255)                                                                                                                                              agagGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTg  >  1:474931/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGAGGCGATCAATCGTCTGAAAAATGCCCGCGATTTACGCACCAATCCGGTGTTGCAGCTCAACCTGGCGAACGCTTATTTGCAAGGCGGTCAACCACAAGAAGCGGCCAATATTCTGAATCGCTACACCTTTAATAATAAAGATGACAGCAACGGCTGGGATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTG  >  NC_000913/2617173‑2617464

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