Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	111,516	C→A	32.4%	Q196K (CAA→AAA)	USA300HOU_0110 →	possible staphylococcal tandem lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	111,516	0	C	A	32.4%	29.4 / 13.9	27	Q196K (CAA→AAA)	USA300HOU_0110	possible staphylococcal tandem lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (3/6);  total (10/17)
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.64e-01

AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT  >  CP000730/111379‑111644                                                                                                                                          |                                                                                                                                 aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                               >  2:80448/1‑141 (MQ=255)          aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAAAGGTGCCGAGTTATTCAGCTAAATATCAATTaa                                                                                                                      >  2:25091/1‑141 (MQ=255)              agaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                                                                  <  1:177008/141‑1 (MQ=255)               aaaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTATTTGCGCGATGTTGCAGTTTGAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                 <  2:212267/140‑1 (MQ=18)                    cATAAAAAAAGAAATCGAAATCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                            <  1:133086/141‑1 (MQ=255)                           aaaGAAATCGAAAACTTTAAGTTTTTTCCGCAATATAGCAATTTTAAATATTTAATTAATTATAAAGATGGAGATATATCATATAATCCGGAGGTGCCGAGTAATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                                                     >  2:157457/1‑141 (MQ=17)                             aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                                                                                                   <  1:40044/141‑1 (MQ=255)                              aaaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                                                                  <  2:156799/140‑1 (MQ=255)                               aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGGGGAATTCGGCTAAATATAGATTAACTAATGATGATGATAATGTaa                                                                                                 >  1:83140/1‑141 (MQ=18)                               aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                 <  1:111707/141‑1 (MQ=255)                               aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                 <  1:111718/141‑1 (MQ=255)                                aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                <  2:146506/141‑1 (MQ=255)                                aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                <  2:106297/141‑1 (MQ=255)                                aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATGTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCGCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                <  1:116287/141‑1 (MQ=255)                                      aaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                          >  1:205994/1‑141 (MQ=255)                                        aCTTTAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                        <  2:43575/141‑1 (MQ=255)                                         cTTTAAGTTTTTTGCGAAATATTGGAATTTAAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                       <  1:208672/141‑1 (MQ=17)                                               gTTTTTTGCGCAATATGGCAATTTTAAAGAGTTAATCAACTATAAAGATGGAGATATATCATATAATCCAGAGGTGGCGAGTTATTCTGCTAAATATCAAGTAACTAATGATGATTATAATGTAAAACAATTCCGCaaacg                                                                                 <  2:221350/141‑3 (MQ=17)                                                ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                <  1:202479/141‑1 (MQ=255)                                                         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                       >  1:72262/1‑141 (MQ=255)                                                            tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGAGGATTATAATGTAAAACAATTACGTAAAAAATATGACATACca                                                                    >  2:83168/1‑141 (MQ=255)                                                                           aaTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCTATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                     <  1:245744/140‑1 (MQ=255)                                                                                                             ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAg                   <  1:73003/141‑1 (MQ=255)                                                                                                                    agagGTGCCGAGTTATTCACCTCGCTTACAAGTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAACAGGTtcatc            >  1:183734/1‑141 (MQ=17)                                                                                                                      agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCCTCAg          >  2:168183/1‑141 (MQ=255)                                                                                                                           ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGa     >  2:162457/1‑141 (MQ=255)                                                                                                                              ggtgatgcAGCTCAATATCAAGTGACTAATGATGATTATAATGTAAAACAATTACGTAAAAGGTATGACATACCAACAAATAAAGGGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat  <  1:174106/134‑1 (MQ=14)                                                                                                                                          |                                                                                                                                 AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT  >  CP000730/111379‑111644

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