Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,856,087	T→A	P80P (CCT→CCA)	USA300HOU_2698 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,856,087	0	T	A	96.9%	229.4 / ‑1.4	65	P80P (CCT→CCA)	USA300HOU_2698	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (34/29);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAT  >  CP000730/2855951‑2856221                                                                                                                                         |                                                                                                                                       agaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAAc                                                                                                                                    >  2:81912/1‑141 (MQ=255)  gaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCACCTACATTTTAAAACGGAAGTGCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   <  1:261161/141‑1 (MQ=255)          gCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACaaa                                                                                                                           <  2:508557/141‑1 (MQ=255)           cGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAAt                                                                                                                          <  1:121509/141‑1 (MQ=255)              tGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       >  1:486393/1‑141 (MQ=255)                 tCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGa                                                                                                                    <  2:444529/141‑1 (MQ=255)                  ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                   <  1:113445/141‑1 (MQ=255)                  ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                   >  1:306240/1‑141 (MQ=255)                    tAAAGCTATATCAATAGTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAttt                                                                                                                 <  1:236404/141‑1 (MQ=255)                            tatCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttg                                                                                                         >  2:551529/1‑141 (MQ=255)                             atCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATAATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTAAACAAATGATAGATTTTATAgttgt                                                                                                        >  2:348514/1‑141 (MQ=255)                              tCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAgttgtt                                                                                                       >  1:310154/1‑141 (MQ=255)                                    ttGGAACAAGCTTTATTTCAGCACCAACTACAGTTTAAAACGGAAGTTCAGTCATCATGACCGTGTAAGATAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGACAACAAAAGATAGATTTTATAGTTGTTAcagca                                                                                                 >  1:630440/1‑141 (MQ=255)                                        aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                             >  1:565614/1‑141 (MQ=255)                                        aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                             >  2:683147/1‑141 (MQ=255)                                           aaGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAg                                                                                          >  2:43141/1‑141 (MQ=255)                                              cTTTATTTCAGCACTAACTACATTTTAAAACGGAAGTTCAGTCATCATTACTTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGata                                                                                       <  1:259318/141‑1 (MQ=255)                                              cTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGata                                                                                       <  1:167919/141‑1 (MQ=255)                                                ttATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATAtt                                                                                     >  2:539970/1‑141 (MQ=255)                                                 tATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACCGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTg                                                                                    <  2:404089/141‑1 (MQ=255)                                                      cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAatta                                                                               >  2:407922/1‑141 (MQ=255)                                                      cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTa                                                                                                      <  1:150418/118‑1 (MQ=255)                                                      cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTa                                                                                                      >  2:150418/1‑118 (MQ=255)                                                      cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAatta                                                                               >  1:657670/1‑141 (MQ=255)                                                       aGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAattat                                                                              >  2:653696/1‑141 (MQ=255)                                                        gCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAattatt                                                                             <  2:539381/141‑1 (MQ=255)                                                         cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   <  1:130828/86‑1 (MQ=255)                                                         cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   >  2:130828/1‑86 (MQ=255)                                                             aaCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGtt                                                                        <  1:349572/141‑1 (MQ=255)                                                                aaagggttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                     <  1:136226/135‑1 (MQ=255)                                                                     tttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGtt                                                                >  1:139546/1‑141 (MQ=255)                                                                         aaaCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaa                                                            >  2:265642/1‑141 (MQ=255)                                                                           aCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaaca                                                          <  2:139546/141‑1 (MQ=255)                                                                           aCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaaca                                                          <  2:649056/141‑1 (MQ=255)                                                                            cGGAAGTTCAGTCATCNTTACCTTTTAAGNTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         >  1:502600/1‑141 (MQ=255)                                                                            cGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         <  1:655275/141‑1 (MQ=255)                                                                            cGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         >  1:288380/1‑141 (MQ=255)                                                                             ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  2:661903/1‑141 (MQ=255)                                                                             ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  1:275153/1‑141 (MQ=255)                                                                             ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  1:257476/1‑141 (MQ=255)                                                                               aaGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAat                                                      <  2:219636/141‑1 (MQ=255)                                                                                 gTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAattt                                                    <  1:159154/141‑1 (MQ=255)                                                                                    cAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttatt                                                 <  2:446585/141‑1 (MQ=255)                                                                                      gTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTg                                               >  2:622056/1‑141 (MQ=255)                                                                                               aCCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTccc                                      <  1:683147/141‑1 (MQ=255)                                                                                                   tttAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGc                                  >  2:194700/1‑141 (MQ=255)                                                                                                           aaCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGa                          >  1:212893/1‑141 (MQ=255)                                                                                                            aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttat                                                  <  2:621793/116‑1 (MQ=255)                                                                                                            aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttat                                                  >  1:621793/1‑116 (MQ=255)                                                                                                            aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGAc                         <  2:698170/141‑1 (MQ=255)                                                                                                              aGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACat                       >  1:256798/1‑141 (MQ=255)                                                                                                               ggAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACata                      >  1:569694/1‑141 (MQ=255)                                                                                                               ggAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACata                      <  2:275080/141‑1 (MQ=255)                                                                                                                 agaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACAtata                    >  1:300994/1‑141 (MQ=255)                                                                                                                  gaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATAc                   >  2:105780/1‑141 (MQ=255)                                                                                                                   aTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACg                  <  1:626111/141‑1 (MQ=255)                                                                                                                    tgtgTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       >  2:296615/1‑39 (MQ=255)                                                                                                                    tgtgTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       <  1:296615/39‑1 (MQ=255)                                                                                                                       gtCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACTGAcc              >  2:261467/1‑141 (MQ=255)                                                                                                                        tCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCt             >  2:266949/1‑141 (MQ=255)                                                                                                                         cATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTc            <  2:143209/141‑1 (MQ=255)                                                                                                                          aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATCGTTGTTATAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa           <  1:245560/141‑1 (MQ=255)                                                                                                                               aTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat      <  2:466229/141‑1 (MQ=255)                                                                                                                                 cATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa    >  2:478258/1‑141 (MQ=255)                                                                                                                                   tGCGCCAAAACGACAACAACAGATAGATTTTACAGTTGTTACAGCAGCATTAGATAATGCATTATTACGTTCACTGGTTACAACAGAATTTATATGGGCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAt  >  2:121908/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                                       AGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAT  >  CP000730/2855951‑2856221

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