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	97.0%	92.8 / NA	33	A155G (GCA→GGA)	ypdF	Aminopeptidase YpdF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (16/16);  minor base A (1/0);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATT  >  USA300TCH1516_ALE/1646051‑1646314                                                                                                                                           |                                                                                                                              atTTTATCACTTGCAACCCCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGcctcct                                                                                                                             <  2:118097/141‑1 (MQ=255)      aTCACTTGCACCACCATGTGGTAATGCACCTCTATGACCAGATGGTACAATCGTATCGAATGATGATCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTCTACTTCTTTTTCACTCATCCATCCTTTTac                                                                                                                        >  2:269850/1‑141 (MQ=255)                   ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATAAATTTCCATATGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACACCACTTACAata                                                                                                           >  2:125151/1‑141 (MQ=255)                      tgtgGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACGGTTCAAatagat                                                                                                        >  1:10116/1‑141 (MQ=255)                              tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                <  2:163561/141‑1 (MQ=255)                                cACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATAttt                                                                                              >  1:282337/1‑141 (MQ=255)                                 aCCTCTATGACCATATGCTACAATCGTATGGAATGATGGTCCATCTTCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtc                                                                                             <  2:10648/141‑1 (MQ=255)                                        tGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaaca                                                                                      <  1:204534/141‑1 (MQ=255)                                          aCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAat                                                                                    >  1:339554/1‑141 (MQ=255)                                           ccAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAata                                                                                   <  2:295335/141‑1 (MQ=255)                                                 gCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTCCCTTTAATTCTTTTTCAGGACTTCCTACTTTTACAACACTCCACCTATATTAATTTCTTGAATCACCAACATTAgtc                                                                             >  2:94780/1‑139 (MQ=255)                                                   tACAATCGTATCGAATGATGGTCCATCTGCGCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATGCTTTTTCAGTCATGCCTCCTTTGACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                           <  1:88779/141‑1 (MQ=255)                                                       aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCGTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACATTTAAAATATATTCATATGGTTCATCACCAATATTAGCTGCtttt                                                                       >  2:216148/1‑141 (MQ=255)                                                                 atgatgGTCCATCTGCTCCTAATTATAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGATGCTTTTTGAATTaaat                                                             >  2:220837/1‑140 (MQ=255)                                                                           aTCTGCTCGTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTc                                                   <  2:163398/141‑1 (MQ=255)                                                                           aTCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTc                                                   <  2:101798/141‑1 (MQ=255)                                                                              tGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCATCAATATTAGCAAGAGTGTTAATTAAAGCAATTTCGTCCGc                                                >  2:138791/1‑141 (MQ=255)                                                                                      aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGa                                        <  1:236914/141‑1 (MQ=255)                                                                                            gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                  <  1:315139/141‑1 (MQ=255)                                                                                              aTTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTa                                <  1:261706/141‑1 (MQ=255)                                                                                               ttttGCTTTCTAATATTGCCTTTAATTCTTGTTCAGTCATGCCTCCTTTTACNACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTaa                               <  2:279384/141‑1 (MQ=255)                                                                                                 ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAACGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                             >  1:83912/1‑141 (MQ=255)                                                                                                 ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                             >  1:21966/1‑141 (MQ=255)                                                                                                  tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAttt                            <  2:29862/141‑1 (MQ=255)                                                                                                    cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAACGCAATTTCGTCAGCATCTTTGACGTCTCTAAttttc                          >  1:138518/1‑140 (MQ=255)                                                                                                    cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTa                          >  1:314344/1‑141 (MQ=255)                                                                                                      ttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTtctc                        >  2:46788/1‑141 (MQ=255)                                                                                                            tatTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGt                  <  1:16285/141‑1 (MQ=255)                                                                                                              ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                <  1:216148/141‑1 (MQ=255)                                                                                                               tGCCTTTAATTCGTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAtt               <  2:308332/141‑1 (MQ=255)                                                                                                                gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa              >  2:44184/1‑141 (MQ=255)                                                                                                                gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATATACAGTATTa              >  2:217972/1‑141 (MQ=255)                                                                                                                      aaTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTGAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATg        <  2:339292/141‑1 (MQ=255)                                                                                                                            ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAAAATCTTTTACTTCTCTAATTTTATCTACCATATTAGAACTGCttatt  >  2:220582/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                              ATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATT  >  USA300TCH1516_ALE/1646051‑1646314

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: ‑