Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	473,673	G→A	14.5%	K104K (AAG→AAA)	USA300HOU_0449 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	473,673	0	G	A	14.5%	64.8 / 2.6	29	K104K (AAG→AAA)	USA300HOU_0449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (2/2);  total (13/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.88e-01

AGGAAAATGACGAGAAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAGCGCAAAAAGAAAATATTCGCAC  >  CP000730/473548‑473813                                                                                                                              |                                                                                                                                             aGGAAAATGACGAGAAATATAGAATCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATaa                                                                                                                               <  1:218257/141‑1 (MQ=255)      aaTGACGAGAAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTgg                                                                                                                          >  1:88145/1‑141 (MQ=255)           cGAGAAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGgaggag                                                                                                                     <  2:196342/141‑1 (MQ=255)           cGAGAAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAACTAAGTGACTGATAATTTGgaggag                                                                                                                     <  1:162424/141‑1 (MQ=255)            gagaAATATAGAAGAGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGgaga                                                                                                                    <  2:123195/141‑1 (MQ=255)             agaAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGAt                                                                                                                   <  2:68939/141‑1 (MQ=255)               aaaTATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAACTAAGTGACTGATAATTTGGAGGAGATTa                                                                                                                 >  1:297669/1‑141 (MQ=255)                    taGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGa                                                                                                            >  2:153568/1‑141 (MQ=255)                       aaGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACCTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGActc                                                                                                         >  2:191594/1‑141 (MQ=255)                        aGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGActct                                                                                                        <  1:338964/141‑1 (MQ=255)                              tGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGa                                                                                                  >  2:294727/1‑141 (MQ=255)                                aGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAAc                                                                                                <  1:122100/141‑1 (MQ=255)                                    aGAAATCCATGAGATAATTAGAAAGCAAGAATAGAGATTAGAAGAATATATACAAGACTTAAAAAAGTCTTATGAAAATCAGTTAAATAACCTAAGTGACTGATAATTTGGAGGAGATTAAATGTATCTAATAGAACCAGa                                                                                            <  2:79709/141‑1 (MQ=255)                                           aaTGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGAcc                                                                                     >  1:439423/1‑141 (MQ=255)                                              gagaTAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTa                                                                                  >  2:300430/1‑141 (MQ=255)                                                    aTTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATg                                                                            >  2:402873/1‑141 (MQ=255)                                                          aaGCAAGAATGGAAATTACAAGAATATAAAGAAGACATAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTAACTGATAATTTGGAGGAGATTAAATGATTCTAATAGAACCAGATATGACCTTAAGAATGCCAGat                                                                      >  2:133475/1‑141 (MQ=255)                                                                gAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGt                                                                <  2:439423/141‑1 (MQ=255)                                                                         ttATAAGAATATAAAGAAGACTTAAAGAAGTCGTATGAAAAGCAGTTAGATAAACTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGAACTTAAGAAGGCCAGATATAAGTACTACAGTa                                                       <  2:238461/141‑1 (MQ=255)                                                                             aagaaTATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGaaa                                                   >  1:155419/1‑141 (MQ=255)                                                                               gaaTATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAaca                                                 <  1:160088/141‑1 (MQ=255)                                                                                 atatAAATAAGTCGTAATAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGTTAATTTGGAGGAGATTAAATGACTCTAATAGTACTAGATATGACCTTGAGAATGCCAGATATAAGTACTACAGTAGAAACACt                                               <  2:32427/141‑1 (MQ=255)                                                                                   ataAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTa                                             >  2:232357/1‑141 (MQ=255)                                                                                            aCTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCatat                                    <  1:165132/141‑1 (MQ=255)                                                                                                   aaaGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAAt                             <  1:144761/141‑1 (MQ=255)                                                                                                         ttATGAAAAGCAGTTAGATAAACTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGAACTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAgc                       >  1:290770/1‑141 (MQ=255)                                                                                                          tATGAAGAGCAGTTAGATACGCTAGGTGACTGATAAGTTGGAGGAGATTAAATGACTCTAATAAAACTAGATATGGCCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTGAATCTCATATCTAATATGGAAgcg                      <  2:206973/141‑1 (MQ=255)                                                                                                                  gCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAGCGCAAAAAGa              >  2:354810/1‑141 (MQ=255)                                                                                                                            aaGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAGCGCAAAAAGAAAATATTCGc    <  1:191594/141‑1 (MQ=255)                                                                                                                              gCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAGCGCAAAAAGAAAATATTCGcac  <  2:304639/141‑1 (MQ=255)                                                                                                                              |                                                                                                                                             AGGAAAATGACGAGAAATATAGAAGCGTATGAGCAAGAAATCAATGAGATAATTAGAAAGCAAGAATGGAAATTAGAAGAATATAAAGAAGACTTAAAAAAGTCTTATGAAAAGCAGTTAGATAAGCTAAGTGACTGATAATTTGGAGGAGATTAAATGACTCTAATAGAACCAGATATGACCTTAAGAATGCCAGATATAAGTACTACAGTAGAAACACTTAATCTCATATCTAAAATGGAAGCGCAAAAAGAAAATATTCGCAC  >  CP000730/473548‑473813

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