Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,517,794 A→T 7.0% A946A (GCT→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,7940AT 7.0% 530.7 / 20.0 157A946A (GCT→GCAnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base A (79/67);  new base T (9/2);  total (88/69)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 158 total aligned reads displayed.

AACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACTTTTTAAT  >  USA300TCH1516_ALE/2517705‑2517884
                                                                                         |                                                                                          
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:906308/91‑1 (MQ=25)
    gCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCTGATGCAATACTAAGTATTTCATCTGCTAc                                                                                       <  1:906258/91‑1 (MQ=25)
     cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                                      >  2:3361103/1‑91 (MQ=255)
     cATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACa                                                                                      >  2:2179606/1‑91 (MQ=255)
      aTTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAt                                                                                     >  1:1700113/1‑91 (MQ=255)
       ttCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACAtt                                                                                    >  1:3423790/1‑91 (MQ=255)
         ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                                  <  2:1905257/91‑1 (MQ=255)
         ccAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                                  >  1:925244/1‑91 (MQ=255)
          cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                                                 >  2:645123/1‑91 (MQ=255)
          cAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACg                                                                                 >  1:594561/1‑91 (MQ=255)
           aGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGTTATTTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGt                                                                                >  2:1481574/1‑91 (MQ=255)
            gTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTg                                                                               <  1:2432913/91‑1 (MQ=255)
             tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                              <  2:936215/91‑1 (MQ=255)
             tttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGc                                                                              <  1:2349036/91‑1 (MQ=255)
              ttGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCt                                                                             >  1:1387906/1‑91 (MQ=255)
               tGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTg                                                                            >  2:1566952/1‑91 (MQ=255)
                gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                                          >  2:2680916/1‑76 (MQ=255)
                gTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGc                                                                                          <  1:2680916/76‑1 (MQ=255)
                  tcttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTAt                                                                         <  2:846315/91‑1 (MQ=255)
                   cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTAc                                                                                  >  2:2763052/1‑81 (MQ=255)
                   cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                                        >  2:2169830/1‑91 (MQ=255)
                   cttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATc                                                                        <  2:389699/91‑1 (MQ=255)
                    ttcGAGATCTTCATATGATTTTTGTTATTTTTTATCATTTGTAGCAGATTAAATACTTAGTATTTCTTCATCTACATTACGTGCTGTATCa                                                                       >  2:3620431/1‑91 (MQ=255)
                    ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                                                       <  1:102840/91‑1 (MQ=255)
                    ttcGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCa                                                                       >  2:3498925/1‑91 (MQ=255)
                        agaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATAc                                                                   <  2:3754042/91‑1 (MQ=255)
                         gaTCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACg                                                                  <  2:2008346/91‑1 (MQ=255)
                            cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                                       <  1:2297688/67‑1 (MQ=255)
                            cTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTAc                                                                                       >  2:2297688/1‑67 (MQ=255)
                               cATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                            <  1:745525/91‑1 (MQ=255)
                               cATATAATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACg                                                            <  2:173125/91‑1 (MQ=255)
                                atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                                           >  1:187602/1‑91 (MQ=255)
                                atatGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGa                                                           <  1:1661531/91‑1 (MQ=255)
                                   tGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTc                                                        <  1:190619/91‑1 (MQ=255)
                                      tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                                                     <  1:2686428/91‑1 (MQ=255)
                                      tttttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGct                                                     >  1:772216/1‑91 (MQ=255)
                                        tttGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                                   >  2:756990/1‑91 (MQ=255)
                                           gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                                >  1:779216/1‑91 (MQ=255)
                                           gtgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcac                                                <  1:982790/91‑1 (MQ=255)
                                            tgATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTcaca                                               >  2:2787411/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)
                                                  tttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatca                                         >  2:3608573/1‑91 (MQ=255)
                                                   ttACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcat                                        <  2:2629801/91‑1 (MQ=255)
                                                    tACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAatcatc                                       <  2:878240/91‑1 (MQ=255)
                                                      ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                                     >  2:2232875/1‑91 (MQ=255)
                                                      ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                                     >  2:207829/1‑91 (MQ=255)
                                                      ccATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATcatcat                                     >  2:780935/1‑91 (MQ=255)
                                                         ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGctct                                                   <  1:1016991/74‑1 (MQ=255)
                                                         ttCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTg                                  <  1:2126599/91‑1 (MQ=255)
                                                           tgTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                                >  2:2239666/2‑91 (MQ=25)
                                                           cGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTa                                <  1:503071/91‑1 (MQ=255)
                                                            gTAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTAt                               >  2:1144664/1‑91 (MQ=25)
                                                             tAGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATc                              >  2:2540245/1‑91 (MQ=25)
                                                              agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                             >  2:3402497/5‑91 (MQ=25)
                                                              agctgATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                             >  1:2749930/5‑91 (MQ=25)
                                                              aGCTGATGCAATACTAAGTATTTCATCTGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                             >  2:312362/1‑91 (MQ=25)
                                                              aGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCa                             >  2:707178/1‑91 (MQ=255)
                                                               gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAc                            >  1:141054/1‑91 (MQ=255)
                                                               gCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCAc                            <  1:1485582/91‑1 (MQ=255)
                                                                cAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACt                           >  2:1074957/1‑91 (MQ=255)
                                                                 aGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTc                          <  1:2502179/91‑1 (MQ=255)
                                                                  gATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTcc                         >  1:472795/1‑91 (MQ=255)
                                                                   aTGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCa                        <  1:2827484/91‑1 (MQ=255)
                                                                     gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg                      >  1:3855087/1‑91 (MQ=255)
                                                                     gAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACg                      >  1:1384204/1‑91 (MQ=255)
                                                                       aaTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGCATCATCATTTGTATCACTCCACGTa                    >  2:3522007/1‑91 (MQ=255)
                                                                       aaTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTa                    <  1:3556926/91‑1 (MQ=255)
                                                                       aaTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTa                    <  2:3352646/91‑1 (MQ=255)
                                                                        aTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAc                   <  1:960759/91‑1 (MQ=255)
                                                                        aTACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCGCTCACAGAATCATCATTTGTATCACTCCACGTAc                   <  2:1115981/91‑1 (MQ=255)
                                                                         tACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAcc                  >  1:1165685/1‑91 (MQ=255)
                                                                         tACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTAcc                  >  1:1155694/1‑91 (MQ=255)
                                                                           cTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTa                >  1:635816/1‑91 (MQ=255)
                                                                            ttAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa               >  1:3801836/1‑91 (MQ=255)
                                                                            ttAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTaa               <  2:3017289/91‑1 (MQ=255)
                                                                              aGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACa             <  1:1900512/91‑1 (MQ=255)
                                                                               gTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACAt            <  2:3688793/91‑1 (MQ=255)
                                                                                  ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt         >  2:3247351/1‑91 (MQ=255)
                                                                                  ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt         <  1:130497/91‑1 (MQ=255)
                                                                                  ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt         <  2:235711/91‑1 (MQ=255)
                                                                                  ttGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACt         <  1:1284715/91‑1 (MQ=255)
                                                                                    gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt       <  1:1854316/91‑1 (MQ=255)
                                                                                    gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt       >  2:815111/1‑91 (MQ=255)
                                                                                    gCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACttt       >  1:1354073/1‑91 (MQ=255)
                                                                                     cATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCt                 <  2:541101/80‑1 (MQ=255)
                                                                                     cATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACtttt      >  1:3733953/1‑91 (MQ=255)
                                                                                     cATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACtttt      >  1:348016/1‑91 (MQ=255)
                                                                                     cATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACtttt      >  1:682425/1‑91 (MQ=255)
                                                                                         aGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACTTTTTAAt  <  2:2086572/91‑1 (MQ=255)
                                                                                         |                                                                                          
AACGGCATTCCAGTTTGTTCTTCGAGATCTTCATATGATTTTTGTGATAATTTACCATTCGTAGCAGATGAAATACTTAGTATTGCATCAGCTACATTACGTGCTGTATCAATACGTGGACGATTCGCTCTCACAGAATCATCATTTGTATCACTCCACGTACCTAACATACTTTTTAAT  >  USA300TCH1516_ALE/2517705‑2517884

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: