Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,770 A→T 7.3% S954S (TCT→TCA 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,7700AT 7.3% 547.2 / 23.6 164S954S (TCT→TCAnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base A (79/73);  new base T (5/7);  total (84/80)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 165 total aligned reads displayed.

GCGCGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTC  >  USA300TCH1516_ALE/2517681‑2517860
                                                                                         |                                                                                          
   cgTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATg                                                                                        <  1:2550565/91‑1 (MQ=255)
   cgTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATg                                                                                        <  2:2537448/91‑1 (MQ=255)
    gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                       <  1:2383326/91‑1 (MQ=255)
    gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                       <  1:23101/91‑1 (MQ=255)
    gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                       <  1:804031/91‑1 (MQ=255)
    gTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGa                                                                                       >  2:745525/1‑91 (MQ=255)
     ttCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGaa                                                                                      >  1:738110/1‑91 (MQ=255)
     ttCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGaa                                                                                      >  1:3352646/1‑91 (MQ=255)
       cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                    <  1:1357895/91‑1 (MQ=255)
       cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                    <  1:2390138/91‑1 (MQ=255)
        gCTGGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGGAATa                                                                                   <  1:3059896/91‑1 (MQ=255)
        gCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATa                                                                                   <  1:1306216/91‑1 (MQ=255)
        gCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATa                                                                                   <  1:735627/91‑1 (MQ=255)
         cTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATa                                                                                   >  1:548817/1‑90 (MQ=255)
          tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                 <  2:2535186/91‑1 (MQ=255)
          tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                 <  1:2893600/91‑1 (MQ=255)
          tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCAGATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                 <  2:2535436/91‑1 (MQ=255)
            gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATCTGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                               >  1:3072563/1‑91 (MQ=255)
            gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                               >  1:2238300/1‑91 (MQ=255)
             aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                              <  2:1909312/91‑1 (MQ=255)
              aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGt                                                                             <  1:3032636/91‑1 (MQ=255)
                 atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                          >  2:3592236/1‑91 (MQ=255)
                 atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                          >  1:2086572/1‑91 (MQ=255)
                 atCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAtt                                                                          >  2:3150191/1‑91 (MQ=255)
                  tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTATCATTCGTAGCAGATGAAATACTTAGTATTg                                                                         >  1:1775618/1‑91 (MQ=255)
                  tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTg                                                                         >  1:732229/1‑91 (MQ=255)
                   cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                        <  2:831438/91‑1 (MQ=255)
                   cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                        <  1:403272/91‑1 (MQ=255)
                        aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctg                                                                   <  1:403613/91‑3 (MQ=37)
                        aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                   >  1:291531/1‑91 (MQ=255)
                        aaCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAg                                                                   <  1:3867079/91‑1 (MQ=255)
                         aCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctgc                                                                  <  1:715999/91‑4 (MQ=25)
                           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)
                             cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTATCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  2:3414611/1‑91 (MQ=255)
                             cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  2:2494173/1‑91 (MQ=255)
                             cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              <  2:3648733/91‑1 (MQ=255)
                             cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  2:3361103/1‑91 (MQ=255)
                             cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  2:2179606/1‑91 (MQ=255)
                              aTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAt                                                             >  1:1700113/1‑91 (MQ=255)
                               ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                            >  1:3423805/1‑91 (MQ=255)
                                tCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTa                                                           >  1:646060/1‑91 (MQ=255)
                                 ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                          <  2:1905257/91‑1 (MQ=255)
                                  cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                         >  2:645123/1‑91 (MQ=255)
                                   aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                        >  2:1481574/1‑91 (MQ=255)
                                     tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                      <  1:2349036/91‑1 (MQ=255)
                                     tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                      <  2:2163068/91‑1 (MQ=255)
                                     tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                      <  2:936215/91‑1 (MQ=255)
                                      ttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCt                                                     >  1:1387906/1‑91 (MQ=255)
                                       tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                    >  2:1566952/1‑91 (MQ=255)
                                        gTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGt                                                   >  2:130497/1‑91 (MQ=255)
                                        gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                  >  2:2680916/1‑76 (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)
                                           cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                >  2:2169830/1‑91 (MQ=255)
                                            ttcGAGATCTTCATATGATTTTTGTTATTTTTTATCATTTGTAGCAGATTAAATACTTAGTATTTCTTCATCTACATTACGTGCTGTATCa                                               >  2:3620431/1‑91 (MQ=255)
                                            ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                               <  1:102840/91‑1 (MQ=255)
                                                 gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                          >  2:2976636/1‑91 (MQ=255)
                                                 gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                          <  2:2008346/91‑1 (MQ=255)
                                                  aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                         >  2:76176/1‑91 (MQ=255)
                                                    cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                               <  1:2297688/67‑1 (MQ=255)
                                                    cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                               >  2:2297688/1‑67 (MQ=255)
                                                      tCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGAc                                     <  1:917378/91‑1 (MQ=255)
                                                       cATATAATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                    <  2:173125/91‑1 (MQ=255)
                                                        atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                   <  1:1661531/91‑1 (MQ=255)
                                                        atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                   >  1:187602/1‑91 (MQ=255)
                                                           tGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTc                                <  1:190619/91‑1 (MQ=255)
                                                              tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                             >  1:772216/1‑91 (MQ=255)
                                                              tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                             <  1:2686428/91‑1 (MQ=255)
                                                                   gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                        <  1:1746951/91‑1 (MQ=255)
                                                                   gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                        <  1:982790/91‑1 (MQ=255)
                                                                   gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                        <  1:982443/91‑1 (MQ=255)
                                                                   gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                        >  1:779216/1‑91 (MQ=255)
                                                                    tgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcaca                       >  2:2787411/1‑91 (MQ=255)
                                                                     gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                      >  2:1527297/1‑91 (MQ=255)
                                                                      aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                     <  2:299207/91‑1 (MQ=255)
                                                                       tAATTTACCATTCGTAGCAGATGAAATACTTCGTATTGCATGAGCTACATTACGTGCTGTATCAATAGATGGACGATTCACTCTCACAGaa                    >  2:2953644/1‑91 (MQ=255)
                                                                       tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAACTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                    <  1:3882591/91‑1 (MQ=255)
                                                                        aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                   <  1:898643/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)
                                                                            tACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcatc               <  2:878240/91‑1 (MQ=255)
                                                                              ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat             >  2:3286130/1‑91 (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:502073/91‑1 (MQ=255)
                                                                                      aGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa     >  2:312362/1‑91 (MQ=25)
                                                                                      aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa     <  1:1399703/91‑1 (MQ=255)
                                                                                       gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAc    <  1:1485582/91‑1 (MQ=255)
                                                                                        cAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt   >  2:1074957/1‑91 (MQ=255)
                                                                                         aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc  <  1:2502179/91‑1 (MQ=255)
                                                                                         |                                                                                          
GCGCGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTC  >  USA300TCH1516_ALE/2517681‑2517860

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: