Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	169,645	A→G	Y721C (TAC→TGC)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	169,645	0	A	G	100.0%	231.2 / NA	63	Y721C (TAC→TGC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/30);  total (33/30)

GCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTAT  >  NC_000913/169521‑169779                                                                                                                             |                                                                                                                                       gCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCg                                                                                                                                      >  2:1011966/1‑127 (MQ=255) gCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgctgct                                                                                                                          >  2:1301460/1‑139 (MQ=255)    aaCTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                       <  1:1390822/139‑1 (MQ=255)      cTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAac                                                                                                                     >  2:1584613/1‑139 (MQ=255)               aGTGGGAAGTTCTACGGTGGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCCCTGCGCGTTAACAACCTGTGCGATCGAGTATGCGTCCACAGCTGCTTTAACACTTATGgc                                                                                                            <  1:858315/139‑1 (MQ=255)               aGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                            <  2:314310/139‑1 (MQ=255)                             cGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcgc                                                                                              <  1:749049/139‑1 (MQ=255)                                tCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgc                                                                                                         >  2:259532/1‑125 (MQ=255)                                tCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgc                                                                                                         <  1:259532/125‑1 (MQ=255)                                          gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGa                                                                                           <  2:240060/129‑1 (MQ=255)                                          gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGa                                                                                           >  1:240060/1‑129 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                               >  1:1055240/1‑83 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                               <  2:1055240/83‑1 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAgg                                                                                   >  1:48842/1‑127 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAgg                                                                                   <  2:48842/127‑1 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       >  2:167846/1‑139 (MQ=255)                                                    tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       >  1:1724846/1‑139 (MQ=255)                                                            ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                              >  1:313130/1‑76 (MQ=255)                                                            ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                              <  2:313130/76‑1 (MQ=255)                                                                cgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                            <  1:1702196/90‑1 (MQ=255)                                                                cgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                            >  2:1702196/1‑90 (MQ=255)                                                                cgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGtt                                                           >  1:1461392/1‑139 (MQ=255)                                                                  aGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaacc                                                                        <  1:1461911/124‑1 (MQ=255)                                                                  aGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaacc                                                                        >  2:1461911/1‑124 (MQ=255)                                                                   gTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCt                                                        >  1:130820/1‑139 (MQ=255)                                                                       gCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                       <  2:1415933/72‑1 (MQ=255)                                                                       gCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                       >  1:1415933/1‑72 (MQ=255)                                                                       gCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAtt                                                    >  1:1135698/1‑139 (MQ=255)                                                                        cATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                   >  2:537106/1‑139 (MQ=255)                                                                               ggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                   >  1:1438714/1‑132 (MQ=255)                                                                               ggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                   <  2:1438714/132‑1 (MQ=255)                                                                               ggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCACAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCtttt                                            >  1:909041/1‑139 (MQ=255)                                                                                gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCATCTGCTTTAACACTTATGGCTGCTTCTGGGACGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTg                                           <  2:1141119/139‑1 (MQ=255)                                                                                gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTg                                           >  1:1141119/1‑139 (MQ=255)                                                                                         tGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGa                                  <  2:216172/139‑1 (MQ=255)                                                                                           gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAt                                                     >  1:859479/1‑118 (MQ=255)                                                                                           gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAt                                                     <  2:859479/118‑1 (MQ=255)                                                                                                 cATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTg                          >  1:1371289/1‑139 (MQ=255)                                                                                                   tGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccgcaa                                                                    >  2:251199/1‑95 (MQ=255)                                                                                                   tGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccgcaa                                                                    <  1:251199/95‑1 (MQ=255)                                                                                                   tGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                        <  1:1584613/139‑1 (MQ=255)                                                                                                      tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcgc                                                                                              <  2:1268856/66‑1 (MQ=255)                                                                                                      tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcgc                                                                                              >  1:1268856/1‑66 (MQ=255)                                                                                                           acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaa                                                                          >  1:285835/1‑81 (MQ=255)                                                                                                           acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaa                                                                          <  2:285835/81‑1 (MQ=255)                                                                                                           acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCa                      <  2:268746/133‑1 (MQ=255)                                                                                                           acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCa                      >  1:268746/1‑133 (MQ=255)                                                                                                            ccTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGc                                       >  2:1398395/1‑115 (MQ=255)                                                                                                            ccTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGc                                       <  1:1398395/115‑1 (MQ=255)                                                                                                              tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTcc                                                              <  2:664195/90‑1 (MQ=255)                                                                                                              tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTcc                                                              >  1:664195/1‑90 (MQ=255)                                                                                                               gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       <  1:1463996/80‑1 (MQ=255)                                                                                                               gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       >  2:1463996/1‑80 (MQ=255)                                                                                                                ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       >  2:463421/1‑79 (MQ=255)                                                                                                                ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                       <  1:463421/79‑1 (MQ=255)                                                                                                                ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGtt           >  2:130907/1‑139 (MQ=255)                                                                                                                  cGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCACAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg         <  2:909041/139‑1 (MQ=255)                                                                                                                      cGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACTTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGt     <  2:1724846/139‑1 (MQ=255)                                                                                                                       gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAt                                                     <  1:943711/90‑1 (MQ=255)                                                                                                                       gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAt                                                     >  2:943711/1‑90 (MQ=255)                                                                                                                        tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg         <  1:889709/133‑1 (MQ=255)                                                                                                                        tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg         >  2:889709/1‑133 (MQ=255)                                                                                                                         gAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTAt  <  2:1135698/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                       GCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTAT  >  NC_000913/169521‑169779

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