Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,178,128	T→G	L5R (CTC→CGC)	tolC →	transport channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,178,128	0	T	G	100.0%	194.0 / NA	53	L5R (CTC→CGC)	tolC	transport channel
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/26);  total (27/26)

ATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGAAAAAATTA  >  NC_000913/3177982‑3178274                                                                                                                                                   |                                                                                                                                                   aTGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGcccc                                                                                                                                                <  2:50732/151‑1 (MQ=255)   gtgtTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCa                                                                                                                                               <  2:92094/150‑1 (MQ=255)          aTTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATc                                                                                                                                       <  2:358286/151‑1 (MQ=255)            tGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCgg                                                                                                                                     >  2:95680/1‑151 (MQ=255)                 cacaTTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGa                                                                                                                                <  2:557024/151‑1 (MQ=255)                 cacaTTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTAACGGCCTg                                                                                                                                 >  1:191320/1‑150 (MQ=255)                       aaCGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCttt                                                                                                                          <  2:117625/151‑1 (MQ=255)                       aaCGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCttt                                                                                                                          <  1:561361/151‑1 (MQ=255)                            ccTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTggg                                                                                                                     >  1:265527/1‑151 (MQ=255)                             cTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGTAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGt                                                                                                                    >  1:693437/1‑151 (MQ=255)                                tGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTc                                                                                                                  >  1:375487/1‑150 (MQ=255)                                   cACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGtt                                                                                                              >  2:468794/1‑151 (MQ=255)                                   cACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGtt                                                                                                              >  2:460377/1‑151 (MQ=255)                                     cGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCg                                                                                                            >  2:484161/1‑151 (MQ=255)                                        aaCGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTg                                                                                                         <  1:611365/151‑1 (MQ=255)                                        aaCGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTg                                                                                                         >  2:229832/1‑151 (MQ=255)                                             cAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCa                                                                                                    <  2:436298/151‑1 (MQ=255)                                                     tGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCGGGcc                                                                                                <  1:18235/147‑1 (MQ=255)                                                     tGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCg                                                                                               >  1:388612/1‑148 (MQ=255)                                                      gCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCgaga                                                                                            <  1:285398/150‑1 (MQ=255)                                                      gCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGaa                                                                                           <  2:289597/151‑1 (MQ=255)                                                        ggTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAAcc                                                                                         <  1:232526/151‑1 (MQ=255)                                                          tAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTg                                                                                       >  2:315879/1‑151 (MQ=255)                                                               ggtgcTGCTAGAATCCGCAATAATTTGACAGTTTGATCGCGCTAAATACTGCTGCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCTAGGCCGAGAACCTGATg                                                                                    <  1:367812/145‑1 (MQ=255)                                                                cTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCa                                                                                  >  2:256956/1‑150 (MQ=255)                                                                   cTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGtt                                                                              >  1:132039/1‑151 (MQ=255)                                                                     gCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGtt                                                                              <  1:205659/149‑1 (MQ=255)                                                                      cTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTAt                                                                           >  1:30551/1‑151 (MQ=255)                                                                        aGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTAt                                                                           <  1:557709/149‑1 (MQ=255)                                                                           aTCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGACCCTGATGCAAGTTTATCagc                                                                       >  2:328252/1‑150 (MQ=255)                                                                               gCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagcaagc                                                                   >  1:443512/1‑150 (MQ=255)                                                                                  ataatTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGc                                                               >  2:400466/1‑151 (MQ=255)                                                                                     atTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCtt                                                            >  2:55294/1‑151 (MQ=255)                                                                                            aGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAAcc                                                     >  1:127274/1‑151 (MQ=255)                                                                                                tGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCgg                                                  <  2:213006/150‑1 (MQ=255)                                                                                                   tCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAAtt                                              >  1:346205/1‑151 (MQ=255)                                                                                                   tCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAAtt                                              >  2:504064/1‑151 (MQ=255)                                                                                                      cgcTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCg                                           >  2:7123/1‑151 (MQ=255)                                                                                                          aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTaa                                        <  1:188591/150‑1 (MQ=255)                                                                                                          aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTaa                                        <  2:657079/150‑1 (MQ=255)                                                                                                          aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTa                                         >  1:405941/1‑149 (MQ=255)                                                                                                           aaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGt                                      >  2:237053/1‑151 (MQ=255)                                                                                                              aCTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCt                                    <  2:233113/150‑1 (MQ=255)                                                                                                              aCTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTg                                   >  2:28679/1‑151 (MQ=255)                                                                                                                       ccacAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGACCCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGAt                           <  1:374733/150‑1 (MQ=255)                                                                                                                       ccacAAGGAATGCAAATGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTGCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATc                          <  1:460481/151‑1 (MQ=255)                                                                                                                           aaGGAATGCAAATGAAGAAATTGCGCCCCATTCTTAACGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGa                      <  2:191341/151‑1 (MQ=255)                                                                                                                                     aaTGAAGAAATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCttt            >  2:220140/1‑151 (MQ=255)                                                                                                                                        gaagaaATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGaa         <  2:131965/151‑1 (MQ=255)                                                                                                                                        gaagaaATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGaa         <  1:522903/151‑1 (MQ=255)                                                                                                                                        gaagaaATTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGa          <  1:322528/150‑1 (MQ=255)                                                                                                                                            aaaTTGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGaaaaaa     <  2:465720/151‑1 (MQ=255)                                                                                                                                                tGCGCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGAAAAAATTa  >  1:19065/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAGTAACCCGGAATTGCGTAAGTCTGCCGCCGATCGTGATGCTGCCTTTGAAAAAATTA  >  NC_000913/3177982‑3178274

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