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%	146.4 / NA	46	T87I (ACC→ATC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/26);  total (20/26)

CACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAG  >  NC_000913/167609‑167877                                                                                                                                       |                                                                                                                                       cACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGAtcgcc                                                                                                                                    <  2:132045/139‑1 (MQ=255) cACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGAtcgcc                                                                                                                                    >  1:263551/1‑139 (MQ=255)      ttACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAaga                                                                                                                               >  2:736416/1‑139 (MQ=255)      ttACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAaga                                                                                                                               >  1:475976/1‑139 (MQ=255)             tGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcg                                                                                                                        <  1:269536/139‑1 (MQ=255)                    gcgcCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATc                                                                                                                 >  2:501239/1‑139 (MQ=255)                            aGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                         <  1:719341/139‑1 (MQ=255)                            aGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAg                                                                                                         <  2:475976/139‑1 (MQ=255)                                            ccTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAgcgc                                                                                         <  2:585119/139‑1 (MQ=255)                                                          tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAgaga                                                                                                                             <  1:424880/89‑1 (MQ=255)                                                          tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAgaga                                                                                                                             >  2:424880/1‑89 (MQ=255)                                                          tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                           <  2:35352/139‑1 (MQ=255)                                                          tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                           >  1:121055/1‑139 (MQ=255)                                                             ggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                        <  2:724024/139‑1 (MQ=255)                                                               cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      <  2:705279/139‑1 (MQ=255)                                                               cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      <  2:141144/139‑1 (MQ=255)                                                                              ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAcgcg                                                       <  1:190491/139‑1 (MQ=255)                                                                                  cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                   <  1:501239/139‑1 (MQ=255)                                                                                  cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                   <  1:745691/139‑1 (MQ=255)                                                                                  cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                   <  2:121055/139‑1 (MQ=255)                                                                                   aCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                  <  2:682641/139‑1 (MQ=255)                                                                                     tAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                                <  2:143741/139‑1 (MQ=255)                                                                                     tAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                                >  2:210792/1‑139 (MQ=255)                                                                                      aaaaCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCAt                                               >  2:357747/1‑139 (MQ=255)                                                                                              tACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACAcc                                       >  2:600869/1‑139 (MQ=255)                                                                                                  ccGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATg                                   <  1:371052/139‑1 (MQ=255)                                                                                                  ccGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATg                                   <  1:600869/139‑1 (MQ=255)                                                                                                   cGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                        >  1:326222/1‑101 (MQ=255)                                                                                                   cGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgt                                                                        <  2:326222/101‑1 (MQ=255)                                                                                                   cGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                  <  1:736416/139‑1 (MQ=255)                                                                                                    gATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGac                                 <  1:512228/139‑1 (MQ=255)                                                                                                     aTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCa                                                                                                                   <  1:542825/56‑1 (MQ=255)                                                                                                     aTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCa                                                                                                                   >  2:542825/1‑56 (MQ=255)                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGa                          <  2:269129/135‑1 (MQ=255)                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGa                          >  1:269129/1‑135 (MQ=255)                                                                                                               tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCAt                      >  1:557210/1‑139 (MQ=255)                                                                                                                gCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCAtt                     >  1:395200/1‑139 (MQ=255)                                                                                                                       tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcgc                                                 <  2:66132/104‑1 (MQ=255)                                                                                                                       tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcgc                                                 >  1:66132/1‑104 (MQ=255)                                                                                                                                 ttgtGCTCGCAGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGATACACGTCGGGTGACTCTGTTGGTACGCGTGGCGCATCCAACACATATGACCACCTGATCAATCGCGGCTTTGCGGCAGa    >  2:239297/1‑139 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                  <  1:367165/93‑1 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                  >  2:367165/1‑93 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                  >  2:33168/1‑109 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                  <  1:33168/109‑1 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgc                  >  1:718529/1‑125 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTcgc                  <  2:718529/125‑1 (MQ=255)                                                                                                                                 ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa    >  1:78515/1‑139 (MQ=255)                                                                                                                                   gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAg  >  2:742376/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       CACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAG  >  NC_000913/167609‑167877

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