Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	396,540	G→A	A208T (GCA→ACA)	efeM →	putative iron uptake system component EfeM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	396,540	0	G	A	100.0%	86.2 / NA	27	A208T (GCA→ACA)	efeM	putative iron uptake system component EfeM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/9);  total (18/9)

ATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGC  >  USA300TCH1516_ALE/396404‑396665                                                                                                                                         |                                                                                                                              aTTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGcacaca                                                                                                                           >  1:167438/1‑141 (MQ=255)      atatCATACCAAAATTAATGTTCCCCGGTTGGGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACacaaacaa                                                                                                                      >  1:118847/1‑141 (MQ=255)        atCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAtt                                                                                                                    >  1:39426/1‑141 (MQ=255)              ccAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAc                                                                                                              >  1:41008/1‑141 (MQ=255)              ccAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAc                                                                                                              >  1:80816/1‑141 (MQ=255)               cAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACt                                                                                                             <  1:257145/141‑1 (MQ=255)                 aaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAGGGCACACAAAAAATTTATGACTTa                                                                                                           >  1:433789/1‑141 (MQ=255)                       aTGTTACAAGGTTCTGTTGACCTATTAAGTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                                                     >  1:323403/1‑141 (MQ=255)                       aTGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                                                     >  1:167879/1‑141 (MQ=255)                        tGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAc                                                                                                    <  1:266733/141‑1 (MQ=255)                        tGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAc                                                                                                    <  1:365777/141‑1 (MQ=255)                         gTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAcc                                                                                                   >  1:272972/1‑141 (MQ=255)                                ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                                                            <  1:156734/141‑1 (MQ=255)                                     tGTTGACCTCTTAAATGCAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACGAACAGATTTATATGATTTTAAAGCGAACGTTTAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                                                                       >  1:6151/1‑141 (MQ=255)                                         gACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaaa                                                                                   >  1:422101/1‑141 (MQ=255)                                              aTTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACCTAAAGATTTATATGATTTTAAAACGAACGTTGAAGGCACACAAGAAATTTATGGCTTATTTAAACCTATTTTAGATAAAAAACATaa                                                                              >  1:155252/1‑141 (MQ=255)                                                 aaaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaa                                                                           <  1:152073/141‑1 (MQ=255)                                                  aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaaa                                                                          <  1:37421/141‑1 (MQ=255)                                                   aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                                                         >  1:399131/1‑141 (MQ=255)                                                                      aaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAAt                                                      <  1:74342/141‑1 (MQ=255)                                                                             ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTc                                               >  1:413772/1‑141 (MQ=255)                                                                                gaagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGAt                                            <  1:493856/141‑1 (MQ=255)                                                                                       aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTGAGTGATGATATCCAAATGAACTTCGATAAAGTGa                                     <  1:311375/141‑1 (MQ=255)                                                                                         aTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAAt                                   >  1:200826/1‑141 (MQ=255)                                                                                          tttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTGGTTAAACCTATTTTCGGCTCTCTCAATAAAAATCTAAGTGATGATCTGCAAGTGAACTTCGATAAAGTGAATc                                  >  1:235541/1‑141 (MQ=255)                                                                                                              tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATaaa              >  1:210568/1‑141 (MQ=255)                                                                                                                        aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg    >  1:270704/1‑141 (MQ=255)                                                                                                                          gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcggc  >  1:29435/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                              ATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGC  >  USA300TCH1516_ALE/396404‑396665

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