Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,770 A→T 6.9% 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 6.9% 564.4 / 22.4 169S954S (TCT→TCAnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base A (63/94);  new base T (2/10);  total (65/104)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

AGCGCGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTC  >  USA300TCH1516_ALE/2517680‑2517860
                                                                                          |                                                                                          
agcGCGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCa                                                                                            >  2:1570199/1‑91 (MQ=255)
 gcgcGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAg                                                                                           >  2:2784982/1‑91 (MQ=255)
      ttCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGaa                                                                                      <  1:2179406/91‑1 (MQ=255)
      ttCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGaa                                                                                      >  1:201698/1‑91 (MQ=255)
        cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                    >  2:225642/1‑91 (MQ=255)
        cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                    >  1:2012488/1‑91 (MQ=255)
           tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                 <  1:2538069/91‑1 (MQ=255)
             gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                               >  2:440056/1‑91 (MQ=255)
              aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACtaag                                                                              <  1:232753/91‑4 (MQ=255)
              aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                              <  1:1847724/91‑1 (MQ=255)
              aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                              <  2:3645257/91‑1 (MQ=255)
              aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                              <  1:3396829/91‑1 (MQ=255)
               aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGt                                                                             >  2:726169/1‑91 (MQ=255)
               aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGt                                                                             <  2:3221922/91‑1 (MQ=255)
                atatCTTTTAACGGCTTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTCCCATTCGTAGCAGATGAAATACTTAGTa                                                                            <  2:827964/91‑1 (MQ=255)
                   tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAttt                                                                         <  1:2892364/91‑2 (MQ=38)
                   tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                         <  1:3709752/91‑2 (MQ=255)
                   tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                         <  2:3046474/91‑2 (MQ=255)
                    cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                        <  1:2342482/91‑3 (MQ=38)
                    cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                        <  2:1006879/91‑1 (MQ=255)
                      tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                      <  2:2714821/91‑1 (MQ=255)
                      tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                      <  1:1338076/91‑1 (MQ=255)
                       ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                     <  1:1165544/91‑1 (MQ=38)
                       ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                     <  1:3257492/91‑1 (MQ=38)
                        tAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATct                                                                    <  1:2331376/91‑2 (MQ=37)
                          aCGGCATTCCAGTTTGTTCTTCGAGATCTTCGTATGATTTTTGTGATAATTTAGCATGCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                  <  1:3157735/91‑1 (MQ=255)
                           cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCt                                                                 <  2:2603314/91‑1 (MQ=25)
                           cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCt                                                                 <  2:2418777/91‑1 (MQ=255)
                            ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                <  2:2677416/91‑1 (MQ=255)
                              cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  1:2157444/1‑91 (MQ=255)
                              cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  1:1672946/1‑91 (MQ=255)
                              cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                              >  2:1078403/1‑91 (MQ=255)
                                ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                            <  1:104777/91‑1 (MQ=255)
                                ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                            <  1:3582862/91‑1 (MQ=255)
                                ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                            >  2:3487059/1‑91 (MQ=255)
                                  ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                          <  2:2542254/91‑1 (MQ=255)
                                   cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                         <  1:3435960/91‑1 (MQ=255)
                                    aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                        >  1:443100/1‑91 (MQ=255)
                                      tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                      >  2:2317227/1‑91 (MQ=255)
                                        tGTTCTTCGGGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                    <  2:1820638/91‑1 (MQ=255)
                                        tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                    >  1:1272400/1‑91 (MQ=255)
                                        tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                    >  1:1745315/1‑91 (MQ=255)
                                          ttcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTa                                                  <  1:1191144/91‑1 (MQ=255)
                                           tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                 <  2:3221637/91‑1 (MQ=255)
                                            cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                <  1:3845333/91‑1 (MQ=255)
                                             ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                               >  2:1742873/1‑91 (MQ=255)
                                              tcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCaa                                              <  1:1364742/91‑1 (MQ=255)
                                                  gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                          <  1:3089714/91‑1 (MQ=255)
                                                  gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTACATCAGCTACATTACGTGCTGTATCAATACg                                          >  1:1743721/1‑91 (MQ=255)
                                                   aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                         <  1:1651330/91‑1 (MQ=255)
                                                   aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                         >  2:3243496/1‑91 (MQ=255)
                                                     cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                       >  1:3542505/1‑91 (MQ=255)
                                                     cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                       >  2:143906/1‑91 (MQ=255)
                                                       tCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGAc                                     <  1:896164/91‑1 (MQ=255)
                                                        cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                    <  2:3388184/91‑1 (MQ=255)
                                                         atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                   <  2:2550137/91‑1 (MQ=255)
                                                          tatGATTTTTTTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                  >  2:1273075/1‑91 (MQ=255)
                                                          tatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCGGCTACATTACGTGCTGTATCAATACGTGGACGAt                                  <  2:3619143/91‑1 (MQ=255)
                                                          tatGATTATTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                  <  1:1953284/91‑1 (MQ=255)
                                                           atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                 <  2:2230261/91‑1 (MQ=255)
                                                           atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                 >  2:2088317/1‑91 (MQ=255)
                                                               tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                             <  1:3523495/91‑1 (MQ=255)
                                                               tttttGTGATAATTTACCAGTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGcc                             >  2:3451747/1‑90 (MQ=255)
                                                                ttttGTGATAATTTACCATTCGTAGCAGATGAAATACTGAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctc                            <  2:2957818/91‑1 (MQ=255)
                                                                 tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTTCATTACGTGCTGTATCAATACGTGGACGATTCGctct                           <  1:3340429/91‑1 (MQ=255)
                                                                 tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                           >  2:3537280/1‑91 (MQ=255)
                                                                 tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                           <  1:666547/91‑1 (MQ=255)
                                                                   tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                         <  2:2926893/91‑1 (MQ=255)
                                                                   tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                         >  1:3638732/1‑91 (MQ=255)
                                                                    gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                        >  1:634149/1‑91 (MQ=255)
                                                                      gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                      >  1:2617120/1‑91 (MQ=255)
                                                                      gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                      <  2:3894786/91‑1 (MQ=255)
                                                                       ttAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCTCAGa                     >  2:3127895/2‑91 (MQ=255)
                                                                       aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                     <  1:2780704/91‑1 (MQ=255)
                                                                        tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                    >  2:3026221/1‑91 (MQ=255)
                                                                        tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                    <  1:1199208/91‑1 (MQ=255)
                                                                        tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                    <  1:1037870/91‑1 (MQ=255)
                                                                         aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                   <  2:3844989/91‑1 (MQ=255)
                                                                           tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                 >  1:2052680/1‑91 (MQ=255)
                                                                           tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                 <  1:2003992/91‑1 (MQ=255)
                                                                            ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                >  2:1289828/1‑91 (MQ=255)
                                                                            ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                >  1:218425/1‑91 (MQ=255)
                                                                            ttACCATGCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAGTACGTGCTGGATCAATACGTGGGCGATTCGCTCGCACAGAatcat                <  2:688112/91‑1 (MQ=255)
                                                                              aCCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAAtcatca              <  2:2506451/91‑1 (MQ=255)
                                                                                cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGTATCattatt            >  2:1640115/1‑91 (MQ=255)
                                                                                cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCAtt            <  2:1147632/91‑1 (MQ=255)
                                                                                  ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg          >  2:1946329/1‑91 (MQ=255)
                                                                                  ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg          <  2:580473/91‑1 (MQ=255)
                                                                                  ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg          >  1:3820440/1‑91 (MQ=255)
                                                                                  ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg          <  2:201698/91‑1 (MQ=255)
                                                                                   tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt         >  1:242829/1‑91 (MQ=255)
                                                                                   tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt         >  1:189993/1‑91 (MQ=255)
                                                                                   tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt         >  1:243572/1‑91 (MQ=255)
                                                                                    cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa        >  1:68692/1‑91 (MQ=255)
                                                                                     gTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt       >  1:1322373/1‑91 (MQ=25)
                                                                                      tAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATc      >  2:362590/1‑91 (MQ=25)
                                                                                       aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa     >  1:1235425/1‑91 (MQ=255)
                                                                                         cAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt   >  2:624995/1‑91 (MQ=255)
                                                                                          aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAGTc  >  2:1202678/1‑91 (MQ=255)
                                                                                          aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc  <  2:901614/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGCGCGTTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTC  >  USA300TCH1516_ALE/2517680‑2517860

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: