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

ATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAG  >  NC_000913/167607‑167877                                                                                                                                         |                                                                                                                                       aTCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGAtcg                                                                                                                                      >  1:866795/1‑139 (MQ=255)            cGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATg                                                                                                                           >  1:276395/1‑139 (MQ=255)              cTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgc                                                                                                                         >  1:597429/1‑139 (MQ=255)                   ccTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGc                                                                                                                    >  1:889887/1‑139 (MQ=255)                       cgccgcAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCa                                                                                                                >  2:1081797/1‑139 (MQ=255)                                                            tGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCg                                                                           >  1:709695/1‑139 (MQ=255)                                                             gcggcgCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCgg                                                                          >  2:271621/1‑139 (MQ=255)                                                                 cgcgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGTCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      <  1:718915/139‑1 (MQ=255)                                                                   cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcgc                                                                                                                       <  2:66467/88‑1 (MQ=255)                                                                   cgACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGgcgc                                                                                                                       >  1:66467/1‑88 (MQ=255)                                                                                ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGcc                                                                            <  2:63426/118‑1 (MQ=255)                                                                                ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGcc                                                                            >  1:63426/1‑118 (MQ=255)                                                                                ccGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTAcgcg                                                       >  1:499051/1‑139 (MQ=255)                                                                                  ggCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTg                                                     <  2:370437/139‑1 (MQ=255)                                                                                   gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACCGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTac                                                                                                  >  2:1092626/1‑92 (MQ=255)                                                                                   gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGc                                                                                                                    >  1:1057414/1‑75 (MQ=255)                                                                                   gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGc                                                                                                                    <  2:1057414/75‑1 (MQ=255)                                                                                   gCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTaa                                                                                                  <  1:1092626/93‑1 (MQ=255)                                                                                    cACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgc                                                   <  1:271621/139‑1 (MQ=255)                                                                                       tAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                                >  1:187975/1‑139 (MQ=255)                                                                                       tAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCa                                                >  1:116080/1‑139 (MQ=255)                                                                                           aCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCa                                            <  2:276395/139‑1 (MQ=255)                                                                                           aCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCa                                            <  2:116080/139‑1 (MQ=255)                                                                                                tACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGat                                                                                                          >  2:188169/1‑71 (MQ=255)                                                                                                tACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGaa                                                                                                          <  1:188169/72‑1 (MQ=255)                                                                                                   gCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                    >  1:944946/1‑107 (MQ=255)                                                                                                   gCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGtct                                                                    <  2:944946/107‑1 (MQ=255)                                                                                                   gCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                    <  2:187975/139‑1 (MQ=255)                                                                                                    ccGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      <  1:1019103/104‑1 (MQ=255)                                                                                                    ccGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      >  2:1019103/1‑104 (MQ=255)                                                                                                       aTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGTCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGacc                                >  2:729245/1‑139 (MQ=255)                                                                                                       aTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                  <  2:314795/137‑1 (MQ=255)                                                                                                       aTTCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                  >  1:314795/1‑137 (MQ=255)                                                                                                        ttCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAg                                                                                                               <  2:660025/59‑1 (MQ=255)                                                                                                        ttCAAAAAGTGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAg                                                                                                               >  1:660025/1‑59 (MQ=255)                                                                                                                 tGCCACAGTCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCAt                      >  1:305897/1‑139 (MQ=255)                                                                                                                         tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                    <  2:965423/117‑1 (MQ=255)                                                                                                                         tCTATTTCTGTTGTGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAt                                    >  1:965423/1‑117 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCt                                                                                        <  1:950111/55‑1 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCt                                                                                        >  2:950111/1‑55 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                  >  2:1023413/1‑93 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGgcg                                                  <  1:1023413/93‑1 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa    <  2:362680/139‑1 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa    >  1:362680/1‑139 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGa    >  2:736121/1‑139 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCAGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      >  2:1042688/1‑73 (MQ=255)                                                                                                                                   ttgtGATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCAGTAAAAGAAGCGCTTAGCTACACGCCGGgtgt                                                                      <  1:1042688/73‑1 (MQ=255)                                                                                                                                     gtgATCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAg  >  2:78216/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       ATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAG  >  NC_000913/167607‑167877

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