Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,775 A→C 11.1% S953A (TCA→GCA)  narG ← Respiratory nitrate reductase 1 alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,517,7750AC11.1% 504.4 / 50.4 163S953A (TCA→GCA) narGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base A (77/67);  new base C (9/9);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

Only 100 of 170 total aligned reads displayed.

GTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACG  >  USA300TCH1516_ALE/2517685‑2517864
                                                                                          |                                                                                         
gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                           <  1:2383326/91‑1 (MQ=255)
gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                           <  1:23101/91‑1 (MQ=255)
gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                           >  2:745525/1‑91 (MQ=255)
 ttCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGaa                                                                                          >  1:737910/1‑91 (MQ=255)
   cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                        <  1:1357895/91‑1 (MQ=255)
   cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                        <  1:2390138/91‑1 (MQ=255)
    gCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATa                                                                                       <  1:735627/91‑1 (MQ=255)
     cTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATa                                                                                       <  2:548817/90‑1 (MQ=255)
      tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                     <  2:2535186/91‑1 (MQ=255)
      tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                     <  1:2893600/91‑1 (MQ=255)
      tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                     <  2:2535787/91‑1 (MQ=255)
       aGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACtt                                                                                    <  1:3259962/91‑1 (MQ=255)
        gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTTTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                   >  2:1073550/1‑91 (MQ=255)
        gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                   <  1:787550/91‑1 (MQ=255)
        gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                   >  1:2238300/1‑91 (MQ=255)
          aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGt                                                                                 <  1:3032636/91‑1 (MQ=255)
            tatCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAt                                                                               <  2:3160228/91‑1 (MQ=255)
             atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                              >  2:3592236/1‑91 (MQ=255)
             atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                              >  1:2086572/1‑91 (MQ=255)
               cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                            <  2:1930537/91‑3 (MQ=38)
               cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                            <  1:403337/91‑1 (MQ=255)
               cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                            <  1:403272/91‑1 (MQ=255)
               cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                            <  2:831438/91‑1 (MQ=255)
               cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                            <  1:363764/91‑1 (MQ=255)
                ttttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCa                                                                           >  1:1323241/1‑91 (MQ=255)
                 tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                          >  1:3069822/1‑91 (MQ=255)
                 tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                          <  1:727170/91‑1 (MQ=255)
                    aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                       >  1:291443/1‑91 (MQ=255)
                     aCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                      >  2:147246/1‑91 (MQ=255)
                       tgCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                    >  2:1031553/2‑91 (MQ=255)
                       ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                    <  2:2593129/91‑1 (MQ=255)
                       ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                    >  1:3295158/1‑91 (MQ=255)
                       ggCATTCCAGTTTGTTCTTCGAGATCATCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                    >  1:1615245/1‑91 (MQ=255)
                        gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCTAc                                                                   <  1:906258/91‑1 (MQ=25)
                        gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                   <  1:1881963/91‑1 (MQ=255)
                         cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTATCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                  >  2:3414611/1‑91 (MQ=255)
                         cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                  >  2:2494173/1‑91 (MQ=255)
                         cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                  >  2:2179606/1‑91 (MQ=255)
                         cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                  >  2:3361103/1‑91 (MQ=255)
                          aTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAt                                                                 >  1:1700113/1‑91 (MQ=255)
                           ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                >  1:3423790/1‑91 (MQ=255)
                            tCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTa                                                               >  1:646060/1‑91 (MQ=255)
                             ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                              >  1:925244/1‑91 (MQ=255)
                             ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                              <  2:1905257/91‑1 (MQ=255)
                              cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                             >  1:594561/1‑91 (MQ=255)
                                 tttGTTCTTCGAGATCTTCATATGATTTTTGTTATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                          <  1:2692310/91‑1 (MQ=255)
                                 tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                          <  1:2349036/91‑1 (MQ=255)
                                 tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                          <  2:2163068/91‑1 (MQ=255)
                                 tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                          <  2:936215/91‑1 (MQ=255)
                                   tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                        >  2:1566952/1‑91 (MQ=255)
                                    gTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                       >  2:130497/1‑91 (MQ=255)
                                    gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                      >  2:2680916/1‑76 (MQ=255)
                                    gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                       >  1:878240/1‑91 (MQ=255)
                                    gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                       <  2:2310775/91‑1 (MQ=255)
                                      tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                     <  2:846315/91‑1 (MQ=255)
                                       cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                              >  2:2763052/1‑81 (MQ=255)
                                       cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                              <  1:2763052/81‑1 (MQ=255)
                                       cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                    <  2:389699/91‑1 (MQ=255)
                                       cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                    >  1:1260937/1‑91 (MQ=255)
                                        ttcGAGATCTTCATATGATTTTTGTTATTTTTTATCATTTGTAGCAGATTAAATACTTAGTATTTCTTCATCTACATTACGTGCTGTATCa                                                   >  2:3620431/1‑91 (MQ=255)
                                        ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                                   <  1:102840/91‑1 (MQ=255)
                                          cGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAAt                                                 >  2:2179167/1‑91 (MQ=255)
                                             gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                              >  2:2976636/1‑91 (MQ=255)
                                             gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                              <  2:2008346/91‑1 (MQ=255)
                                                  tCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGAc                                         <  1:917378/91‑1 (MQ=255)
                                                   cATATGATTTTTTTTTTTTTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAAGTGCTGTATCAATACGTGGACg                                        >  2:3121937/1‑91 (MQ=255)
                                                   cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                        <  1:745525/91‑1 (MQ=255)
                                                   cATATAATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                        <  2:173125/91‑1 (MQ=255)
                                                    atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                       >  1:187602/1‑91 (MQ=255)
                                                     tatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                      >  1:1699592/1‑91 (MQ=255)
                                                          tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                                 <  1:2686428/91‑1 (MQ=255)
                                                           ttttGTGATAATTTACCATTCGTAGCAAATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctc                                >  1:2284681/1‑91 (MQ=255)
                                                            tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                               >  2:756990/1‑91 (MQ=255)
                                                               gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                            <  1:1746951/91‑1 (MQ=255)
                                                               gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                            >  1:779216/1‑91 (MQ=255)
                                                               gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                            <  1:982790/91‑1 (MQ=255)
                                                                tgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcaca                           >  2:2787411/1‑91 (MQ=255)
                                                                 gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                          >  2:1527297/1‑91 (MQ=255)
                                                                   tAATTTACCATTCGTAGCAGATGAAATACTTCGTATTGCATGAGCTACATTACGTGCTGTATCAATAGATGGACGATTCACTCTCACAGaa                        >  2:2953644/1‑91 (MQ=255)
                                                                   tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAACTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                        <  1:3882591/91‑1 (MQ=255)
                                                                     aTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatc                      >  2:91245/1‑91 (MQ=255)
                                                                     aTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatc                      <  2:1736957/91‑1 (MQ=255)
                                                                       ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                    <  2:2629801/91‑1 (MQ=255)
                                                                          ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                 >  2:780935/1‑91 (MQ=255)
                                                                             ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                               <  1:1016991/74‑1 (MQ=255)
                                                                             ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg              <  2:653932/91‑1 (MQ=255)
                                                                              tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt             >  2:1095903/1‑91 (MQ=255)
                                                                               tgTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa            >  2:2239666/2‑91 (MQ=25)
                                                                               cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa            <  1:503071/91‑1 (MQ=255)
                                                                                  agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa         >  2:3402483/5‑91 (MQ=25)
                                                                                  agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa         >  2:3554662/5‑91 (MQ=25)
                                                                                  aGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa         >  2:312362/1‑91 (MQ=25)
                                                                                  aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa         <  1:1399703/91‑1 (MQ=255)
                                                                                   gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAc        >  1:141054/1‑91 (MQ=255)
                                                                                     aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc      <  1:2502179/91‑1 (MQ=255)
                                                                                      gATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc     <  1:1488426/91‑1 (MQ=255)
                                                                                      gATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc     >  1:472795/1‑91 (MQ=255)
                                                                                         gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg  >  2:3743122/1‑91 (MQ=255)
                                                                                         gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg  >  1:3855087/1‑91 (MQ=255)
                                                                                         gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg  >  1:1384204/1‑91 (MQ=255)
                                                                                          |                                                                                         
GTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACG  >  USA300TCH1516_ALE/2517685‑2517864

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