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%	208.2 / NA	58	Y721C (TAC→TGC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/25);  total (33/25)

GGTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTA  >  NC_000913/169512‑169778                                                                                                                                      |                                                                                                                                      ggTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCg                                                                                                                                  >  1:450052/1‑139 (MQ=255)  gTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGc                                                                                                                                 >  1:599078/1‑139 (MQ=255)              aCTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTaa                                                                                                                     <  2:599078/139‑1 (MQ=255)               cTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAac                                                                                                                    >  2:21578/1‑139 (MQ=255)                     tAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTAt                                                                                                              >  2:561982/1‑139 (MQ=255)                        aGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                           <  1:561982/139‑1 (MQ=255)                                         tCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGa                                                                                          >  1:234989/1‑139 (MQ=255)                                         tCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATAGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGa                                                                                          >  1:869311/1‑139 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCa                                                                                            <  2:372670/135‑1 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCa                                                                                            >  1:372670/1‑135 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAAc                                                                                        >  1:1066409/1‑139 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAAc                                                                                        >  2:1136556/1‑139 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAAc                                                                                        >  2:503066/1‑139 (MQ=255)                                           gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAAc                                                                                        >  2:497521/1‑139 (MQ=255)                                                  cGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGt                                                                                 >  2:1156180/1‑139 (MQ=255)                                                   gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTc                                                                                >  1:695696/1‑139 (MQ=255)                                                      aGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                              <  2:422598/90‑1 (MQ=255)                                                      aGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                              >  1:422598/1‑90 (MQ=255)                                                                    tGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                               <  2:398733/139‑1 (MQ=255)                                                                      gcgcGAGTCGACATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgc                                                                                                        <  1:166480/96‑1 (MQ=255)                                                                      gcgcGAGTCGACATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgc                                                                                                        >  2:166480/1‑96 (MQ=255)                                                                           aGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTc                                                        >  1:828297/1‑139 (MQ=255)                                                                             tCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTa                                                      >  1:426364/1‑139 (MQ=255)                                                                                      ctggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCtt                                             >  2:1276789/1‑139 (MQ=255)                                                                                         gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgggg                                                                                                >  1:856763/1‑85 (MQ=255)                                                                                         gctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgggg                                                                                                <  2:856763/85‑1 (MQ=255)                                                                                             cAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGc                                      <  1:497521/139‑1 (MQ=255)                                                                                               aCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCAc                                    >  1:999779/1‑139 (MQ=255)                                                                                                  tGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGa                                 <  2:450052/139‑1 (MQ=255)                                                                                                    gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAtt                               >  1:237151/1‑139 (MQ=255)                                                                                                     cgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAttt                              >  1:984926/1‑139 (MQ=255)                                                                                                          cATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTACTCTTTTGGGGCACGGATTTCCGTg                         <  1:1136556/139‑1 (MQ=255)                                                                                                            tGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                       <  2:984926/139‑1 (MQ=255)                                                                                                            tGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                       <  1:1276789/139‑1 (MQ=255)                                                                                                             gTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGccc                      <  2:1066409/139‑1 (MQ=255)                                                                                                             gTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGccc                      <  1:21578/139‑1 (MQ=255)                                                                                                              ttAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCa                     <  1:822236/139‑1 (MQ=255)                                                                                                               tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAt                    <  2:426364/139‑1 (MQ=255)                                                                                                                    acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                           <  1:339198/47‑1 (MQ=255)                                                                                                                    acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                           >  2:339198/1‑47 (MQ=255)                                                                                                                    acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTcac               >  1:419082/1‑139 (MQ=255)                                                                                                                       tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                       >  2:1106401/1‑128 (MQ=255)                                                                                                                       tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                       <  1:1106401/128‑1 (MQ=255)                                                                                                                       tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg            >  1:989198/1‑139 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaa                                                                         >  2:1269456/1‑77 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaa                                                                         <  1:1269456/77‑1 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                               <  2:166729/87‑1 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                               >  1:166729/1‑87 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg            <  2:289917/138‑1 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg            >  1:289917/1‑138 (MQ=255)                                                                                                                        gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGt           >  1:902236/1‑139 (MQ=255)                                                                                                                           cgttCGTGAAAGCATCGCCAGCAGCTTTAACACTTATGGCTGATTCTCGGGCGCAGAACGTCAGGTCATTGCAACCGGAACCTTCCGTTTCTAAATTCTCTTTTGGGGAACGGATTTCCGGGCCCATTTCGCAAGTTgg        <  2:514324/136‑1 (MQ=255)                                                                                                                                 tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCa                                     >  2:1279894/1‑104 (MQ=255)                                                                                                                                 tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCa                                     <  1:1279894/104‑1 (MQ=255)                                                                                                                                 tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa  <  2:1082217/139‑1 (MQ=255)                                                                                                                                 tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa  >  1:1082217/1‑139 (MQ=255)                                                                                                                                   aaTGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTca                >  2:697540/1‑123 (MQ=255)                                                                                                                                   aaTGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTCGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTca                <  1:697540/123‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      GGTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTA  >  NC_000913/169512‑169778

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