Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,866 A→T 7.6% intergenic (+52/‑198) 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,8660AT 7.6% 746.6 / 36.7 222intergenic (+52/‑198)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base A (106/99);  new base T (5/12);  total (111/111)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 224 total aligned reads displayed.

ATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAA  >  USA300TCH1516_ALE/770776‑770954
                                                                                          |                                                                                        
    cTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGAATAGGAGTGTGACAGAGATGAAAtttt                                                                                      <  2:2247788/91‑1 (MQ=255)
    cTAAAGGCAAAATTGGAGCAGATATCAATGAATGATTAATAACTACAAAAAGTATGGAGCAAAATTAGGAGAGGGACAGAAAGGAAAttta                                                                                      <  2:240985/91‑2 (MQ=255)
     tAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAACTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAAttttt                                                                                     <  2:3245791/91‑1 (MQ=255)
       aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                                   >  2:3032506/1‑91 (MQ=255)
       aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                                   >  1:2342358/1‑91 (MQ=255)
       aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                                   <  2:1155863/91‑1 (MQ=255)
         ggCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaa                                                                                 <  1:3911078/91‑1 (MQ=255)
          gCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaa                                                                                <  1:895385/91‑1 (MQ=255)
          gCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATaaa                                                                                <  2:3020681/91‑1 (MQ=255)
            aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                              <  2:55494/91‑1 (MQ=255)
             caaTTGGAGCAGTTATCCATGAATGATTAATAACAATAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                             <  2:3914379/90‑1 (MQ=255)
             aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                             <  1:3468789/91‑1 (MQ=255)
             aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                             >  2:1362291/1‑91 (MQ=255)
             aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                             <  1:20045/91‑1 (MQ=255)
              aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                            >  2:330971/1‑91 (MQ=255)
                ttGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAt                                                                          >  1:2140966/1‑91 (MQ=255)
                 tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                         <  2:296291/91‑1 (MQ=255)
                 tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                         >  1:3372830/1‑91 (MQ=255)
                 tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                         >  1:3276541/1‑91 (MQ=255)
                   gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                       >  2:2150906/1‑91 (MQ=255)
                   gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                       >  2:2262701/1‑91 (MQ=255)
                    aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                                      <  1:3631791/91‑1 (MQ=255)
                    aGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCt                                                                      <  2:2571981/91‑1 (MQ=255)
                     gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTCGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                     >  2:3490973/1‑91 (MQ=255)
                     gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                     >  1:1642694/1‑91 (MQ=255)
                     gCAGTTATCCATGAATGATTAAAAAAAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTg                                                                     <  2:2534497/91‑1 (MQ=12)
                      cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgt                                                                    <  1:3587273/91‑2 (MQ=17)
                        tttATCCATGAGTTAGTAATAACAACAAAAAGTATTGGGCAAGGTTAGGAGTGTGACAGAAATGAATTTTATAAAAATTTCATCTCTGTca                                                                  <  2:1631853/90‑1 (MQ=9)
                        gTTGTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAGTTCTGCca                                                                  <  2:3694322/91‑1 (MQ=255)
                        gTTATCCATGGATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                  <  2:2633975/91‑4 (MQ=17)
                        gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                  >  2:2378494/1‑91 (MQ=255)
                        gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                                  <  1:1520340/91‑1 (MQ=255)
                         ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcac                                                                 <  1:1076918/91‑1 (MQ=17)
                         ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcac                                                                 <  1:1852230/91‑1 (MQ=14)
                         ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                                 <  1:3296049/91‑1 (MQ=255)
                         ttATCCATGAATGACTAATAACAGAAAAAAGTATGGGGCAAGAGTAGGAGTGTGTCAGAGATGAAATTTTTATAAAATTCATTGCTGCcac                                                                 <  2:3675936/91‑1 (MQ=255)
                          tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                                <  2:3738775/91‑1 (MQ=17)
                          tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                <  1:2479858/91‑1 (MQ=255)
                          tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcaca                                                                <  2:3133026/91‑1 (MQ=255)
                           agccATGAATGATTAATAACAACAAAAAGAATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCacac                                                               <  2:3628458/89‑1 (MQ=9)
                           aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                               <  1:844407/91‑1 (MQ=255)
                           aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                               >  1:641845/1‑91 (MQ=255)
                            tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                              <  1:2909432/91‑1 (MQ=255)
                            tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                              >  2:3794286/1‑91 (MQ=255)
                            tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACt                                                              <  2:1261130/91‑1 (MQ=255)
                             ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                             <  2:2867210/91‑1 (MQ=255)
                             ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTc                                                             <  1:3533149/91‑1 (MQ=255)
                              cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                            <  1:2348513/91‑1 (MQ=255)
                              cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                            <  1:3648154/91‑1 (MQ=255)
                               atgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCt                                                           >  2:234570/1‑91 (MQ=255)
                                tgtatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGGGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtt                                                          <  2:484326/88‑1 (MQ=255)
                                tgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtt                                                          <  2:23725/91‑1 (MQ=255)
                                 gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGAGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                         <  1:3135355/91‑1 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                        <  2:1006792/91‑1 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                        <  1:2535084/91‑1 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                        <  2:327287/91‑1 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                        <  2:1725175/91‑1 (MQ=255)
                                   atgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttttt                                                       <  1:2518447/91‑1 (MQ=255)
                                    tgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttg                                                      >  1:2218893/1‑91 (MQ=255)
                                      aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                    <  1:2367196/91‑1 (MQ=255)
                                        taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                  >  1:2721731/1‑91 (MQ=255)
                                        taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattg                                                  >  1:2721705/1‑91 (MQ=255)
                                         aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                                 >  1:54090/1‑91 (MQ=255)
                                          ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                                >  1:722449/1‑91 (MQ=255)
                                             acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                             >  2:2486191/1‑91 (MQ=255)
                                             acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                             >  2:1725319/1‑91 (MQ=255)
                                               aacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                           <  2:1972726/91‑1 (MQ=255)
                                                acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                          >  1:1256131/1‑91 (MQ=255)
                                                acaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCa                                          >  1:2041758/1‑91 (MQ=255)
                                                 caaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAt                                         >  2:1421366/1‑91 (MQ=255)
                                                  aaaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtt                                        >  1:421088/1‑91 (MQ=255)
                                                   aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                                       >  2:2217555/1‑91 (MQ=255)
                                                    aaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAtttt                                      >  2:1050433/1‑91 (MQ=255)
                                                     aaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTa                                     <  1:2450082/91‑1 (MQ=255)
                                                      aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                    >  1:2764868/1‑91 (MQ=255)
                                                          tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                >  1:3212691/1‑91 (MQ=255)
                                                           ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGTTTGAATTAGCATTTTACGATCa                               >  2:1868317/1‑91 (MQ=255)
                                                           ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                               <  2:2107305/91‑1 (MQ=255)
                                                           ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                               <  1:2711333/91‑1 (MQ=255)
                                                            gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                              >  2:742524/1‑91 (MQ=255)
                                                             ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                             >  2:1692751/1‑91 (MQ=255)
                                                             ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                             <  1:3762164/91‑1 (MQ=255)
                                                              gCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                            >  2:663760/1‑91 (MQ=255)
                                                              gCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                            >  2:2014918/1‑91 (MQ=255)
                                                              gCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaa                            <  1:1577788/91‑1 (MQ=255)
                                                               cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                           >  1:1328888/1‑91 (MQ=255)
                                                                aaGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAAc                          >  1:3096859/1‑91 (MQ=255)
                                                                   aTTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGt                       >  1:1112409/1‑91 (MQ=255)
                                                                          gtgtGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATAGCATTTTACGATCATAAACAGTCATTATaa               >  1:1981509/1‑91 (MQ=2)
                                                                            gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt              <  1:2589346/91‑1 (MQ=255)
                                                                                 agaaaTGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATCAACAGTCATTATAATTGAGt         >  1:3793370/5‑91 (MQ=18)
                                                                                 agaaATGAATTTTATAAAAATTTCATATCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt         >  2:3655346/1‑91 (MQ=21)
                                                                                   agaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAt       >  1:1922027/1‑91 (MQ=255)
                                                                                   agaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAt       >  2:1914428/1‑91 (MQ=255)
                                                                                    gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt      >  1:612239/1‑91 (MQ=255)
                                                                                      tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg    <  1:179251/91‑1 (MQ=255)
                                                                                      tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg    >  1:1567629/1‑91 (MQ=255)
                                                                                      tGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg    <  1:671429/91‑1 (MQ=255)
                                                                                       gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa   >  1:1283089/1‑91 (MQ=255)
                                                                                        aaaTTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa  <  1:1664147/91‑1 (MQ=255)
                                                                                          |                                                                                        
ATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAA  >  USA300TCH1516_ALE/770776‑770954

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: