Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,883 T→C 7.5% intergenic (+69/‑181) ybaK → / → glcR Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE770,8830TC 7.5% 644.9 / 31.5 202intergenic (+69/‑181)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base T (105/81);  new base C (4/11);  total (110/92)
Fisher's exact test for biased strand distribution p-value = 3.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 217 total aligned reads displayed.

TGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAT  >  USA300TCH1516_ALE/770793‑770973
                                                                                          |                                                                                          
tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            >  1:3276541/1‑91 (MQ=255)
tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            <  2:296291/91‑1 (MQ=255)
tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            >  1:1943071/1‑91 (MQ=255)
aggAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                                            <  2:2361419/90‑1 (MQ=255)
  gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                                          >  2:54323/1‑91 (MQ=255)
   aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                                                         >  1:544407/1‑91 (MQ=255)
   aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                                                         <  1:3631791/91‑1 (MQ=255)
    gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                                        >  2:1851005/1‑91 (MQ=255)
    gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                                        >  2:3033289/1‑91 (MQ=255)
     cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgt                                                                                       <  1:3587273/91‑2 (MQ=17)
       tttATCCATGAGTTAGTAATAACAACAAAAAGTATTGGGCAAGGTTAGGAGTGTGACAGAAATGAATTTTATAAAAATTTCATCTCTGTca                                                                                     <  2:1631853/90‑1 (MQ=9)
       gTTATCCATGGATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                                     <  2:2633975/91‑4 (MQ=17)
       gTTATCCATGAATGATTAATAACACCAAAACGTCGGGGGGACGATTAGGAGCGTGACAGACATGAAATTTTTCTAAAGTCCAATTCTGCca                                                                                     >  2:2381751/1‑91 (MQ=255)
       gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                                     >  2:2378319/1‑91 (MQ=255)
       gTTATCCATGAATGATTAAAAAAAAAAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                                     <  2:3693871/91‑1 (MQ=255)
        ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcac                                                                                    <  1:1076918/91‑1 (MQ=17)
        ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                                    <  1:3813290/91‑1 (MQ=255)
        ttATCCATGAATGACTAATAACAGAAAAAAGTATGGGGCAAGAGTAGGAGTGTGTCAGAGATGAAATTTTTATAAAATTCATTGCTGCcac                                                                                    <  2:3675936/91‑1 (MQ=255)
         tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                                                   <  2:3738775/91‑1 (MQ=17)
         tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                                   <  2:3133026/91‑1 (MQ=255)
          agccATGAATGATTAATAACAACAAAAAGAATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCacac                                                                                  <  2:3628458/89‑1 (MQ=9)
          aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                                                  >  1:644699/1‑91 (MQ=255)
           tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                                                 <  2:275782/91‑1 (MQ=255)
           tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                                                 <  2:1259754/91‑1 (MQ=255)
            ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCTACACTc                                                                                >  1:2400741/1‑91 (MQ=255)
            ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                                                <  1:3533149/91‑1 (MQ=255)
            ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                                                <  2:2867210/91‑1 (MQ=255)
            ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                                                <  2:3488643/91‑1 (MQ=255)
             cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                                               <  1:3648154/91‑1 (MQ=255)
             cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                                               >  1:2340554/1‑91 (MQ=255)
               tgtatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGGGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtt                                                                             <  2:484326/88‑1 (MQ=255)
                gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGAGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                                            <  1:3135355/91‑1 (MQ=255)
                 aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                                           >  2:581727/1‑91 (MQ=255)
                 aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                                           <  1:613139/91‑1 (MQ=255)
                 aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                                           <  2:327287/91‑1 (MQ=255)
                 aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                                           <  2:1725190/91‑1 (MQ=255)
                  atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttttt                                                                          <  1:2518447/91‑1 (MQ=255)
                       taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                                     >  1:2721731/1‑91 (MQ=255)
                       taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                                     >  1:2721705/1‑91 (MQ=255)
                        aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                                                    <  2:1893630/91‑1 (MQ=255)
                         ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                                                   >  1:722449/1‑91 (MQ=255)
                            acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                                                >  2:2486191/1‑91 (MQ=255)
                              aacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                                              <  2:1972726/91‑1 (MQ=255)
                               acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                                             >  1:1256131/1‑91 (MQ=255)
                                caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                                            <  1:1317215/91‑1 (MQ=255)
                                caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                                            >  1:1824744/1‑91 (MQ=255)
                                 aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                                                           >  1:421088/1‑91 (MQ=255)
                                  accaGTATGGGGCAAGATTAGGAGTCTGACAGAGATGAAATTTTTATAACATTCATTTCTGCAACACTCCCTTTTGATTGCATTATCAttt                                                          >  2:1023494/4‑91 (MQ=255)
                                  aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                                                          <  2:1096286/91‑1 (MQ=255)
                                    aaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTa                                                        <  1:2450082/91‑1 (MQ=255)
                                    aaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTa                                                        >  2:2326445/1‑91 (MQ=255)
                                     aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                       >  1:2764868/1‑91 (MQ=255)
                                         tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                   >  1:256968/1‑91 (MQ=255)
                                         tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                   >  1:3212691/1‑91 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  <  1:2711333/91‑1 (MQ=255)
                                          ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                  <  2:2107356/91‑1 (MQ=255)
                                           gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                 >  2:742524/1‑91 (MQ=255)
                                            ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                >  2:2536479/1‑91 (MQ=255)
                                             gCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                                               >  2:663760/1‑91 (MQ=255)
                                             gCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                                               >  2:2014918/1‑91 (MQ=255)
                                              cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                              <  1:2891777/91‑1 (MQ=255)
                                               aaGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAAc                                             >  1:3096859/1‑91 (MQ=255)
                                                aGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACa                                            <  2:3068357/91‑1 (MQ=255)
                                                  aTTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGt                                          >  1:1112409/1‑91 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                         <  2:1539059/91‑1 (MQ=255)
                                                   ttAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCGTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                         <  1:346640/91‑1 (MQ=255)
                                                     aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                       >  1:3241729/1‑91 (MQ=255)
                                                      ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTGGCATTTTACGATCATAAACAGTCAtt                                      <  2:2240707/91‑1 (MQ=255)
                                                      ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                      <  1:2022146/91‑1 (MQ=255)
                                                      ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACCCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                                      >  1:1635303/1‑91 (MQ=255)
                                                        aGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                                    >  1:682169/1‑91 (MQ=255)
                                                           gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                 <  1:2589346/91‑1 (MQ=255)
                                                           gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                 >  1:876459/1‑91 (MQ=255)
                                                            tgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCTCACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt                                >  1:3318339/1‑91 (MQ=255)
                                                             gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                               <  1:2358012/91‑1 (MQ=255)
                                                              gcAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGa                              >  1:1860126/2‑91 (MQ=17)
                                                              aCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGa                              >  1:1338816/1‑91 (MQ=255)
                                                               cAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                             >  1:5295/1‑91 (MQ=18)
                                                                agagATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGTTTGTATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                            >  2:1712145/1‑91 (MQ=255)
                                                                agaaATGAATTTTATAAAAATTTCATATCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                            >  2:3655346/1‑91 (MQ=21)
                                                                 gagaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTa                           >  1:1644313/1‑91 (MQ=255)
                                                                  agaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAt                          >  2:1914428/1‑91 (MQ=255)
                                                                     tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg                       <  1:3779243/91‑1 (MQ=255)
                                                                     tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg                       >  1:1567629/1‑91 (MQ=255)
                                                                      gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa                      >  1:2525411/1‑91 (MQ=255)
                                                                      gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa                      >  1:1283089/1‑91 (MQ=255)
                                                                       aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                     <  1:1664147/91‑1 (MQ=255)
                                                                             ttATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa               <  2:3897757/91‑1 (MQ=255)
                                                                                 aaaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATg           >  2:1295926/1‑91 (MQ=255)
                                                                                 aaaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATg           <  1:1975189/91‑1 (MQ=255)
                                                                                  aaaTTCATTTCTGCCACACTCCTTTTTGTTTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTCTTTGAACATAAACATGt          >  2:2623081/1‑91 (MQ=255)
                                                                                  aaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGt          <  1:14159/91‑1 (MQ=255)
                                                                                  aaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGt          <  1:3866788/91‑1 (MQ=255)
                                                                                  aaaTTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGt          >  2:1680732/1‑91 (MQ=255)
                                                                                     ttCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAt       >  2:2153832/1‑91 (MQ=255)
                                                                                       cATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAttt     >  2:3214124/1‑91 (MQ=255)
                                                                                       cATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAttt     >  1:2097794/1‑91 (MQ=255)
                                                                                        aTTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAtttt    >  1:841860/1‑91 (MQ=255)
                                                                                         tttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTa   >  1:3826614/1‑91 (MQ=255)
                                                                                          ttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAt  >  1:1553442/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAT  >  USA300TCH1516_ALE/770793‑770973

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