Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,860 G→A 8.6% intergenic (+46/‑204) 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,8600GA 8.6% 376.5 / 17.9 116intergenic (+46/‑204)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base G (58/48);  new base A (3/7);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 122 total aligned reads displayed.

ATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATT  >  USA300TCH1516_ALE/770770‑770950
                                                                                          |                                                                                          
aTTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGTAGTGTGACagag                                                                                            <  1:3181258/91‑1 (MQ=255)
 ttACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAgaga                                                                                           <  1:3791931/91‑1 (MQ=255)
  tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                          <  1:1800984/91‑1 (MQ=255)
  tACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGAt                                                                                          <  1:174729/91‑1 (MQ=255)
   aCAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATg                                                                                         >  1:3952550/1‑91 (MQ=255)
   aCAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATg                                                                                         >  1:82282/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)
             aaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTtat                                                                               >  1:2768853/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:2722207/1‑91 (MQ=255)
                  aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                          <  1:2874325/91‑1 (MQ=255)
                  aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                          <  1:2874804/91‑1 (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:1010887/91‑1 (MQ=255)
                    aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                        >  1:550457/1‑91 (MQ=255)
                     aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                       <  1:1926359/91‑1 (MQ=255)
                     aTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCa                                                                       <  1:3431038/91‑1 (MQ=255)
                      ttGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAt                                                                      <  1:485071/91‑1 (MQ=255)
                       tGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCAtt                                                                     >  1:746115/1‑91 (MQ=255)
                         gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                   >  1:930279/1‑91 (MQ=255)
                         gAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTc                                                                   >  1:121657/1‑91 (MQ=255)
                           gCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTg                                                                 >  1:1993397/1‑91 (MQ=255)
                            cAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGc                                                                >  1:2112570/1‑91 (MQ=255)
                              gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                              <  1:3099896/91‑4 (MQ=17)
                              gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCca                                                              <  1:934096/91‑1 (MQ=255)
                               ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCcac                                                             <  1:378467/91‑1 (MQ=255)
                                tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTcaca                                                            <  1:1753057/91‑1 (MQ=14)
                                 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)
                                    cATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTcc                                                        <  1:4104396/91‑1 (MQ=255)
                                      tgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtt                                                      <  1:2934184/91‑1 (MQ=255)
                                       gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                     <  1:2326469/91‑1 (MQ=255)
                                       gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCttt                                                     <  1:4027861/91‑1 (MQ=255)
                                        aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCtaat                                                    <  1:2582593/91‑4 (MQ=2)
                                        aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                    <  1:3189800/91‑1 (MQ=255)
                                        aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                    >  1:243119/1‑91 (MQ=255)
                                        aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                    >  1:1852017/1‑91 (MQ=255)
                                           gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                 <  1:1640639/91‑1 (MQ=255)
                                           gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                 >  1:3818351/1‑91 (MQ=255)
                                            aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                <  1:1399048/91‑1 (MQ=255)
                                            aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                >  1:3571537/1‑91 (MQ=255)
                                               aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                             >  1:1018175/1‑91 (MQ=255)
                                                ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                            >  1:1972601/1‑91 (MQ=255)
                                                ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                            >  1:3935935/1‑91 (MQ=255)
                                                ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                            >  1:1970447/1‑91 (MQ=255)
                                                 taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAt                                           <  1:3411534/91‑1 (MQ=255)
                                                  aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                          <  1:3112305/91‑1 (MQ=255)
                                                  aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                          <  1:337186/91‑1 (MQ=255)
                                                  aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                          <  1:336878/91‑1 (MQ=255)
                                                   acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                         >  1:496747/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)
                                                            aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                <  1:2246148/91‑1 (MQ=255)
                                                            aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                <  1:2245523/91‑1 (MQ=255)
                                                            aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                <  1:3160809/91‑1 (MQ=255)
                                                               aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                             >  1:2104325/1‑91 (MQ=255)
                                                                tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                            >  1:1571222/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)
                                                                 ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                           >  1:2336528/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:1999169/91‑1 (MQ=255)
                                                                            aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                <  1:3772122/91‑1 (MQ=255)
                                                                            aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                >  1:2663868/1‑91 (MQ=255)
                                                                             ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt               >  1:1139591/1‑91 (MQ=255)
                                                                             ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt               <  1:911528/91‑1 (MQ=255)
                                                                              gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta              >  1:322715/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:4020385/91‑1 (MQ=255)
                                                                                  gtgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt          <  1:4020387/91‑1 (MQ=255)
                                                                                   tgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt         <  1:3200977/91‑1 (MQ=255)
                                                                                   tgACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAtt         <  1:3198840/91‑1 (MQ=255)
                                                                                    gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg        <  1:2960081/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)
                                                                                      cAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg      >  1:10473/1‑91 (MQ=255)
                                                                                      cAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg      >  1:3009693/1‑91 (MQ=18)
                                                                                          gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt  >  1:140868/1‑91 (MQ=255)
                                                                                          gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt  >  1:2607177/1‑91 (MQ=255)
                                                                                          gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt  >  1:2184089/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTACAATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATT  >  USA300TCH1516_ALE/770770‑770950

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: