Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,866 A→T 8.8% 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 8.8% 382.0 / 22.3 113intergenic (+52/‑198)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base A (58/45);  new base T (3/7);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 114 total aligned reads displayed.

ATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGA  >  USA300TCH1516_ALE/770776‑770953
                                                                                          |                                                                                       
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)
      aaaGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTta                                                                                   >  1:1053133/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:2874325/91‑1 (MQ=255)
            aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                             >  1:2722207/1‑91 (MQ=255)
            aaaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAt                                                                             >  1:2721877/1‑91 (MQ=255)
             aaaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAAtt                                                                            >  1:2059504/1‑91 (MQ=255)
              aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                           >  1:550457/1‑91 (MQ=255)
              aaTTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTc                                                                           <  1:3084588/91‑1 (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)
                       aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtc                                                                  <  1:2886246/91‑3 (MQ=17)
                        gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTgtca                                                                 <  1:3099755/91‑4 (MQ=17)
                        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:1133654/1‑91 (MQ=255)
                           aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCacac                                                              >  1:2232223/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)
                                 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)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCtaat                                                       <  1:2582593/91‑4 (MQ=2)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                       >  1:243119/1‑91 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                       >  1:1852017/1‑91 (MQ=255)
                                  aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCtttt                                                       <  1:3189800/91‑1 (MQ=255)
                                     gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                    <  1:1640639/91‑1 (MQ=255)
                                     gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttga                                                    >  1:3818351/1‑91 (MQ=255)
                                      aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                   >  1:2520084/1‑91 (MQ=255)
                                      aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                   >  1:3571537/1‑91 (MQ=255)
                                      aTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgat                                                   <  1:1399048/91‑1 (MQ=255)
                                         aataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTttgattga                                                >  1:1018175/1‑91 (MQ=255)
                                          ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                               >  1:3935935/1‑91 (MQ=255)
                                          ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                               >  1:1970447/1‑91 (MQ=255)
                                          ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGaa                                               >  1:1972601/1‑91 (MQ=255)
                                           taaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAt                                              <  1:3411534/91‑1 (MQ=255)
                                            aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                             <  1:337186/91‑1 (MQ=255)
                                            aacaacAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                             <  1:336878/91‑1 (MQ=255)
                                             acaacaAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTa                                            >  1:497326/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:3160809/91‑1 (MQ=255)
                                                      aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                   <  1:2245523/91‑1 (MQ=255)
                                                       gTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACg                                  >  1:1273361/1‑91 (MQ=255)
                                                         aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                                >  1:2104325/1‑91 (MQ=255)
                                                          tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                               <  1:2182310/91‑1 (MQ=255)
                                                          tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                               >  1:1571222/1‑91 (MQ=255)
                                                           ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                              >  1:510250/1‑91 (MQ=255)
                                                           ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                              >  1:510064/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:2663868/1‑91 (MQ=255)
                                                                      aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                   <  1:1999169/91‑1 (MQ=255)
                                                                      aGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                   <  1:3772122/91‑1 (MQ=255)
                                                                       ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                  <  1:911528/91‑1 (MQ=255)
                                                                       ggAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAtt                  >  1:575465/1‑91 (MQ=255)
                                                                        gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                 <  1:1178726/91‑1 (MQ=255)
                                                                        gAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                 >  1:322715/1‑91 (MQ=255)
                                                                         aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                >  1:1802078/1‑91 (MQ=14)
                                                                         aGTGTGGCAGAAATGAATTTTATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                >  1:83481/1‑91 (MQ=14)
                                                                          gtgtGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtata               >  1:843355/1‑91 (MQ=255)
                                                                            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:1381638/91‑1 (MQ=255)
                                                                              gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg           <  1:2960081/91‑1 (MQ=255)
                                                                              gACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg           <  1:1381512/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)
                                                                                 agagATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt        >  1:3580204/1‑91 (MQ=255)
                                                                                    gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt     >  1:2607177/1‑91 (MQ=255)
                                                                                    gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt     >  1:2184089/1‑91 (MQ=255)
                                                                                    gaTGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt     >  1:140868/1‑91 (MQ=255)
                                                                                       gAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGa  >  1:2156992/1‑91 (MQ=255)
                                                                                          |                                                                                       
ATAACTAAAGGCAAAATTGGAGCAGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAATTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGA  >  USA300TCH1516_ALE/770776‑770953

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