Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,775 A→C 12.0% 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,7750AC12.0% 507.4 / 58.7 174S953A (TCA→GCA) narGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base A (64/89);  new base C (5/16);  total (69/105)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

TTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGT  >  USA300TCH1516_ALE/2517686‑2517865
                                                                                         |                                                                                          
  cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                         <  2:2534016/91‑1 (MQ=255)
  cGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAAt                                                                                         >  1:2012488/1‑91 (MQ=255)
     tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                      <  1:2538069/91‑1 (MQ=255)
     tAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACt                                                                                      <  1:1669081/91‑1 (MQ=255)
       gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                    <  2:3525487/91‑1 (MQ=255)
       gAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTa                                                                                    >  2:440056/1‑91 (MQ=255)
        aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACtaag                                                                                   <  1:232753/91‑4 (MQ=255)
        aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                                   <  2:3645257/91‑1 (MQ=255)
        aaaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAg                                                                                   <  1:3396829/91‑1 (MQ=255)
         aaTATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGt                                                                                  <  2:3394889/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)
           tatCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAt                                                                                <  2:1260840/91‑1 (MQ=255)
           tatCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTAt                                                                                >  2:2864403/1‑91 (MQ=255)
             tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAttt                                                                              <  1:2892364/91‑2 (MQ=38)
             tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                              <  1:3709752/91‑2 (MQ=255)
             tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                              <  1:1960131/91‑2 (MQ=255)
             tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGCAATACTAAGTAttt                                                                              <  2:3046473/91‑2 (MQ=255)
             tCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTg                                                                              <  1:157706/91‑1 (MQ=255)
              cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTAtttc                                                                             <  1:2342482/91‑3 (MQ=38)
              cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                             <  2:1006879/91‑1 (MQ=255)
                tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                           <  1:1338076/91‑1 (MQ=255)
                 ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                          <  1:1283951/91‑1 (MQ=38)
                 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)
                    aCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctgc                                                                       <  2:3142803/91‑4 (MQ=37)
                     cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCt                                                                      <  2:2603314/91‑1 (MQ=25)
                     cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCt                                                                      <  2:2418777/91‑1 (MQ=255)
                       gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                    <  1:1783855/91‑1 (MQ=255)
                        cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                   >  1:2157444/1‑91 (MQ=255)
                        cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                   >  2:1078403/1‑91 (MQ=255)
                          ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                 <  1:3582912/91‑1 (MQ=255)
                          ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                 <  1:104777/91‑1 (MQ=255)
                          ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                 <  1:427555/91‑1 (MQ=255)
                          ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                 >  2:3487059/1‑91 (MQ=255)
                            ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                               >  1:2046530/1‑91 (MQ=255)
                                tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                           >  2:2317227/1‑91 (MQ=255)
                                  tGTTCTTCGGGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                         <  2:1820638/91‑1 (MQ=255)
                                  tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                         >  1:1745315/1‑91 (MQ=255)
                                    ttcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTa                                                       <  1:1191144/91‑1 (MQ=255)
                                     tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                      <  1:3455753/91‑1 (MQ=255)
                                       ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                                    >  2:1742873/1‑91 (MQ=255)
                                        tcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCaa                                                   <  1:1364718/91‑1 (MQ=255)
                                          gagaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATa                                                 <  1:3663749/91‑1 (MQ=255)
                                           agaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATAc                                                <  2:2681522/91‑1 (MQ=255)
                                            gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                               <  1:3089714/91‑1 (MQ=255)
                                            gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                               <  2:1328909/91‑1 (MQ=255)
                                            gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTACATCAGCTACATTACGTGCTGTATCAATACg                                               >  1:1743721/1‑91 (MQ=255)
                                             aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAACTACATTACGTGCTGTATCAATACGt                                              <  2:2200769/91‑1 (MQ=255)
                                               cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                            >  1:3542505/1‑91 (MQ=255)
                                               cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                            <  2:1467379/91‑1 (MQ=255)
                                               cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                            <  2:644491/91‑1 (MQ=255)
                                                 tCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGAc                                          <  1:896164/91‑1 (MQ=255)
                                                  cATATGATTTTTGTGATAATTTACCATTCGTATCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                         <  1:1234099/91‑1 (MQ=255)
                                                  cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                         <  2:3231986/91‑1 (MQ=255)
                                                   atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                        <  2:2550137/91‑1 (MQ=255)
                                                    tatGATTTTTTTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                       >  2:1273075/1‑91 (MQ=255)
                                                     atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                      <  2:2230261/91‑1 (MQ=255)
                                                       gATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCg                                    >  1:460081/1‑91 (MQ=255)
                                                         tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                                  <  1:3523495/91‑1 (MQ=255)
                                                          ttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctc                                 <  2:2957790/91‑1 (MQ=255)
                                                           tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTTCATTACGTGCTGTATCAATACGTGGACGATTCGctct                                <  1:3340429/91‑1 (MQ=255)
                                                           tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                <  1:2071518/91‑1 (MQ=255)
                                                           tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                <  1:666547/91‑1 (MQ=255)
                                                              gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                             >  1:634149/1‑91 (MQ=255)
                                                              gtaaTAATTCACCATTCGTAGCAGATGAAATACTTAGTATTTCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGTTCTcac                             >  2:2623957/4‑91 (MQ=255)
                                                                gATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAg                           >  1:2617120/1‑91 (MQ=255)
                                                                 aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                          >  2:3127790/1‑91 (MQ=255)
                                                                 aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                          <  1:2780704/91‑1 (MQ=255)
                                                                  tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                         <  1:1037870/91‑1 (MQ=255)
                                                                  tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                         >  2:357573/1‑91 (MQ=255)
                                                                  tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                         >  2:3026221/1‑91 (MQ=255)
                                                                   aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                        <  2:3845930/91‑1 (MQ=255)
                                                                   aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                        <  2:3844989/91‑1 (MQ=255)
                                                                     tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                      <  1:2003992/91‑1 (MQ=255)
                                                                      ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                     >  1:218425/1‑91 (MQ=255)
                                                                          cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGTATCattatt                 >  2:1640115/1‑91 (MQ=255)
                                                                          cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGCATCATCAtt                 >  2:2481983/1‑91 (MQ=255)
                                                                          cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCAtt                 <  2:1147632/91‑1 (MQ=255)
                                                                            ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg               >  1:3820440/1‑91 (MQ=255)
                                                                            ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg               <  2:580473/91‑1 (MQ=255)
                                                                            ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg               <  2:201698/91‑1 (MQ=255)
                                                                             tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt              >  1:243572/1‑91 (MQ=255)
                                                                             tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt              >  1:242829/1‑91 (MQ=255)
                                                                             tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt              >  1:189993/1‑91 (MQ=255)
                                                                               gTAGCAGATGAAATACTTAGTATTGCCTCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt            >  2:2242061/1‑91 (MQ=255)
                                                                                tAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATc           >  2:362590/1‑91 (MQ=25)
                                                                                 aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa          >  1:1235425/1‑91 (MQ=255)
                                                                                  gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTTCTGTATCAATACGTTGGCGATTCGCTCTCACAGAATCATCATTTGTATccc         >  2:3280950/1‑89 (MQ=255)
                                                                                   ctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt        >  2:206735/3‑91 (MQ=37)
                                                                                   cAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt        >  2:624995/1‑91 (MQ=255)
                                                                                    aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAGTc       >  2:1202678/1‑91 (MQ=255)
                                                                                    aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc       <  2:901614/91‑1 (MQ=255)
                                                                                     gATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc      >  1:2177161/1‑91 (MQ=255)
                                                                                      aTGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCa     <  1:697701/91‑1 (MQ=255)
                                                                                       tGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCAc    <  1:1272326/91‑1 (MQ=255)
                                                                                         aaaTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGt  >  1:193708/1‑91 (MQ=255)
                                                                                         |                                                                                          
TTCGCTAGAAATATCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGT  >  USA300TCH1516_ALE/2517686‑2517865

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