Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,789 G→T 8.1% A948E (GCA→GAA)  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,7890GT 8.1% 618.1 / 27.2 187A948E (GCA→GAA) narGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base G (71/99);  new base T (8/8);  total (80/107)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

TCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACTTT  >  USA300TCH1516_ALE/2517699‑2517879
                                                                                          |                                                                                          
 cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                           <  2:1006879/91‑1 (MQ=255)
 cTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGAAAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGc                                                                                           <  2:846491/91‑1 (MQ=255)
   tttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCAt                                                                                         >  2:1758167/1‑91 (MQ=255)
    ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                                        <  1:1165544/91‑1 (MQ=38)
    ttAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATc                                                                                        <  1:3210991/91‑1 (MQ=38)
     tAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATct                                                                                       <  1:2331376/91‑2 (MQ=37)
       aCGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATctgc                                                                                     <  2:3142803/91‑4 (MQ=37)
        cGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCt                                                                                    <  2:2418777/91‑1 (MQ=255)
         ggCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTa                                                                                   <  1:1902760/91‑1 (MQ=255)
          gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                                  <  1:1783855/91‑1 (MQ=255)
             ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                               >  2:3487059/1‑91 (MQ=255)
             ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                               <  2:470859/91‑1 (MQ=255)
             ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                               <  1:3582862/91‑1 (MQ=255)
             ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                               <  1:104777/91‑1 (MQ=255)
               ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                             >  1:2046530/1‑91 (MQ=255)
                cAGTTGGTTCTTCGAGATCTTCATATGATTGTTGTGATGATTTAGCATTGGTAGCAGATGAGATACTTAGTATTGCATCAGCTACATTACg                                                                            <  2:591350/91‑1 (MQ=255)
                 aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                                           >  1:443100/1‑91 (MQ=255)
                   tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                         <  1:2010395/91‑1 (MQ=255)
                   tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                         <  1:2446355/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)
                        tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                                    <  1:3455753/91‑1 (MQ=255)
                         cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                                   <  1:3845333/91‑1 (MQ=255)
                           tcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCaa                                                                 <  1:1364742/91‑1 (MQ=255)
                             gagaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATa                                                               <  1:2622849/91‑1 (MQ=255)
                              agaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATAc                                                              <  2:2681522/91‑1 (MQ=255)
                               gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                                             <  1:3089714/91‑1 (MQ=255)
                               gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTACATCAGCTACATTACGTGCTGTATCAATACg                                                             >  1:1743721/1‑91 (MQ=255)
                                aTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGt                                                            >  2:3243496/1‑91 (MQ=255)
                                 tCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTg                                                           >  1:2163110/1‑91 (MQ=255)
                                  cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTgg                                                          >  2:143906/1‑91 (MQ=255)
                                    tCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGAc                                                        <  1:896164/91‑1 (MQ=255)
                                     cATATGATTTTTGTGATAATTTACCATTCGTATCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                       <  1:1234099/91‑1 (MQ=255)
                                      atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                                      <  2:2550137/91‑1 (MQ=255)
                                       tatGATTATTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAt                                                     <  1:1953284/91‑1 (MQ=255)
                                        atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                                    >  2:2088317/1‑91 (MQ=255)
                                        atGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGAtt                                                    <  2:2230261/91‑1 (MQ=255)
                                          gATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCg                                                  >  1:460081/1‑91 (MQ=255)
                                            tttttGTGATAATTTACCAGTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGcc                                                >  2:3451747/1‑90 (MQ=255)
                                             ttttGTGATAATTTACCATTCGTAGCAGATGAAATACTGAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctc                                               <  2:2957818/91‑1 (MQ=255)
                                              tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTTCATTACGTGCTGTATCAATACGTGGACGATTCGctct                                              <  1:3340319/91‑1 (MQ=255)
                                              tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                              <  1:2953973/91‑1 (MQ=255)
                                              tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                              <  1:2071518/91‑1 (MQ=255)
                                                tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                                            <  2:2926893/91‑1 (MQ=255)
                                                tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                                            <  1:1235559/91‑1 (MQ=255)
                                                tgtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTca                                            >  1:3638732/1‑91 (MQ=255)
                                                 gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                           <  1:3535159/91‑1 (MQ=255)
                                                 gtaaTAATTCACCATTCGTAGCAGATGAAATACTTAGTATTTCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGTTCTcac                                           >  2:2623957/4‑91 (MQ=255)
                                                    aTAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGa                                        >  2:3127790/1‑91 (MQ=255)
                                                     tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                                       <  1:1037870/91‑1 (MQ=255)
                                                     tAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGaa                                       >  2:357573/1‑91 (MQ=255)
                                                      aaTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAat                                      <  2:3845930/91‑1 (MQ=255)
                                                         ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                                   >  2:1289828/1‑91 (MQ=255)
                                                         ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                                   >  1:218425/1‑91 (MQ=255)
                                                         ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                                   <  2:2594743/91‑1 (MQ=255)
                                                         ttACCATGCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAGTACGTGCTGGATCAATACGTGGGCGATTCGCTCGCACAGAatcat                                   <  2:688112/91‑1 (MQ=255)
                                                           aCCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAAtcatca                                 <  2:2506451/91‑1 (MQ=255)
                                                             cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGTATCattatt                               >  2:1640115/1‑91 (MQ=255)
                                                             cATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCAtt                               >  1:658521/1‑91 (MQ=255)
                                                               ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg                             <  2:201698/91‑1 (MQ=255)
                                                                tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt                            >  1:242929/1‑91 (MQ=255)
                                                                tCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGt                            >  1:189993/1‑91 (MQ=255)
                                                                 cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                           >  1:68692/1‑91 (MQ=255)
                                                                  gTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt                          >  1:1322373/1‑91 (MQ=25)
                                                                  gTAGCAGATGAAATACTTAGTATTGCCTCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt                          >  2:2242061/1‑91 (MQ=255)
                                                                   tAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATc                         >  2:362590/1‑91 (MQ=25)
                                                                    agctgATGCAATACTAAGTATTTCATCTGTTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                        >  2:1213026/5‑91 (MQ=25)
                                                                     gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTTCTGTATCAATACGTTGGCGATTCGCTCTCACAGAATCATCATTTGTATccc                       >  2:3280950/1‑89 (MQ=255)
                                                                      cAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt                      >  2:624995/1‑91 (MQ=255)
                                                                       aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc                     <  2:901614/91‑1 (MQ=255)
                                                                       aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc                     >  2:1202659/1‑91 (MQ=255)
                                                                           gcaaTACTAAGTATTTCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg                 >  1:3754757/3‑91 (MQ=255)
                                                                              aTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAc              >  1:2985567/1‑91 (MQ=255)
                                                                                aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGCACCATCATTTGTATCACTCCACGTAccc            >  2:2318946/1‑90 (MQ=255)
                                                                                aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt            >  1:299057/1‑91 (MQ=255)
                                                                                aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt            <  2:3763016/91‑1 (MQ=255)
                                                                                aCTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCACACAGAATCATCATTTGTATCACACCACGTACCt            <  2:1107046/91‑1 (MQ=255)
                                                                                  ttAGTATTGGATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa          <  1:3787223/91‑1 (MQ=255)
                                                                                  ttAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa          >  1:2541822/1‑91 (MQ=255)
                                                                                   tAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAc         >  2:2138097/1‑91 (MQ=255)
                                                                                   tAGCATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAAc         <  1:3062860/91‑1 (MQ=255)
                                                                                    aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTCCCTcaca        >  2:341414/1‑91 (MQ=255)
                                                                                    aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACa        >  1:3268603/1‑91 (MQ=255)
                                                                                    aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACa        <  1:1671496/91‑1 (MQ=255)
                                                                                    aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACa        >  2:3816499/1‑91 (MQ=255)
                                                                                    aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACa        <  1:2912822/91‑1 (MQ=255)
                                                                                     gTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACAt       >  1:3399240/1‑91 (MQ=255)
                                                                                     gTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACAt       <  2:1828195/91‑1 (MQ=255)
                                                                                     gTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACAt       <  1:2812951/91‑1 (MQ=255)
                                                                                      tATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATa      <  1:239184/91‑1 (MQ=255)
                                                                                        ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt    <  1:1610488/91‑1 (MQ=255)
                                                                                        ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt    <  2:3655258/91‑1 (MQ=255)
                                                                                        ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt    <  1:1353795/91‑1 (MQ=255)
                                                                                         tGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACtt   <  1:3199425/91‑1 (MQ=255)
                                                                                         tGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACtt   <  1:1456116/91‑1 (MQ=255)
                                                                                          gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt  <  1:2169510/91‑1 (MQ=255)
                                                                                          gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt  <  1:965302/91‑1 (MQ=255)
                                                                                          gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt  >  2:2414753/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCTTTTAACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACTTT  >  USA300TCH1516_ALE/2517699‑2517879

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