Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	168,090	A→G	T203A (ACT→GCT)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	168,090	0	A	G	100.0%	90.5 / NA	26	T203A (ACT→GCT)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/15);  total (11/15)

AACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAGAGCA  >  NC_000913/167962‑168214                                                                                                                                 |                                                                                                                             aaCGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGa                                                                                                                          >  2:623190/1‑133 (MQ=255) aaCGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGa                                                                                                                          >  2:438422/1‑133 (MQ=255) aaCGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACAACCGAACCACTGGAAGGAGTTCAGTATAAAGCCGGTGCTGa                                                                                                                          >  2:210444/1‑133 (MQ=255)                tGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGAc                                                                                                           <  1:623190/133‑1 (MQ=255)                       cccGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGtttt                                                                                                    <  1:468447/133‑1 (MQ=255)                           gTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACt                                                                                                >  2:381508/1‑133 (MQ=255)                                cGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGc                                                                                           <  2:423887/133‑1 (MQ=255)                                      ttACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCg                                                                                     <  1:381508/133‑1 (MQ=255)                                           gTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGga                                                                                <  1:544965/133‑1 (MQ=255)                                                aGCAGTCCTGGCGGCCTGCTGACTATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGatgatg                                                                           <  1:236443/133‑1 (MQ=255)                                                 gCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAtgatga                                                                          <  1:438422/133‑1 (MQ=255)                                                  cAGTCCTGGCGGCCTATTGAATTTGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGAc                                                                         <  1:210444/133‑1 (MQ=255)                                                   aGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACg                                                                        <  1:374928/133‑1 (MQ=255)                                                   aGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACg                                                                        <  2:42410/133‑1 (MQ=255)                                                     tCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTGTTGCCTTTAGCGATTCATTGGGTGATGACGgt                                                                      >  2:236443/1‑133 (MQ=255)                                                          gcggcCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGAATTTTGCGTTTCGTTGGATGATGACGGTGTTTa                                                                 >  2:259307/1‑133 (MQ=255)                                                               cTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCtt                                                            <  1:573135/133‑1 (MQ=255)                                                                            tCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCgg                                               >  1:577293/1‑133 (MQ=255)                                                                               gCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCt                                            >  1:466838/1‑133 (MQ=255)                                                                                cAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCtt                                           >  1:400155/1‑133 (MQ=255)                                                                                       ccGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGtt                                    <  2:466838/133‑1 (MQ=255)                                                                                       ccGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGtt                                    <  2:577293/133‑1 (MQ=255)                                                                                                 aaCCGCTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGc                          >  1:488433/1‑133 (MQ=255)                                                                                                      cTGAAAGAAGTTCAGTTTAAAGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCcagc                     >  1:395047/1‑133 (MQ=255)                                                                                                             aagTTCAGATTAAAGCCGATGCTGACAGCCAGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAgg              <  1:259307/133‑1 (MQ=255)                                                                                                                        aaaGCCGGTGCTGACTGCCTGTTCCAGATTGGCTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAGagc   <  2:171285/133‑1 (MQ=255)                                                                                                                        aaaGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAGagc   <  2:400155/133‑1 (MQ=255)                                                                                                                         aaGCCGGTGCTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAGagca  <  2:488433/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                             AACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAGAGCA  >  NC_000913/167962‑168214

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