Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,743	C→T	T87I (ACC→ATC)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,743	0	C	T	100.0%	178.8 / NA	56	T87I (ACC→ATC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/28);  total (28/28)

TTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGA  >  NC_000913/167614‑167875                                                                                                                                  |                                                                                                                                     ttACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAaga                                                                                                                             >  1:1196843/1‑139 (MQ=255)   aCCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAgaga                                                                                                                           >  2:357152/1‑139 (MQ=255)   aCCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAgaga                                                                                                                           <  1:357152/139‑1 (MQ=255)     cacTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAGAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATg                                                                                                                         >  2:775188/3‑139 (MQ=255)               gcgcCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATc                                                                                                               >  1:742895/1‑139 (MQ=255)                               cATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATgg                                                                                                                        <  1:77330/114‑1 (MQ=255)                               cATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATgg                                                                                                                        >  2:77330/1‑114 (MQ=255)                                       ccTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAgcgc                                                                                       <  2:1002677/139‑1 (MQ=255)                                        ctgctgCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCt                                                                                      <  1:433783/139‑1 (MQ=255)                                              cAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTa                                                                                >  2:795589/1‑139 (MQ=255)                                                aCTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                              >  2:1218362/1‑139 (MQ=255)                                                     tGCGGCGCGACAGTCTGGTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTACTTCTGTTGTGATCGCCGAAGAGATGGCGCTGTATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGcca                                                                         <  1:795589/139‑2 (MQ=255)                                                      gcggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCgg                                                                        <  1:973387/139‑1 (MQ=255)                                                        ggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                      <  1:151814/139‑1 (MQ=255)                                                         gcgcGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtg                                                                     >  1:169991/1‑139 (MQ=255)                                                         gcgcGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtg                                                                     >  2:425120/1‑139 (MQ=255)                                                          cgcgACAGTCTGCTAGCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTAATTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAACGAAGCGCTTAGCTGCACGCCGGgtgt                                                                    <  2:303319/139‑1 (MQ=255)                                                          cgcgACAGTCTGCTACCGGCACTAAAACGGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    <  1:775188/139‑1 (MQ=255)                                                          cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    >  2:1295004/1‑139 (MQ=255)                                                          cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGCGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                    >  1:784926/1‑139 (MQ=255)                                                              aCAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTctct                                                                <  2:1196843/139‑1 (MQ=255)                                                                        aCCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAcgc                                                      <  1:425120/139‑1 (MQ=255)                                                                           ggCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTg                                                   >  1:441526/1‑139 (MQ=255)                                                                            gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtg                                                                     <  2:935061/120‑1 (MQ=255)                                                                            gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtg                                                                     >  1:935061/1‑120 (MQ=255)                                                                            gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTgg                                                  <  2:441526/139‑1 (MQ=255)                                                                             cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                 <  2:148989/139‑1 (MQ=255)                                                                                tAAAACCGATACGCCGATTCAAAAAGCGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                              <  2:784926/139‑1 (MQ=255)                                                                                     ccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCaa                                         <  1:1295004/139‑1 (MQ=255)                                                                                       gATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAaca                                       <  2:58070/139‑1 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGaa                                                                                                                               <  2:475086/50‑1 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGaa                                                                                                                               >  1:475086/1‑50 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGTGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCgg                                                                        >  2:956637/1‑105 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGTGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCgg                                                                        <  1:956637/105‑1 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGc                                                                                  >  1:3545/1‑95 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGc                                                                                  <  2:3545/95‑1 (MQ=255)                                                                                        aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAacac                                      >  1:878500/1‑139 (MQ=255)                                                                                             ccGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATg                                 <  1:822728/139‑1 (MQ=255)                                                                                               gATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATcc                                           <  1:654540/127‑1 (MQ=255)                                                                                               gATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATcc                                           >  2:654540/1‑127 (MQ=255)                                                                                                          tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCAt                    >  2:1129520/1‑139 (MQ=255)                                                                                                                          tgttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAcgc                                                      >  2:134303/1‑89 (MQ=255)                                                                                                                          tgttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAcgc                                                      <  1:134303/89‑1 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                              >  1:863967/1‑63 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTaca                                                                              <  2:863967/63‑1 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTa                                                                                >  2:863264/1‑61 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                         >  1:792409/1‑68 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                         <  2:792409/68‑1 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGacc                              >  2:241965/1‑111 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGacc                              <  1:241965/111‑1 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa  >  2:945772/1‑139 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa  >  2:951916/1‑139 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa  >  1:1275434/1‑139 (MQ=255)                                                                                                                            ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGATAAAAGAAGCGCTTAGCTa                                                                                <  1:863264/61‑1 (MQ=255)                                                                                                                              gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                  <  2:396769/105‑1 (MQ=255)                                                                                                                              gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                  >  1:396769/1‑105 (MQ=255)                                                                                                                                  |                                                                                                                                     TTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGA  >  NC_000913/167614‑167875

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