Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,646,189	G→C	A155G (GCA→GGA)	ypdF ←	Aminopeptidase YpdF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,646,189	0	G	C	89.7%	95.4 / 2.4	30	A155G (GCA→GGA)	ypdF	Aminopeptidase YpdF
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base C (12/14);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTAT  >  USA300TCH1516_ALE/1646046‑1646313                                                                                                                                                |                                                                                                                             cAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                          <  1:213119/148‑1 (MQ=255) cAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                          <  1:269405/148‑1 (MQ=255)         ttATCACTTGCAACACCATGTGGTAATGCGCCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                  >  2:75650/1‑148 (MQ=255)         ttATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                  >  2:351297/1‑148 (MQ=255)                   cAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatata                                                                                                        <  1:183050/148‑1 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:118642/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  2:277660/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  2:102149/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:65445/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:134546/1‑148 (MQ=255)                        ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                                                                   <  2:118683/148‑1 (MQ=255)                              ggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGttt                                                                                             >  1:350482/1‑148 (MQ=255)                                            aTGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa                                                                               <  2:317028/148‑1 (MQ=255)                                                  aGATGCTACACTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCATTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                                                         >  2:388755/1‑148 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:350126/148‑1 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:148388/148‑1 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:229991/148‑1 (MQ=255)                                                                             tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTGAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:205908/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  1:320497/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  1:84488/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:345106/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:315027/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCATTTTCAGTCATGTCGACCTTTACAACAGTTAAAATATAGTCATATGTTTAATCACCAATATCAGCTGCTTTTTGAAATAAAGCAATTTCGTCAGCa                                              <  2:114177/148‑1 (MQ=255)                                                                                 tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCtt                                          >  1:332028/1‑148 (MQ=255)                                                                                               tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                            <  2:63356/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:91349/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:378593/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:139934/148‑1 (MQ=255)                                                                                                              tAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa             >  2:86100/1‑148 (MQ=255)                                                                                                                         ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttat  >  1:111642/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                             CAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTAT  >  USA300TCH1516_ALE/1646046‑1646313

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