Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,855 A→G 7.4% intergenic (+41/‑209) 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,8550AG 7.4% 396.8 / 15.2 121intergenic (+41/‑209)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base A (60/52);  new base G (3/6);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 122 total aligned reads displayed.

TTTGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGA  >  USA300TCH1516_ALE/770766‑770945
                                                                                         |                                                                                          
tttGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGAc                                                                                           <  1:3298727/91‑1 (MQ=255)
tttGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGAc                                                                                           >  1:1778221/1‑91 (MQ=255)
  tGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACag                                                                                         >  1:1004651/1‑91 (MQ=255)
  tGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACag                                                                                         >  1:3778787/1‑91 (MQ=255)
   gATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACaga                                                                                        >  1:3901102/1‑91 (MQ=255)
   gATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACaga                                                                                        <  1:444180/91‑1 (MQ=255)
    aTTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGTAGTGTGACagag                                                                                       <  1:3181258/91‑1 (MQ=255)
     ttACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAgaga                                                                                      <  1:3791931/91‑1 (MQ=255)
      tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                     <  1:1800984/91‑1 (MQ=255)
      tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                     <  1:174941/91‑1 (MQ=255)
      tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                     <  1:174729/91‑1 (MQ=255)
       aCAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATg                                                                                    >  1:82282/1‑91 (MQ=255)
       aCAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATg                                                                                    >  1:3952550/1‑91 (MQ=255)
        cAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGa                                                                                   >  1:798652/1‑91 (MQ=255)
         aaTAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGaa                                                                                  >  1:2382181/1‑91 (MQ=255)
          aTAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGaaa                                                                                 >  1:3881393/1‑91 (MQ=255)
            aaCTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtt                                                                               <  1:3092150/91‑1 (MQ=255)
              cTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAtttt                                                                             >  1:1678237/1‑91 (MQ=255)
               tAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAttttt                                                                            >  1:1043052/1‑91 (MQ=255)
                aaaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTta                                                                           >  1:1053117/1‑91 (MQ=255)
                   ggCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaa                                                                        <  1:1068958/91‑1 (MQ=255)
                     cAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaaa                                                                      <  1:2715617/91‑1 (MQ=255)
                      aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGGGTGTGACAGAGATGAAATTTTTATAAAAt                                                                     >  1:2346438/1‑91 (MQ=255)
                      aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                     <  1:2874804/91‑1 (MQ=255)
                      aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                     <  1:2874325/91‑1 (MQ=255)
                      aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                     >  1:2722207/1‑91 (MQ=255)
                      aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                     >  1:2721877/1‑91 (MQ=255)
                       aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                    >  1:3750028/1‑91 (MQ=255)
                       aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAAATGAATTTTATAAAAAttt                                                                   <  1:579528/91‑2 (MQ=37)
                        aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                   <  1:3084588/91‑1 (MQ=255)
                        aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                   >  1:550457/1‑91 (MQ=255)
                        aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                   <  1:1010887/91‑1 (MQ=255)
                         aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                  <  1:3431038/91‑1 (MQ=255)
                         aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                  <  1:1926359/91‑1 (MQ=255)
                          ttGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAt                                                                 <  1:485071/91‑1 (MQ=255)
                           tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                >  1:746115/1‑91 (MQ=255)
                             gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                              >  1:121657/1‑91 (MQ=255)
                             gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                              >  1:930279/1‑91 (MQ=255)
                               gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                            >  1:3197752/1‑91 (MQ=255)
                               gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                            >  1:1993397/1‑91 (MQ=255)
                                cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGc                                                           >  1:2112570/1‑91 (MQ=255)
                                 aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtc                                                          <  1:2886246/91‑3 (MQ=17)
                                  gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                         <  1:3099896/91‑4 (MQ=17)
                                  gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                         <  1:3099755/91‑4 (MQ=17)
                                  gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                         <  1:934096/91‑1 (MQ=255)
                                   ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                        <  1:378467/91‑1 (MQ=255)
                                     aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                      >  1:2232223/1‑91 (MQ=255)
                                     aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                      >  1:1133654/1‑91 (MQ=255)
                                       ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                    <  1:3239499/91‑1 (MQ=255)
                                        cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTcc                                                   <  1:497038/91‑1 (MQ=14)
                                          tgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtt                                                 <  1:2934184/91‑1 (MQ=255)
                                           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATGAAATTCATTTCTGCCACACTCCttt                                                >  1:4052960/1‑91 (MQ=255)
                                           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                <  1:4027861/91‑1 (MQ=255)
                                           gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                <  1:2326469/91‑1 (MQ=255)
                                            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                               >  1:2332863/1‑91 (MQ=255)
                                            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                               >  1:1852017/1‑91 (MQ=255)
                                            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                               <  1:3189800/91‑1 (MQ=255)
                                            aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                               >  1:243119/1‑91 (MQ=255)
                                               gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                            <  1:1640639/91‑1 (MQ=255)
                                                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                           <  1:1399048/91‑1 (MQ=255)
                                                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                           >  1:3571537/1‑91 (MQ=255)
                                                aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                           >  1:2520084/1‑91 (MQ=255)
                                                   aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                        >  1:1018175/1‑91 (MQ=255)
                                                    ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                       >  1:1970447/1‑91 (MQ=255)
                                                    ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                       >  1:3935935/1‑91 (MQ=255)
                                                    ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                       >  1:1972601/1‑91 (MQ=255)
                                                     taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAt                                      <  1:3411534/91‑1 (MQ=255)
                                                      aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                     <  1:336878/91‑1 (MQ=255)
                                                      aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                     <  1:3112305/91‑1 (MQ=255)
                                                       acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                    >  1:496747/1‑91 (MQ=255)
                                                        caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATCTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                   >  1:3286612/1‑91 (MQ=255)
                                                           caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                >  1:969945/1‑91 (MQ=255)
                                                           caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                >  1:1737287/1‑91 (MQ=255)
                                                            aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                               <  1:4020534/91‑1 (MQ=255)
                                                               aaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTa                            >  1:3789443/1‑91 (MQ=255)
                                                                aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                           <  1:2245523/91‑1 (MQ=255)
                                                                aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                           <  1:2246148/91‑1 (MQ=255)
                                                                   aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                        >  1:2104325/1‑91 (MQ=255)
                                                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                      >  1:2593080/1‑91 (MQ=255)
                                                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                      >  1:510064/1‑91 (MQ=255)
                                                                     ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                      >  1:510250/1‑91 (MQ=255)
                                                                      gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                     <  1:2166830/91‑1 (MQ=255)
                                                                      gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                     <  1:2166780/91‑1 (MQ=255)
                                                                       ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                    >  1:513040/1‑91 (MQ=255)
                                                                              ttACGAGTGCGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc             <  1:3727446/91‑1 (MQ=255)
                                                                                aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt           >  1:2663868/1‑91 (MQ=255)
                                                                                aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt           <  1:1999169/91‑1 (MQ=255)
                                                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt          >  1:1452989/1‑91 (MQ=255)
                                                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt          <  1:911528/91‑1 (MQ=255)
                                                                                 ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt          >  1:1139591/1‑91 (MQ=255)
                                                                                  gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta         <  1:1178726/91‑1 (MQ=255)
                                                                                   aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat        >  1:1802078/1‑91 (MQ=14)
                                                                                   aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat        >  1:83481/1‑91 (MQ=14)
                                                                                      gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt     <  1:4020387/91‑1 (MQ=255)
                                                                                      gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt     <  1:3953783/91‑1 (MQ=255)
                                                                                       tgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt    <  1:3198840/91‑1 (MQ=255)
                                                                                       tgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt    <  1:3200977/91‑1 (MQ=255)
                                                                                        gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg   <  1:1381512/91‑1 (MQ=255)
                                                                                        gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg   <  1:1381638/91‑1 (MQ=255)
                                                                                         aCAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGa  <  1:4112789/91‑1 (MQ=255)
                                                                                         |                                                                                          
TTTGATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGA  >  USA300TCH1516_ALE/770766‑770945

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