Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 428,608 C→G 5.6% G228A (GGT→GCT) ‡ ahpF ← Alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE428,6080CG 5.6% 613.3 / 19.0 179G228A (GGT→GCT) ‡ahpFAlkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (82/87);  new base G (3/7);  total (85/94)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 182 total aligned reads displayed.

TTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTCAG  >  USA300TCH1516_ALE/428519‑428698
                                                                                         |                                                                                          
ttgttTCTTTAGCAGTAATGAAGATCTCAATACCAGCAGTATCATTAACTTGGCCACCGTTACGATCAGCAACTATACCAGTACGTAAAcc                                                                                           <  2:2965766/91‑1 (MQ=255)
 tgttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCt                                                                                          <  2:2885332/91‑1 (MQ=255)
 tgttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCt                                                                                          >  1:2504178/1‑91 (MQ=255)
  gttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCtt                                                                                         >  1:2454483/1‑91 (MQ=255)
       tttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGt                                                                                    >  2:3026336/1‑91 (MQ=255)
        ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTg                                                                                   >  1:2452541/1‑91 (MQ=255)
        ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGta                                                                                   <  2:1887953/91‑2 (MQ=255)
         tAACAGTAGTGGAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGtac                                                                                  <  2:2752381/91‑3 (MQ=255)
         tAACAGTAATGAAGTTCTCAATACCAGCGGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGtac                                                                                  <  2:2807069/91‑3 (MQ=255)
         tAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGc                                                                                  <  2:1287920/91‑1 (MQ=255)
          aaCAGTAATGAAGTTCTCAATGCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCt                                                                                 <  2:1341974/91‑1 (MQ=255)
           aCAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATTCGATCAGCAACTATACCAGCACGTAAAGGTTTACGta                                                                                   <  1:3279780/88‑2 (MQ=255)
           aCAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtg                                                                                >  1:1048641/1‑91 (MQ=255)
            cAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgt                                                                               >  1:409795/1‑91 (MQ=255)
             aGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgtg                                                                              <  1:1407645/91‑1 (MQ=255)
             aGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgtg                                                                              <  1:1358139/91‑1 (MQ=255)
              gTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                             >  1:398662/1‑91 (MQ=255)
              gTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                             <  1:3528770/91‑1 (MQ=255)
               tAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTa                                                                            >  2:2507221/1‑91 (MQ=255)
                aaTGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTaa                                                                           >  2:145895/1‑91 (MQ=255)
                 aTGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTaaa                                                                          >  2:3040203/1‑91 (MQ=255)
                  tGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAAt                                                                         >  1:448645/1‑91 (MQ=255)
                  tGAAGTTCTCAATACCAGCAGTATCATTAACCTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAAt                                                                         <  1:2048445/91‑1 (MQ=255)
                    aaGTTCTCAATACCAGCAGTATCATTAACTTGGTCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCg                                                                       <  1:3172212/91‑1 (MQ=255)
                     aGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGc                                                                      >  1:1509059/1‑91 (MQ=255)
                     aGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGc                                                                      <  2:3390862/91‑1 (MQ=255)
                      gTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCt                                                                     >  1:1622959/1‑91 (MQ=255)
                      gTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCt                                                                     >  1:1622908/1‑91 (MQ=255)
                          tcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCAc                                                                 <  1:592227/91‑1 (MQ=255)
                          tcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCAc                                                                 >  1:2903174/1‑91 (MQ=255)
                           cAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACt                                                                <  1:3354063/91‑1 (MQ=255)
                           cAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACt                                                                <  2:3297149/91‑1 (MQ=255)
                               aCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCa                                                            <  1:1617200/91‑1 (MQ=255)
                               aCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCa                                                            <  2:3154836/91‑1 (MQ=255)
                                 cagcagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACt                                                          <  2:3529741/91‑1 (MQ=255)
                                 cagcagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACt                                                          <  1:1981815/91‑1 (MQ=255)
                                     agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                             >  2:3790146/1‑68 (MQ=255)
                                     agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCa                                                      >  1:782247/1‑91 (MQ=255)
                                     agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCa                                                      >  1:3689324/1‑91 (MQ=255)
                                     agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCa                                                      >  1:3805279/1‑91 (MQ=255)
                                      gTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAg                                                     <  2:664690/91‑1 (MQ=255)
                                      gTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAg                                                     <  1:1848393/91‑1 (MQ=255)
                                         tCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGac                                                  >  2:3820437/1‑91 (MQ=255)
                                          cATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGacc                                                 >  1:2282409/1‑91 (MQ=255)
                                            ttAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccac                                               >  2:831018/1‑91 (MQ=255)
                                             tAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacc                                              >  2:921248/1‑91 (MQ=255)
                                             tAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacc                                              <  2:434237/91‑1 (MQ=255)
                                              aaCTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacca                                             <  1:2724536/91‑1 (MQ=255)
                                               aCTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccac                                            <  1:1490534/91‑1 (MQ=255)
                                               aCTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccac                                            >  2:1428154/1‑91 (MQ=255)
                                                cTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccacc                                           <  1:2846295/91‑1 (MQ=255)
                                                cTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccacc                                           <  1:245860/91‑1 (MQ=255)
                                                cTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccacc                                           <  1:2373747/91‑1 (MQ=255)
                                                  tGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCaa                                         >  2:1557625/1‑91 (MQ=255)
                                                   ggCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAAc                                        <  2:2399645/91‑1 (MQ=255)
                                                    gCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACg                                       >  1:2090663/1‑91 (MQ=255)
                                                    gCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACg                                       >  2:2615863/1‑91 (MQ=255)
                                                       aCCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGCGGAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGAtt                                    >  2:3661236/1‑91 (MQ=255)
                                                        ccGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTa                                   >  2:3535506/1‑91 (MQ=255)
                                                          gATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAg                                 >  2:3052617/1‑91 (MQ=255)
                                                          gATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAg                                 >  1:485519/1‑91 (MQ=255)
                                                            tACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAc                               <  1:660230/91‑1 (MQ=255)
                                                             aCGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAca                              <  2:1187442/91‑1 (MQ=255)
                                                              cGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGCCCACCACCAACGATTAAGAcat                             >  2:789386/1‑91 (MQ=255)
                                                               gATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatc                            >  1:1043974/1‑91 (MQ=255)
                                                               gATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatc                            >  2:2121105/1‑91 (MQ=255)
                                                                aTCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatca                           >  1:3576259/1‑91 (MQ=255)
                                                                 tCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCTCTACCACTAGCAGGACCACCACCAACGATTAAGAcatcat                          <  2:3182944/91‑1 (MQ=255)
                                                                  cAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATa                         <  2:84736/91‑1 (MQ=255)
                                                                  cAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATa                         <  1:3402924/91‑1 (MQ=255)
                                                                    gCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAg                       <  2:3752021/91‑1 (MQ=255)
                                                                    gCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAg                       <  2:3532209/91‑1 (MQ=255)
                                                                     cAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAgg                      <  1:959205/91‑1 (MQ=255)
                                                                     cAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAgg                      <  1:1776270/91‑1 (MQ=255)
                                                                      aaCTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATACGGt                     >  2:591614/1‑91 (MQ=255)
                                                                        cTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTc                   <  1:1253382/91‑1 (MQ=255)
                                                                        cTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTc                   >  2:3280448/1‑91 (MQ=255)
                                                                          ataCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCtt                 >  1:2415298/1‑91 (MQ=255)
                                                                           taCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCttt                <  1:1386811/91‑1 (MQ=255)
                                                                            aCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTa               <  2:782247/91‑1 (MQ=255)
                                                                            aCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTa               <  2:3689324/91‑1 (MQ=255)
                                                                            aCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTa               >  1:3638793/1‑91 (MQ=255)
                                                                             ccAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAt              >  2:2189736/1‑91 (MQ=255)
                                                                              cAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAtt             >  2:1737266/1‑91 (MQ=255)
                                                                               agcaCGTAAAGGTTTACGTACTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAttt            >  1:744396/4‑91 (MQ=255)
                                                                               aGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAttt            >  2:3668567/1‑91 (MQ=255)
                                                                               aGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAttt            <  1:3755958/91‑1 (MQ=255)
                                                                                 tACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTc          <  2:1028402/91‑1 (MQ=255)
                                                                                  aCGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCa         <  1:1750793/91‑1 (MQ=255)
                                                                                  aCGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCa         <  2:828352/91‑1 (MQ=255)
                                                                                   cGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCaa        <  2:2026974/91‑1 (MQ=255)
                                                                                     tAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAAc      <  1:1723177/91‑1 (MQ=255)
                                                                                     tAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGCTTCTTTATTTTCAAAc      >  2:3713476/1‑91 (MQ=255)
                                                                                      aaaCCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACt     <  1:2528073/91‑1 (MQ=255)
                                                                                       aaCCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTc    <  1:1506857/91‑1 (MQ=255)
                                                                                       aaCCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTc    >  1:1408217/1‑91 (MQ=255)
                                                                                       aaCCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATAAAGACATCATAAGGTTCTTTATTTTCAAACTc    <  2:1771582/91‑1 (MQ=255)
                                                                                        aCCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTCa   <  2:741358/91‑1 (MQ=255)
                                                                                         ccTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACtagt  <  2:3214013/91‑4 (MQ=255)
                                                                                         ccTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTCAg  >  1:3228251/1‑91 (MQ=255)
                                                                                         |                                                                                          
TTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTCAAACTCAG  >  USA300TCH1516_ALE/428519‑428698

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: