Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,532,298	A→T	100%	T15T (ACA→ACT)	femA_3 →	Aminoacyltransferase FemA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,532,298	0	A	T	100.0%	54.4 / NA	19	T15T (ACA→ACT)	femA_3	Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (7/12);  total (7/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  USA300TCH1516_ALE/2532173‑2532433                                                                                                                              |                                                                                                                                        aTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTca                                                                                                                          >  1:187916/1‑141 (MQ=255)           gATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacacac                                                                                                                <  1:139240/141‑1 (MQ=255)              aaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAt                                                                                                             >  1:122615/1‑141 (MQ=255)                 ataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTACAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTAGACGTCGGATGAGTTCAATGCGTCTACTACAAAGCATTTTTCACATTACACACAATCAg                                                                                                          >  2:83946/1‑141 (MQ=255)                        aaTAGCATGATTAATAAAGTCAAAAAGTATAGTTCAATTCAGAGACGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCa                                                                                                   <  2:140160/141‑1 (MQ=255)                            gCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  1:137034/141‑1 (MQ=255)                                    aataaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGCTGACTTCGGATGAGTTCAATGCCTTTACTACAAAGCATTTTTTACATCACACACAATGCTCTATTCATTACaacagtaa                                                                                       >  2:210175/1‑135 (MQ=255)                                             aaaaaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAttt                                                                              <  1:189902/141‑1 (MQ=255)                                                               cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            <  2:138002/141‑1 (MQ=255)                                                                          gggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTa                                                 >  2:69718/1‑141 (MQ=255)                                                                          gggatAAAATGAACGTTGTAACGTTAACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACGCACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTa                                                 <  2:119394/141‑1 (MQ=255)                                                                                aaaTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATg                                           <  1:189880/141‑1 (MQ=255)                                                                                           ttAATGTTGAGTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCaa                                <  2:200693/140‑1 (MQ=255)                                                                                            tAACGTTGACGTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAg                               <  1:57497/141‑1 (MQ=255)                                                                                                       tCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTGCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGa                    <  2:133900/141‑1 (MQ=255)                                                                                                             gAGTTCAATGCGTTTACTACAAAGCATGTTGCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTa              <  2:219430/141‑1 (MQ=255)                                                                                                              aGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAt             <  2:187916/141‑1 (MQ=255)                                                                                                                   aaTGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGGCAATGGTAAAGTGATTGCAGGATGCTTATTGACa        >  2:13228/1‑141 (MQ=255)                                                                                                                         tttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGCTAAAGATGACAATGGTCAAGCGATTGCAGGATGCTTATTGACAGAAgca  >  1:158320/1‑141 (MQ=255)                                                                                                                              |                                                                                                                                        ATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  USA300TCH1516_ALE/2532173‑2532433

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