Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 873,442 A→G 5.1% E116G (GAA→GGA)  clfA → Clumping factor A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE873,4420AG 5.1% 230.3 / 3.2 78E116G (GAA→GGA) clfAClumping factor A
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

CGACACAAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCAC  >  USA300TCH1516_ALE/873302‑873582
                                                                                                                                             |                                                                                                                                            
cgacACAAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGa                                                                                                                                              <  1:53645/141‑1 (MQ=255)
 gacACAAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGaa                                                                                                                                             <  1:252237/141‑1 (MQ=255)
    acaAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑GTCATCATCAACAACTGCAACTACGGAAGAAGCGCCGGTAACTGGTGAAGct                                                                                                                                          >  1:599456/1‑141 (MQ=255)
     caAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGcta                                                                                                                                         >  2:387551/1‑141 (MQ=255)
      aaaCGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACCCC‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGctac                                                                                                                                        <  2:86926/141‑1 (MQ=255)
      aaaCGTGAGGTATACGAAAACATCGCCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGctac                                                                                                                                        <  2:346563/141‑1 (MQ=255)
        aCGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGctacta                                                                                                                                      <  1:523403/141‑1 (MQ=255)
             agtgaTACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAAATACGGTAGCAACGCCGGGAAGTAGTGAAGCTGATACTACg                                                                                                                                 >  1:759565/1‑141 (MQ=255)
              gtgaTACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGa                                                                                                                                <  1:358313/141‑1 (MQ=255)
                 aTACTAACACATCACCAAACACCAATAATGGCGAAACGAGTGTTGCGCAAAATCCACCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACaa                                                                                                                             <  1:598/141‑1 (MQ=255)
                   aCTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACg                                                                                                                           <  1:214900/141‑1 (MQ=255)
                   aCTAAAACATCGGCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAAGGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACg                                                                                                                           <  2:330973/141‑1 (MQ=255)
                          cATCGTAAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACATGAAACGACACA‑ATCATCATCAACAAATGCAACGACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAg                                                                                                                    >  2:486986/1‑141 (MQ=255)
                              ggcAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTaa                                                                                                                <  1:348072/139‑1 (MQ=255)
                              gTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTaa                                                                                                                >  1:769229/1‑141 (MQ=255)
                                cAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATa                                                                                                              >  1:239018/1‑141 (MQ=255)
                                   aaccTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACCGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATacac                                                                                                           <  2:477292/138‑1 (MQ=255)
                                        aataatGGCGAAACGAGTGAGACGCAAAATCCAGCCAGAAAGGAAGTTCCCCA‑TCCGTCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGca                                                                                                      <  1:363789/141‑1 (MQ=255)
                                          taatGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGcaac                                                                                                    <  1:63911/141‑1 (MQ=255)
                                          taatGGCGAAACGAGTGTGGCGCAAAATCCACCACAACAGTAAACGACACA‑ATCATCAACCACAACTCCAACTACGGAAGAAACACCGGAAAATGGTAAAGCTAATACTACGACAACCAATCAAGCAAATGCACCGGcagc                                                                                                    >  2:229654/1‑139 (MQ=255)
                                            atGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGcaacaa                                                                                                  <  2:114624/141‑1 (MQ=255)
                                            atGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAAACAAGCTAATACACAGGaaaaaa                                                                                                  >  2:462037/1‑141 (MQ=255)
                                            acggCGAAACCAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACAAA‑ATCATCATCAACAAATGCACCTACGAAAGAAACGCCGGTAACTCGTGAAGCTACTACTACGATAACGAATCAAGCTAATACACCGGcaacaa                                                                                                  <  2:221377/139‑1 (MQ=255)
                                             tGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGCACCG‑CTCACAGCATCAGCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCaacaac                                                                                                 <  1:250981/141‑1 (MQ=255)
                                             tGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGTATCAAGCTATTACACCGGCaacaac                                                                                                 <  1:387418/141‑1 (MQ=255)
                                              ggCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACt                                                                                                <  1:387551/141‑1 (MQ=255)
                                              ggCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCAAATACACCGGCaacaacc                                                                                                >  1:135723/1‑140 (MQ=255)
                                              ggCGAAACGAGTGTGGCGCAAAACCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACCACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACt                                                                                                <  2:252635/141‑1 (MQ=255)
                                                   aaCGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGGAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACtcaatc                                                                                           >  1:329505/1‑141 (MQ=255)
                                                   aaCGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGGAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACtcaatc                                                                                           >  2:328643/1‑141 (MQ=255)
                                                   aaCGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGATAACGAATCAAGCTAATACACCGGCAACAACtcaatc                                                                                           >  2:143928/1‑141 (MQ=255)
                                                   aaCGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCACCAACtcaata                                                                                           >  2:641861/1‑140 (MQ=255)
                                                     cGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACtcaatcaa                                                                                         >  1:260765/1‑141 (MQ=255)
                                                      gAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATAAACCGGCAACAACTCAATCAAg                                                                                        >  2:158956/1‑141 (MQ=255)
                                                        gtgtGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACCACTCAATCAAGCa                                                                                      >  1:42012/1‑141 (MQ=255)
                                                            ggCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATAc                                                                                  <  1:124428/141‑1 (MQ=255)
                                                            ggCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATAc                                                                                  <  1:51572/141‑1 (MQ=255)
                                                             gcgcAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACa                                                                                 <  1:501090/141‑1 (MQ=255)
                                                             gcgcAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGCAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACa                                                                                 <  2:746160/141‑1 (MQ=255)
                                                               gcAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGGAACTACGGAAGAACCGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAGAAACTCAATCGAGCAGTACaaa                                                                               >  1:217968/1‑141 (MQ=255)
                                                               gcAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACaaa                                                                               >  1:477180/1‑141 (MQ=255)
                                                                cAAAATCCAGCACAACAGGAAAAGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAAGCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAAt                                                                              <  2:340207/141‑1 (MQ=255)
                                                                   aaTCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCg                                                                           <  1:12767/141‑1 (MQ=255)
                                                                     tCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCgga                                                                         <  1:509569/141‑1 (MQ=255)
                                                                      ccAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCggag                                                                        >  1:560569/1‑141 (MQ=255)
                                                                        aGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCggagga                                                                      <  1:272629/141‑1 (MQ=255)
                                                                              acaGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAAtaggg                                                                >  2:222163/1‑137 (MQ=255)
                                                                              acaGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAAGGCGGAGGAATAAgg                                                                >  1:230481/1‑140 (MQ=255)
                                                                               caGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATACGGAGGAATAAGTg                                                               >  2:586735/1‑141 (MQ=255)
                                                                                  gAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAAt                                                            <  1:143004/141‑1 (MQ=255)
                                                                                    aaCGACACA‑AACAACAGCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCa                                                          <  2:169762/141‑1 (MQ=255)
                                                                                       gACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAAc                                                       <  1:443767/141‑1 (MQ=255)
                                                                                       gACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAAc                                                       <  1:316845/141‑1 (MQ=255)
                                                                                         caca‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACaa                                                     <  2:560569/141‑1 (MQ=255)
                                                                                         caca‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACaa                                                     >  1:271989/1‑141 (MQ=255)
                                                                                         caca‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAACACAAATGCGGAGGAATTAGTGAGTCAAACaa                                                     >  1:86335/1‑141 (MQ=255)
                                                                                         caca‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGGAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACaa                                                     >  2:386340/1‑141 (MQ=255)
                                                                                                  ccatcaACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAATCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTATTGAATCAAACAAGTAATGaa                                             <  1:543522/140‑1 (MQ=255)
                                                                                                       aacaaATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACAAATCAAGCTAATACACCGACAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAAACAAACAAGTAATGAAACGaa                                        >  2:83917/1‑140 (MQ=255)
                                                                                                       aacaaACGCAACTACGGAAGAAACGCAGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGAc                                        <  1:319485/141‑1 (MQ=255)
                                                                                                         caaATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACtt                                      <  1:536600/141‑1 (MQ=255)
                                                                                                          aaaTGCAACTACGGAAGAAACGCCGGTAACTGGTGGAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTc                                     <  1:328643/141‑1 (MQ=255)
                                                                                                          aaaTGCAACTACGGAAGAAACGCCGGTAACTGGTGGAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTc                                     <  2:329505/141‑1 (MQ=255)
                                                                                                           aaTGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATGCAACGGCAACAACACAAGCCATTAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCt                                    <  1:180681/141‑1 (MQ=255)
                                                                                                            aTGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTa                                   <  2:769229/141‑1 (MQ=255)
                                                                                                             tGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCCAGCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTaa                                  <  1:488455/141‑1 (MQ=255)
                                                                                                             tGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTaa                                  >  1:13388/1‑141 (MQ=255)
                                                                                                               cAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATg                                <  2:697331/141‑1 (MQ=255)
                                                                                                               cAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAAAGAAACGACTTCTAATg                                >  1:609829/1‑141 (MQ=255)
                                                                                                                 aCTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGGATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGAAATTAGTGAAGCAAACAAGTAATGAACCGACTTCTAAtcaa                              >  1:285501/1‑138 (MQ=255)
                                                                                                                  cTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATa                             <  1:683644/141‑1 (MQ=255)
                                                                                                                        aagaaACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATa                       <  1:286037/141‑1 (MQ=255)
                                                                                                                         agaaACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGTCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATCCTAATAc                      <  1:487804/141‑1 (MQ=255)
                                                                                                                           aaaCGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATTATACTAATaacg                    >  2:361777/1‑138 (MQ=255)
                                                                                                                           aaaCGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATCCTAGTAAAg                    >  2:76599/1‑141 (MQ=255)
                                                                                                                             aCGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTa                  >  2:435426/1‑141 (MQ=255)
                                                                                                                               cccGGTAACTGGTGAAGCTAATACTACGACAACGAATCAAGCTAATCCACTGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTatc                <  2:505444/140‑1 (MQ=255)
                                                                                                                                 cGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTatcat              >  1:666105/1‑141 (MQ=255)
                                                                                                                                    tAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTg           >  2:328456/1‑141 (MQ=255)
                                                                                                                                          gTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAAtt     >  2:206643/1‑141 (MQ=255)
                                                                                                                                             aaGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCAc  <  2:477180/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CGACACAAACGTGAGTGATACTAAAACATCGTCAAACACTAATAATGGCGAAACGAGTGTGGCGCAAAATCCAGCACAACAGGAAACGACACA‑ATCATCATCAACAAATGCAACTACGGAAGAAACGCCGGTAACTGGTGAAGCTACTACTACGACAACGAATCAAGCTAATACACCGGCAACAACTCAATCAAGCAATACAAATGCGGAGGAATTAGTGAATCAAACAAGTAATGAAACGACTTCTAATGATACTAATACAGTATCATCTGTAAATTCAC  >  USA300TCH1516_ALE/873302‑873582

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