Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,894,258	T→C	6.9%	N19S (AAT→AGT)	USA300HOU_1755 ←	lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,894,258	0	T	C	6.9%	169.6 / 4.0	58	N19S (AAT→AGT)	USA300HOU_1755	lysophospholipase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

ACGTTTGAAAGTTTTCAATTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGTATTATTGT  >  CP000730/1894118‑1894398                                                                                                                                             |                                                                                                                                             aCGTTTGAAAGTTTTCAATTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAatatat                                                                                                                                              >  2:4999/1‑141 (MQ=255)      tGAAAGTTTTCAATTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGa                                                                                                                                         >  2:372027/1‑141 (MQ=255)            ttttCAATTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTGTTCTTAAGATCGTAATAACATATGGATATCTTCCCGTATGTTCTAATATATTATGACGAATa                                                                                                                                   >  2:401740/1‑141 (MQ=255)                  aTTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAacg                                                                                                                             >  1:26245/1‑141 (MQ=255)                       tCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGcc                                                                                                                        <  1:106689/141‑1 (MQ=255)                        ccTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATACTATGAGCAATAACAACGACGCCt                                                                                                                       >  1:62418/1‑141 (MQ=255)                        ccTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATACTATGAGCAATAACAACGACGCCt                                                                                                                       >  1:62503/1‑141 (MQ=255)                            tattaGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCg                                                                                                                   <  1:219237/141‑1 (MQ=255)                             attaGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGc                                                                                                                  <  1:562576/141‑1 (MQ=255)                             attaGCGCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTGCTTAACAGCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGc                                                                                                                  <  2:725364/141‑1 (MQ=255)                               taGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGcat                                                                                                                <  1:695240/141‑1 (MQ=255)                                aGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGcatc                                                                                                               >  1:571606/1‑141 (MQ=255)                                  ctctTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGcatcat                                                                                                             >  2:537695/1‑141 (MQ=255)                                         aGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGc                                                                                                      <  1:305263/141‑1 (MQ=255)                                         aGTTTGACATTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATAGATTATGTGCAATAACAACGACGCCTTTCGCATCATTTTCAGc                                                                                                      <  1:676750/141‑1 (MQ=255)                                          gTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCt                                                                                                     >  1:534515/1‑141 (MQ=255)                                              gaccttgaccAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCAttt                                                                                                           <  1:284955/131‑1 (MQ=255)                                              gaccttgaccAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCAttt                                                                                                           >  2:284955/1‑131 (MQ=255)                                              gaccttgaccAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTccc                                                                                                 >  1:432974/1‑141 (MQ=255)                                                ccttgaccAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCAc                                                                                               <  1:636963/141‑1 (MQ=255)                                                    gaccAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTcc                                                                                           <  1:13540/141‑1 (MQ=255)                                                           cAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACattat                                                                                    <  1:785095/141‑1 (MQ=255)                                                           cAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACattat                                                                                    <  2:322053/141‑1 (MQ=255)                                                                   ccATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATACCAAAGACGAATTTAGCATCATTTTCAGATTCCAACTTCCACATTATTATACTGc                                                                            >  1:625823/1‑141 (MQ=255)                                                                    cATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGcc                                                                           <  2:550069/141‑1 (MQ=255)                                                                    cATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGcc                                                                           <  2:472987/141‑1 (MQ=255)                                                                        aTCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCtt                                                                       >  1:629496/1‑141 (MQ=255)                                                                         acacaTGATAGCCATTTCTGCTTAACCACGTAATAACATATGAATATCTTCCCGTATGTTATAACATATTATGAGCAATAACAACGGCCCCTTTCGCATCACTTTCAGCTTCCCACGTCCACATTATTATACTGCCCCttt                                                                      <  2:219857/140‑1 (MQ=255)                                                                            caTGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTc                                                                   <  2:707158/141‑1 (MQ=255)                                                                             aTGATAGCCATTTCTTCTTAACATCGTAATAACATATGTATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCa                                                                  <  2:51470/141‑1 (MQ=255)                                                                              tGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCAt                                                                 >  2:235171/1‑141 (MQ=255)                                                                               gATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCAtt                                                                <  1:76227/141‑1 (MQ=255)                                                                               gATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCAtt                                                                <  1:76505/141‑1 (MQ=255)                                                                               gATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATACTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCAtt                                                                <  2:62503/141‑1 (MQ=255)                                                                               gATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATACTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCAtt                                                                <  2:62418/141‑1 (MQ=255)                                                                                  aGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCttttatgaaaca                                                             >  2:301107/1‑133 (MQ=255)                                                                                  aGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAAt                                                             >  2:424170/1‑141 (MQ=255)                                                                                  aGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAAt                                                             <  2:433501/141‑1 (MQ=255)                                                                                   gCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATc                                                            <  1:77128/141‑1 (MQ=255)                                                                                   gCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATAATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATc                                                            >  2:704622/1‑141 (MQ=255)                                                                                    ccATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCt                                                           <  1:380980/141‑1 (MQ=255)                                                                                     cATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATCACTCtt                                                          >  2:424596/1‑141 (MQ=255)                                                                                      aTTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTc                                                         >  2:114727/1‑141 (MQ=255)                                                                                         tcttctTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAAt                                                      >  2:348008/1‑141 (MQ=255)                                                                                             cttAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACa                                                  >  2:416750/1‑141 (MQ=255)                                                                                              ttAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACAt                                                 >  1:28552/1‑141 (MQ=255)                                                                                              ttAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACAt                                                 <  1:605433/141‑1 (MQ=255)                                                                                               tAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATa                                                >  1:258778/1‑141 (MQ=255)                                                                                                    tCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATtat                                           >  1:605080/1‑141 (MQ=255)                                                                                                     cGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATGACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATtata                                          <  2:9725/141‑1 (MQ=255)                                                                                                       taataaCATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGc                                        >  1:487388/1‑141 (MQ=255)                                                                                                            aCATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGcc                                                                           >  2:660636/1‑101 (MQ=255)                                                                                                            aCATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGcc                                                                           <  1:660636/101‑1 (MQ=255)                                                                                                                     tatCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTa                          >  2:697124/1‑141 (MQ=255)                                                                                                                      atCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAg                         <  2:645106/141‑1 (MQ=255)                                                                                                                             cGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTa                  >  2:121605/1‑141 (MQ=255)                                                                                                                                  gTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTtat             <  1:348451/141‑1 (MQ=255)                                                                                                                                   ttCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTATTCTTCAATAACATAATTATAGCAAATTCACTATGTAGTTTTCTATTtata            >  2:174118/1‑141 (MQ=255)                                                                                                                                     cTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGt          <  2:124509/141‑1 (MQ=255)                                                                                                                                           tattatGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGTattatt    <  2:571606/141‑1 (MQ=255)                                                                                                                                            attatGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGTATTAttg   <  1:522138/141‑1 (MQ=255)                                                                                                                                             ttatGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGTATTAttgt  >  1:656064/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                             ACGTTTGAAAGTTTTCAATTTGTCCTTTATTAGCTCTTGAAGTTTGACCTTGACCAGGCAAATCTCCCATAATCACATGATAGCCATTTCTTCTTAACATCGTAATAACATATGCATATCTTCCCGTATGTTCTAATATATTATGAGCAATAACAACGACGCCTTTCGCATCATTTTCAGCTTCCCACTTCCACATTATTATACTGCCCCTTTTTCATTAATCTTCAATAACATAATTATAGCAAATTCACTATGTAGATTTCTATTTATAGTATTATTGT  >  CP000730/1894118‑1894398

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