Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,342,380	A→C	L78V (TTG→GTG)	USA300HOU_2199 ←	possible HAD superfamily hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,342,380	0	A	C	100.0%	187.4 / NA	49	L78V (TTG→GTG)	USA300HOU_2199	possible HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (25/24);  total (25/24)

TTCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCAAATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTCTCTTA  >  CP000730/2342240‑2342526                                                                                                                                             |                                                                                                                                                   ttCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAg                                                                                                                                          <  2:217843/151‑1 (MQ=255) ttCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAg                                                                                                                                          <  2:227541/151‑1 (MQ=255)  tCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGa                                                                                                                                         <  2:220244/151‑1 (MQ=255)  tCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGa                                                                                                                                         <  2:362447/151‑1 (MQ=255)   cGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGat                                                                                                                                        <  1:313294/151‑1 (MQ=255)         tataTTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTaa                                                                                                                                  >  2:132658/1‑151 (MQ=255)           tatTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATg                                                                                                                                >  2:441579/1‑151 (MQ=255)            atTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATgg                                                                                                                               >  2:419645/1‑151 (MQ=255)            atTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATgg                                                                                                                               >  2:187716/1‑151 (MQ=255)            atTGCGCAGAAATATGACGTAATATATCTTCAATAACTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCTCATCAGAAAAGATATGTAATgg                                                                                                                               >  2:419839/1‑151 (MQ=255)                cgcAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTa                                                                                                                           >  2:255372/1‑151 (MQ=255)                 gcAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTaa                                                                                                                          >  1:145525/1‑151 (MQ=255)                 gcAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGGAAGGGTTTaa                                                                                                                          >  1:193185/1‑151 (MQ=255)                        tatGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCat                                                                                                                   <  1:102365/151‑1 (MQ=255)                          tGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCatat                                                                                                                 >  2:196352/1‑151 (MQ=255)                          tGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCatat                                                                                                                 >  1:450925/1‑151 (MQ=255)                              gTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAAc                                                                                                             <  1:306067/151‑1 (MQ=255)                                   atatcttcaatatcttcaAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGAt                                                                                                        <  1:36605/151‑1 (MQ=255)                                    tatcttcaatatcttcaaATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATc                                                                                                       >  1:224962/1‑151 (MQ=255)                                     atcttcaatatcttcaaATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATcc                                                                                                      >  1:236827/1‑151 (MQ=255)                                                     aaTAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGacac                                                                                      >  2:439731/1‑151 (MQ=255)                                                         aaaaTGATTCACTCATCACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGCTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTg                                                                                  <  2:25071/151‑1 (MQ=255)                                                          aaaTGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGa                                                                                 <  1:382533/151‑1 (MQ=255)                                                             tGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAAt                                                                              >  1:25633/1‑151 (MQ=255)                                                                    cTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTc                                                                       <  2:399274/151‑1 (MQ=255)                                                                          aCTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCa                                                                 >  1:224428/1‑151 (MQ=255)                                                                           ctttctttTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCaa                                                                >  1:247287/1‑151 (MQ=255)                                                                           ctttctttTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCaa                                                                <  1:91356/151‑1 (MQ=255)                                                                           ctttctttTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCaa                                                                <  1:381693/151‑1 (MQ=255)                                                                                 ttttATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATAcagc                                                          <  2:157042/151‑1 (MQ=255)                                                                                 ttttATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATAcagc                                                          <  2:147061/151‑1 (MQ=255)                                                                                    tataTGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGt                                                       >  1:213607/1‑151 (MQ=255)                                                                                     atatGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTg                                                      <  2:226908/151‑1 (MQ=255)                                                                                         gTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCa                                                  >  2:199352/1‑151 (MQ=255)                                                                                         gTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCa                                                  >  2:71474/1‑151 (MQ=255)                                                                                         gTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCa                                                  >  2:233726/1‑151 (MQ=255)                                                                                          tttCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCAc                                                 <  2:47433/151‑1 (MQ=255)                                                                                           ttCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACt                                                >  2:124280/1‑151 (MQ=255)                                                                                                  tCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGGTCACTTCattat                                         <  1:285806/151‑1 (MQ=255)                                                                                                  tCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCattat                                         <  1:332939/151‑1 (MQ=255)                                                                                                                     tCCAAGTGAACATTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCcaa                      >  2:236854/1‑151 (MQ=255)                                                                                                                                aTTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCa           <  1:181735/151‑1 (MQ=255)                                                                                                                                aTTTTAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCa           <  1:199081/151‑1 (MQ=255)                                                                                                                                    tAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTtc       >  2:185955/1‑151 (MQ=255)                                                                                                                                    tAGAACTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTCGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTtc       >  2:71001/1‑151 (MQ=255)                                                                                                                                       aaCTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTctct    <  1:254005/151‑1 (MQ=255)                                                                                                                                       aaCTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTctct    <  1:216583/151‑1 (MQ=255)                                                                                                                                       aaCTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTctct    <  1:263138/151‑1 (MQ=255)                                                                                                                                         cTCACATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTCTCTTa  >  1:314835/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                   TTCGACATTATATTGCGCAGAAATATGACGTAATATATCTTCAATATCTTCAAATAAAAATGGTTCGCTCAAGACTTTCTTTTTATATGTTTCTGATTCAACATATTCGCCATTATTCCAAGTGAACATTTTAGAACTCAAATCAGAAAAGATATGTAATGGTTTAAACGTCATATTAACTTGATCCATTTTATCCATGACACTTTGAAATGTACTTCCTGTCAATACAGCAGTGATCACTTCATTATTATTAGTTAGTTTTTCCAACAATTTTTCCATTTCTCTTA  >  CP000730/2342240‑2342526

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