Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	873,632	A→T	6.0%	S179S (TCA→TCT)	clfA →	Clumping factor A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	873,632	0	A	T	6.0%	213.5 / 2.3	68	S179S (TCA→TCT)	clfA	Clumping factor A
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base T (2/2);  total (34/34)
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.32e-01

CAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCAGATGCAC  >  USA300TCH1516_ALE/873493‑873771                                                                                                                                            |                                                                                                                                            caaGCAATACAAATGCGGAGGAATTAGTGAATCAAAGAAAATAGGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCaa                                                                                                                                            <  1:590545/141‑1 (MQ=255) caaGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCaa                                                                                                                                            <  1:121786/141‑1 (MQ=255) caaGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCaa                                                                                                                                            <  2:488080/141‑1 (MQ=255)             aTGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAacac                                                                                                                                <  1:690905/141‑1 (MQ=255)               gCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCt                                                                                                                              <  1:761368/141‑1 (MQ=255)               gCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCt                                                                                                                              <  1:732452/141‑1 (MQ=255)                   aggaATTAGGGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATTCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCaa                                                                                                                          <  2:312435/141‑1 (MQ=255)                    ggaATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCaaa                                                                                                                         <  1:628896/141‑1 (MQ=255)                       aTTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACaa                                                                                                                      >  2:417100/1‑141 (MQ=255)                          aGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGa                                                                                                                   <  1:296817/141‑1 (MQ=255)                           ttGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGaa                                                                                                                  <  1:136779/140‑1 (MQ=255)                           gTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACTCCTTCAAACAAGGaa                                                                                                                  <  2:761685/141‑1 (MQ=255)                           gTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGaa                                                                                                                  >  2:433424/1‑141 (MQ=255)                                    aCAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAAAACCTTCAAACAATGAATCAGCTccc                                                                                                         >  2:476702/1‑140 (MQ=255)                                    aCAAGTAATGAAACGACTGCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCca                                                                                                         <  2:406285/141‑1 (MQ=255)                                       aGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATTTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACag                                                                                                      <  1:601574/141‑1 (MQ=255)                                         aaaTTAAACGCCTTCTAATGAAACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGGTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACagag                                                                                                    <  1:204818/140‑1 (MQ=255)                                         aaaTGAGACGACTTCTAATGACACTAATACAGTATCACCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATTTTTCTTAAATTAAAGTTACTACCGTTGTATCAACACCTTCAAACAATGAATCAGCTCCACagag                                                                                                    <  1:181371/140‑1 (MQ=255)                                           acgAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTa                                                                                                  <  2:415472/139‑1 (MQ=255)                                           aTGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTa                                                                                                  <  1:27316/141‑1 (MQ=255)                                                  gACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGGACAGATGc                                                                                           >  1:45692/1‑141 (MQ=255)                                                  caCTTCTAATGATACTACTACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGc                                                                                           <  2:389545/140‑1 (MQ=255)                                                   aCTTCTAATGACACTAATACATTATCATCTGTAAATTCACCGCAAAATTCTACAAAGGCGGAAAAGGTTTCAACAACGCAAGATACGTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCa                                                                                          <  1:157001/141‑1 (MQ=255)                                                       cTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCacgaa                                                                                      >  1:274385/1‑137 (MQ=255)                                                              tACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGa                                                                               >  1:116560/1‑141 (MQ=255)                                                              tACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGGACAGATGCAAATAATAAAGa                                                                               >  1:489433/1‑141 (MQ=255)                                                                  taTACAGGATCATCAGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTa                                                                           <  2:724712/140‑1 (MQ=255)                                                                    tACAGTATCATCTGTAAATTCACCTCAAAATCCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGt                                                                         <  2:699291/141‑1 (MQ=255)                                                                     ataGTCTCAACTGTCACTACACCCCAAAGTTCTACAAATGCGTAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGtt                                                                        <  2:148557/139‑1 (MQ=255)                                                                     aCAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGtt                                                                        >  1:37435/1‑141 (MQ=255)                                                                      cAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTa                                                                       >  1:66989/1‑141 (MQ=255)                                                                       aGTAGCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTaa                                                                      <  2:98541/141‑1 (MQ=255)                                                                           tcatcTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCTACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCaa                                                                  >  2:117771/1‑141 (MQ=255)                                                                           tcatcTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCTACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCaa                                                                  >  2:332742/1‑141 (MQ=255)                                                                            catcTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAg                                                                 >  2:113667/1‑141 (MQ=255)                                                                               cTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCAACAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCgg                                                              >  2:145132/1‑141 (MQ=255)                                                                                tGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGt                                                             >  1:616330/1‑141 (MQ=255)                                                                                  tAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTa                                                           >  1:583834/1‑141 (MQ=255)                                                                                    aaTTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATAAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAAt                                                         >  2:539760/1‑141 (MQ=255)                                                                                    aaTTCACCTCAAAATTCAACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACGGCAGCAACACCTTCAAACAAGGAAACAGCCCCACAGAGTACAGATGCAAGCAATAAAGATGTAGTCAATCAAGCGGTTaaa                                                         >  1:197946/1‑140 (MQ=255)                                                                                     aTTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATa                                                        >  2:403596/1‑141 (MQ=255)                                                                                      ttCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATAc                                                       >  2:524756/1‑141 (MQ=255)                                                                                       tCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACa                                                      >  2:71514/1‑141 (MQ=255)                                                                                              aaaaTTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATTGTCCGACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAATTAATCAAGAGGTTAATCCAAGTccgc                                               >  1:338086/1‑141 (MQ=255)                                                                                              aaaaTTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTgcgc                                               <  1:414726/141‑1 (MQ=255)                                                                                                aaTTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCt                                             <  1:248321/141‑1 (MQ=255)                                                                                                     tACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGCAt                                        >  2:744900/1‑141 (MQ=255)                                                                                                       cAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCCAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAACTGCGCCTAGGATga                                      >  2:36373/1‑141 (MQ=255)                                                                                                       aaaaTGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATga                                      <  1:493117/140‑1 (MQ=255)                                                                                                         aaTGCGGAAAATGTCTCAACAACGCAAGATACTCCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATgaga                                    <  2:301980/141‑1 (MQ=255)                                                                                                          aTGCGGAGAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGCCTAATACAAGTGCGCCTAGAATGagag                                   <  2:737278/141‑1 (MQ=255)                                                                                                          aTGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGagag                                   <  2:116560/141‑1 (MQ=255)                                                                                                              ggAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCAtt                               >  2:260314/1‑141 (MQ=255)                                                                                                                  aaTGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCAtttagt                           <  2:582416/141‑1 (MQ=255)                                                                                                                      tttCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCAtttagtttag                       >  2:245875/1‑141 (MQ=255)                                                                                                                      tttCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCAtttagtttag                       <  1:305262/141‑1 (MQ=255)                                                                                                                      tttCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGAGGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCAtttagtttag                       <  1:323830/141‑1 (MQ=255)                                                                                                                       ttCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGc                      <  1:286339/141‑1 (MQ=255)                                                                                                                       ttCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGATCAACAGAGTACAGATGCAAGTAATAAAGGTGTAGTTAATCAAGCGGTTAATGAAAGTACGCCTAGAATGAGAGCAtttagtttagg                      >  2:302793/1‑140 (MQ=255)                                                                                                                               cGCAAGATACTTCTACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAg              <  1:332742/141‑1 (MQ=255)                                                                                                                               cGCAAGATACTTCTACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAg              <  1:117771/141‑1 (MQ=255)                                                                                                                               cGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAg              >  2:393149/1‑141 (MQ=255)                                                                                                                                 cAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCt            >  1:438702/1‑141 (MQ=255)                                                                                                                                 cAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATAAGAGCATTTAGTTTAGCGGCAGTAGcg            >  2:331740/1‑140 (MQ=255)                                                                                                                                   aGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGACTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGc          >  2:504771/1‑141 (MQ=255)                                                                                                                                    gATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCa         <  2:11096/141‑1 (MQ=255)                                                                                                                                    gATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCa         <  1:656266/141‑1 (MQ=255)                                                                                                                                          tCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCAGATGCa   >  1:32517/1‑141 (MQ=255)                                                                                                                                          tCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCAGATGCa   >  1:592166/1‑141 (MQ=255)                                                                                                                                          tCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCAGATGCa   >  1:724857/1‑141 (MQ=255)                                                                                                                                           cAACTGAAGCAACACCATCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGCGAGCATTTAGTTTAGCGGAAGTAGCTGCAGATGCAc  >  1:383499/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                            CAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCACCTCAAAATTCTACAAATGCGGAAAATGTTTCAACAACGCAAGATACTTCAACTGAAGCAACACCTTCAAACAATGAATCAGCTCCACAGAGTACAGATGCAAGTAATAAAGATGTAGTTAATCAAGCGGTTAATACAAGTGCGCCTAGAATGAGAGCATTTAGTTTAGCGGCAGTAGCTGCAGATGCAC  >  USA300TCH1516_ALE/873493‑873771

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