Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,958	A→C	Q616P (CAA→CCA)	tamB →	translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,958	0	A	C	100.0%	211.2 / NA	55	Q616P (CAA→CCA)	tamB	translocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/29);  total (26/29)

CGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACG  >  NC_000913/4445812‑4446105                                                                                                                                                   |                                                                                                                                                    cgGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACATACGCGTTGAGCCAAtt                                                                                                                                                 >  2:658845/1‑151 (MQ=255)    cACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCg                                                                                                                                              <  1:856576/151‑1 (MQ=255)       ggTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGc                                                                                                                                             >  2:62209/1‑149 (MQ=255)             aGCGCCACAACTACTGGCAGATATCACCGCGCGTGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGAt                                                                                                                                     >  2:860962/1‑151 (MQ=255)              gcgcCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATg                                                                                                                                    <  1:859735/151‑1 (MQ=255)                gcCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGt                                                                                                                                   <  1:454015/150‑1 (MQ=255)                    caACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAat                                                                                                                               <  1:390624/150‑1 (MQ=255)                            ggCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTc                                                                                                                      >  1:922027/1‑151 (MQ=255)                            ggCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTc                                                                                                                      >  2:2090/1‑151 (MQ=255)                                gATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCa                                                                                                                  <  2:393546/151‑1 (MQ=255)                                  tatCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAcc                                                                                                                >  2:596314/1‑151 (MQ=255)                                         gcgcgcgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATg                                                                                                         >  1:439779/1‑151 (MQ=255)                                              cgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaaa                                                                                                    <  1:454482/151‑1 (MQ=255)                                              cgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaaa                                                                                                    <  2:907136/151‑1 (MQ=255)                                                   tGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGc                                                                                                 >  1:305332/1‑149 (MQ=255)                                                     cgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCg                                                                                             >  1:48364/1‑151 (MQ=255)                                                     cgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCg                                                                                             >  2:52181/1‑151 (MQ=255)                                                           cAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagc                                                                                       >  2:887832/1‑151 (MQ=255)                                                             ggAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagcag                                                                                     <  2:17800/151‑1 (MQ=255)                                                               aaCTTTCCGTGGCGCAGGTTCGCGTGGAAGGGGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGGACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCCGAATGCCAAAAGCAGCGAAAag                                                                                        >  1:414330/1‑146 (MQ=255)                                                                   ttCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACg                                                                                 >  2:445373/1‑149 (MQ=255)                                                                   ttCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACg                                                                                 >  2:459255/1‑149 (MQ=255)                                                                     ccGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCt                                                                             <  1:288701/151‑1 (MQ=255)                                                                      cGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCt                                                                             >  2:124343/1‑150 (MQ=255)                                                                              aGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCg                                                                    >  2:220542/1‑151 (MQ=255)                                                                               ggTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCgg                                                                   <  1:509475/151‑1 (MQ=255)                                                                                 ttCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGAt                                                                 >  1:327248/1‑151 (MQ=255)                                                                                   cgcgTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTc                                                               >  2:485104/1‑151 (MQ=255)                                                                                      gTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAgg                                                            >  1:797563/1‑151 (MQ=255)                                                                                              cGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCt                                                    >  2:677296/1‑151 (MQ=255)                                                                                                 tatCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGtc                                                 >  2:338905/1‑151 (MQ=255)                                                                                                    cAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTcc                                              <  2:815286/151‑1 (MQ=255)                                                                                                     aaaTCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCg                                             >  1:343604/1‑151 (MQ=255)                                                                                                          cACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAg                                        <  1:653401/151‑1 (MQ=255)                                                                                                              gATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGACGGATGTAAATATCAACCTCGTAACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGCCAgaat                                     <  1:729705/150‑4 (MQ=255)                                                                                                               aTCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTaa                                   <  2:183579/151‑1 (MQ=255)                                                                                                                tCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAAt                                  >  1:373930/1‑151 (MQ=255)                                                                                                                   gATCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTg                               >  1:390189/1‑151 (MQ=255)                                                                                                                    aTCGCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTg                               <  1:294907/150‑1 (MQ=255)                                                                                                                       gCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAg                           <  2:356279/151‑1 (MQ=255)                                                                                                                       gCCGGGAAACTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAg                           <  2:857434/151‑1 (MQ=255)                                                                                                                              aaCTCGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGttt                    <  1:199601/151‑1 (MQ=255)                                                                                                                                  cGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTg                  <  1:793566/149‑1 (MQ=255)                                                                                                                                  cGACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGAt                <  2:932372/151‑1 (MQ=255)                                                                                                                                   gACGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATc               <  1:231660/151‑1 (MQ=255)                                                                                                                                   gACATACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATc               <  1:658825/151‑1 (MQ=255)                                                                                                                                     cGTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGc             >  2:632521/1‑151 (MQ=255)                                                                                                                                      gTACGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCa            <  2:96289/151‑1 (MQ=255)                                                                                                                                        aCGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCaaa          <  2:459613/151‑1 (MQ=255)                                                                                                                                        aCGCGTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCaa           <  2:636057/150‑1 (MQ=255)                                                                                                                                          gcgTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAaaga        <  2:368316/151‑1 (MQ=255)                                                                                                                                          gcgTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAaaga        <  1:935767/151‑1 (MQ=255)                                                                                                                                            gTTGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAaaga        >  1:895272/1‑149 (MQ=255)                                                                                                                                              tGAGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAagaaga     <  2:48366/150‑1 (MQ=255)                                                                                                                                                aGCCAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACg  <  2:747509/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACG  >  NC_000913/4445812‑4446105

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