Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,028,269	G→C	P126A (CCA→GCA)	USA300HOU_RS10230 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,028,269	0	G	C	75.0%	138.3 / 51.8	65	P126A (CCA→GCA)	USA300HOU_RS10230	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (6/10);  new base C (19/29);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACA  >  USA300TCH1516_ALE/2028129‑2028403                                                                                                                                             |                                                                                                                                       tataACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTgg                                                                                                                                        >  2:281192/1‑141 (MQ=255)  ctaACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGgt                                                                                                                                       >  1:935765/2‑141 (MQ=255)      ccTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTAc                                                                                                                                   >  1:801073/1‑141 (MQ=255)          cACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTg                                                                                                                               <  1:233305/141‑1 (MQ=255)           aCTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGc                                                                                                                              <  1:410559/141‑1 (MQ=255)            cTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCt                                                                                                                             >  2:413917/1‑141 (MQ=255)                aCCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCAcc                                                                                                                         <  2:264200/141‑1 (MQ=255)                 ccTAATAGGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGTATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCa                                                                                                                        <  1:15467/141‑1 (MQ=255)                  cTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCaa                                                                                                                       <  1:389435/141‑1 (MQ=255)                   tacaaTGTGTCTAACAACAATGGTATTATTTTCCTTAATGTTTAAACAATTTGTAAAGTAAAACAACAAACCGACTATTATACATTCAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCCCCaaa                                                                                                                      <  2:357682/138‑1 (MQ=255)                   tAAAATGTGTCTAACAACAATTGTATTATTTTACTTAATATTTTCATCTTGATTTAATTCAAGATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCaaa                                                                                                                      <  1:89527/141‑1 (MQ=255)                            tCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTg                                                                                                             >  1:93717/1‑141 (MQ=255)                             cTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGt                                                                                                            <  1:869175/141‑1 (MQ=255)                             cTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGt                                                                                                            <  1:360664/141‑1 (MQ=255)                             cTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGt                                                                                                            <  1:921723/141‑1 (MQ=255)                             cTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGt                                                                                                            <  2:132684/141‑1 (MQ=255)                                    caaTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCCTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTtaataa                                                                                                     <  2:111984/141‑1 (MQ=255)                                    aaaTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTtaataa                                                                                                     >  1:935048/2‑141 (MQ=255)                                     aaTGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATaaa                                                                                                    >  1:962775/1‑141 (MQ=255)                                       tGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAg                                                                                                  >  2:278098/1‑141 (MQ=255)                                        ggTATTATTGTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGa                                                                                                 >  1:867745/1‑141 (MQ=255)                                            ttattTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGt                                                                                             >  1:608479/1‑141 (MQ=255)                                                       aaTATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAAc                                                                                  >  1:954347/1‑141 (MQ=255)                                                         tatTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCt                                                                                <  1:933029/141‑1 (MQ=255)                                                            ttGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCGTAACCTCAt                                                                             <  1:869173/141‑1 (MQ=255)                                                            ttGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCAt                                                                             <  2:49503/141‑1 (MQ=255)                                                             tGCAGCTTTTGTTAAGTTAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATc                                                                            >  1:314611/1‑141 (MQ=255)                                                             tGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATc                                                                            <  1:720928/141‑1 (MQ=255)                                                              gCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAACAGATGGTGAAGTCATAACCTCATcc                                                                           >  1:458225/1‑141 (MQ=255)                                                              gCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATcc                                                                           >  2:18735/1‑141 (MQ=255)                                                               cAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCg                                                                          >  1:88768/1‑141 (MQ=255)                                                               cAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCg                                                                          <  2:334719/141‑1 (MQ=255)                                                                     ttGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCaa                                                                    <  2:137090/141‑1 (MQ=255)                                                                        ttAAGTCAATATGTAGTTCTCCAATTTCATGTTGAACTTTTGCTTATGTTTATTTATATACATTTGCTTTTTTTCTTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAacta                                                                 <  1:391392/141‑1 (MQ=255)                                                                         tAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactac                                                                <  2:158595/141‑1 (MQ=255)                                                                         tAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactac                                                                <  1:451925/141‑1 (MQ=255)                                                                         tAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactac                                                                <  1:649709/141‑1 (MQ=255)                                                                         tAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAactac                                                                >  2:336872/1‑141 (MQ=255)                                                                               aaTATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtttttt                                                          <  2:393326/141‑1 (MQ=255)                                                                               aaTATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACtttttt                                                          <  1:691590/141‑1 (MQ=255)                                                                                atatGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACttttttt                                                         >  1:132725/1‑141 (MQ=255)                                                                                atatGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTATTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACttttttt                                                         >  1:77842/1‑141 (MQ=255)                                                                                    gTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGtt                                                     <  1:908058/141‑1 (MQ=255)                                                                                    gTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGtt                                                     >  1:283431/1‑141 (MQ=255)                                                                                    gTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTCGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGACGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGtt                                                     <  1:712624/141‑1 (MQ=255)                                                                                     tAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGttt                                                    <  2:250833/141‑1 (MQ=255)                                                                                      aGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTa                                                   >  2:69661/1‑141 (MQ=255)                                                                                          cGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACt                                               >  1:765217/1‑141 (MQ=255)                                                                                           gCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACtt                                              <  1:657136/141‑1 (MQ=255)                                                                                              aaTTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTtgtg                                           <  1:675538/141‑1 (MQ=255)                                                                                               aTTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGTTCACCAAACCACTTTGGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTgtgt                                          <  2:81128/141‑1 (MQ=255)                                                                                               aTTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTgtgt                                          <  2:174868/141‑1 (MQ=255)                                                                                                  ccATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTAtt                                       >  1:125348/1‑141 (MQ=255)                                                                                                             ttGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGcac                            <  1:899340/141‑1 (MQ=255)                                                                                                              tGTACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGcaca                           <  2:169588/141‑1 (MQ=255)                                                                                                                tACATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAAt                         <  2:129673/141‑1 (MQ=255)                                                                                                                  cATCTTCCCACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATaa                       >  1:522754/1‑141 (MQ=255)                                                                                                                         ccACTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaataa                <  1:570543/141‑1 (MQ=255)                                                                                                                            cTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaataataa             >  1:857648/1‑141 (MQ=255)                                                                                                                            cTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaataataa             >  1:360286/1‑141 (MQ=255)                                                                                                                            cTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaataataa             <  1:163098/141‑1 (MQ=255)                                                                                                                            cTTATATACATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCCTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGtaataataa             >  1:857521/1‑141 (MQ=255)                                                                                                                                   aCATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACg      <  1:957494/141‑1 (MQ=255)                                                                                                                                   aCATAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACg      <  1:911186/141‑1 (MQ=255)                                                                                                                                     aTAAACTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGac    <  1:936600/141‑2 (MQ=255)                                                                                                                                       aaaCTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACa  >  2:328255/1‑141 (MQ=255)                                                                                                                                       aaaCTGCTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAAATAGTGAACTAAAACCCTCTCCCGACGTAACTACATTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACa  <  1:288471/141‑1 (MQ=255)                                                                                                                                             |                                                                                                                                       TATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTAATATTTGCAGCTTTTGTTAAGTCAATATGTAGGTCGCCAATTCCATGTTGAATTTGTACATCTTCCCACTTATATACATAAACTGGTGTACGTTGCTCACCAAACCACTTTTGTTTAATAAAAGATGGTGAAGTCATAACCTCATCCGTCGCAACTACTTTTTTTGGTTTAAACTTGTGTATTAAATAACGCACAATAAGTAATAATAAAAAGACGAACA  >  USA300TCH1516_ALE/2028129‑2028403

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