Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,845,292 C→T 7.9% M734I (ATG→ATA isdH ← Iron‑regulated surface determinant protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,845,2920CT 7.9% 712.7 / 35.5 227M734I (ATG→ATAisdHIron‑regulated surface determinant protein H
Reads supporting (aligned to +/- strand):  ref base C (98/111);  new base T (3/15);  total (101/126)
Fisher's exact test for biased strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 235 total aligned reads displayed.

GGCATCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATCACA  >  USA300TCH1516_ALE/1845202‑1845382
                                                                                           |                                                                                          
ggCATCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGAc                                                                                            >  2:2271926/1‑91 (MQ=255)
ggCATCTTTATCCACATTTCTATCAGTATCTTTGGCAATCTGAGTATCCATTTCTTTTAAATCATAGTTATCGAAGTGATTATTATCCGAc                                                                                            >  2:3836111/1‑91 (MQ=255)
  cATCTTTATCCACATCGCTATCAGTATCTTTGGCGCTTTGAGTATCCATTTCTTTTAGATCATATTTATCGAAGTGATTATTATCCGCCat                                                                                          >  2:3136633/1‑91 (MQ=255)
  cATCTTGATCCACATTTCTATCAGTATCTTTGGCCATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACat                                                                                          >  2:2761591/1‑91 (MQ=255)
   aTCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACata                                                                                         <  2:2286063/91‑1 (MQ=255)
    tCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACatat                                                                                        <  2:835672/91‑1 (MQ=255)
    tCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACatat                                                                                        <  2:2137883/91‑1 (MQ=255)
    tCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACatat                                                                                        >  1:2894643/1‑91 (MQ=255)
      tttATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCg                                                                                      >  1:382992/1‑91 (MQ=255)
       ttATCCACATTTCTATCAGTCTCTTTGGCAATTTGAGTATCCATTTCTCTTAAATCATATTTATCGAAGTGATTATTCTCCGCCATATAgg                                                                                     >  2:356480/1‑91 (MQ=255)
       ttATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCgg                                                                                     <  1:1030131/91‑1 (MQ=255)
       ttATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCgg                                                                                     >  2:2808961/1‑91 (MQ=255)
       ttATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCgg                                                                                     <  2:2389229/91‑1 (MQ=255)
        tATCCACATTTCTTTCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCGGa                                                                                    <  2:1005279/91‑1 (MQ=255)
        tATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCGGa                                                                                    <  1:3394246/91‑1 (MQ=255)
        tATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGa                                                                                    <  1:1915577/91‑1 (MQ=255)
         aTCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATGATATTTATCGAAGTGATTATTATCCGATATGTCGGAt                                                                                   <  1:3331952/91‑1 (MQ=255)
         aTCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCGGAt                                                                                   <  1:557616/91‑1 (MQ=255)
         aTCCACATTTCTATCAGCATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCGGAt                                                                                   <  1:3386700/91‑1 (MQ=255)
            cacaTTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGATATGTCGGATaat                                                                                <  2:2344804/91‑2 (MQ=255)
             acaTTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATTGGATAAAt                                                                               >  2:1426185/1‑91 (MQ=255)
              caTTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATg                                                                              >  1:983096/1‑91 (MQ=255)
                tttCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGAt                                                                            >  1:2935646/1‑91 (MQ=255)
                tttCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATCATTTATCCGATATGTCGGATaataat                                                                             <  1:1048980/91‑5 (MQ=25)
                  tCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCa                                                                          <  2:3377657/91‑1 (MQ=255)
                   cTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCaa                                                                         <  2:1054028/91‑1 (MQ=255)
                    tATCAGTATCTTTGGCAATTTGAGTATCCATTCCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAAc                                                                        >  2:2530389/1‑91 (MQ=255)
                      tCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAAcat                                                                      <  1:3741664/91‑1 (MQ=255)
                      tCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAAcat                                                                      >  1:485306/1‑91 (MQ=255)
                       aaGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAAGCATATTTATCGAAGTGATTATTATCCGACATGTCGGATAAATGAGCTAcatc                                                                     <  2:2108672/90‑1 (MQ=255)
                          tATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATg                                                                  <  1:1531479/91‑1 (MQ=255)
                           aTCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGt                                                                 <  2:189017/91‑1 (MQ=255)
                           aTCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGt                                                                 <  2:2434082/91‑1 (MQ=255)
                             cTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTg                                                               <  1:81289/91‑1 (MQ=255)
                             cTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTg                                                               <  1:3850385/91‑1 (MQ=255)
                               ttGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGca                                                             >  2:1407835/1‑91 (MQ=255)
                               ttGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGca                                                             >  2:2757187/1‑91 (MQ=255)
                               ttGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGca                                                             >  2:928575/1‑91 (MQ=255)
                               ttGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGACTATTACCCGACATATCGGATAAATGATCAACATCATGTTGca                                                             >  2:1278062/1‑91 (MQ=255)
                                tGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGcac                                                            >  1:2735610/1‑91 (MQ=255)
                                 ggCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGcaca                                                           >  2:716011/1‑91 (MQ=255)
                                  gCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTTTTATCCGACATATCGGATAAATGATCAACATCATGTTGCACAt                                                          <  2:3059028/91‑1 (MQ=255)
                                  gCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACAt                                                          >  2:3307548/1‑91 (MQ=255)
                                   cAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATc                                                         <  1:1655414/91‑1 (MQ=255)
                                     aTTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCtt                                                       >  2:2703632/1‑91 (MQ=255)
                                     aTTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCtt                                                       >  2:744541/1‑91 (MQ=255)
                                      tttGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCttt                                                      >  1:1287124/1‑91 (MQ=255)
                                      tttGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCttt                                                      >  2:2623953/1‑91 (MQ=255)
                                      tttGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCttt                                                      >  1:2059361/1‑91 (MQ=255)
                                       ttGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTa                                                     <  1:1807365/91‑1 (MQ=255)
                                        tGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTAt                                                    >  1:2853109/1‑91 (MQ=255)
                                         gAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATc                                                   <  2:2961669/91‑1 (MQ=255)
                                         gAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATc                                                   <  1:1232921/91‑1 (MQ=255)
                                           gTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCaa                                                 >  1:759162/1‑91 (MQ=255)
                                               ccATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAGtatt                                             <  2:1373738/91‑1 (MQ=255)
                                               ccATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAAtatt                                             <  2:1796911/91‑1 (MQ=255)
                                               ccATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAAtatt                                             <  2:2063966/91‑1 (MQ=255)
                                                 aTTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAAtattat                                           <  1:3607207/91‑1 (MQ=255)
                                                 aTTTCTTTTAAATCAGATTTATCGAGATGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAAtattat                                           <  2:3422877/91‑1 (MQ=255)
                                                  tttCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAAtattatt                                          <  2:521056/91‑1 (MQ=255)
                                                       tttAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAg                                     <  1:1180520/91‑1 (MQ=255)
                                                       tttAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAg                                     <  2:599410/91‑1 (MQ=255)
                                                         tAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCa                                   <  2:464458/91‑1 (MQ=255)
                                                         tAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCa                                   >  1:3515334/1‑91 (MQ=255)
                                                          aaaTCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCAt                                  >  2:3692213/1‑91 (MQ=255)
                                                           aaTCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATc                                 <  1:3459578/91‑1 (MQ=255)
                                                           aaTCATATTTATCGAAGTGATGATTATCCGACATATCGGATAAATGATCAACAGCATGTTGAACATCTTTATCAATATTATTATCAGCATc                                 <  2:3827190/91‑1 (MQ=255)
                                                            aTCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCGTCt                                >  1:2199924/1‑91 (MQ=255)
                                                             tCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCtt                               <  1:1023770/91‑1 (MQ=255)
                                                              cATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATCTTATTATCAGCATCttt                              >  1:1131783/1‑91 (MQ=255)
                                                              cATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCttt                              >  1:1279442/1‑91 (MQ=255)
                                                               atatTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTa                             >  1:458641/1‑91 (MQ=255)
                                                                     aTCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACg                       >  2:3608331/1‑91 (MQ=255)
                                                                      tCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGt                      <  1:3698736/91‑1 (MQ=255)
                                                                        gAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCa                    >  2:3933209/1‑91 (MQ=255)
                                                                        gAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCa                    >  2:3582828/1‑91 (MQ=255)
                                                                        gAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCa                    >  1:156066/1‑91 (MQ=255)
                                                                         aaGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAg                   <  1:870634/91‑1 (MQ=255)
                                                                         aaGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAg                   <  2:1674748/91‑1 (MQ=255)
                                                                         aaGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAg                   >  1:1479258/1‑91 (MQ=255)
                                                                         aaGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAg                   <  2:3230645/91‑1 (MQ=255)
                                                                          aGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGa                  >  1:3814903/1‑91 (MQ=255)
                                                                           gTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGAc                 <  2:1095878/91‑1 (MQ=255)
                                                                              attattATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGActct              >  1:3748587/1‑91 (MQ=255)
                                                                              attattATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGActct              <  1:2676274/91‑1 (MQ=255)
                                                                              attattATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGActct              <  2:1484763/91‑1 (MQ=255)
                                                                                tattatCCGACATATCGGATAAGTGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTgg            <  2:1629167/91‑1 (MQ=255)
                                                                                   ttatCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGtt          >  1:1165759/1‑91 (MQ=255)
                                                                                   ttatCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGtt          >  1:2524995/1‑91 (MQ=255)
                                                                                   ttatCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGtt          >  1:1725945/1‑91 (MQ=255)
                                                                                    tatCCGATATGTCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTc         >  1:2615393/1‑91 (MQ=255)
                                                                                    tatCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTc         >  1:2326703/1‑91 (MQ=255)
                                                                                      tCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTa       <  2:2947161/91‑1 (MQ=255)
                                                                                       ccGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTAt      <  1:2248222/91‑1 (MQ=255)
                                                                                       ccGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTAt      <  2:2038966/91‑1 (MQ=255)
                                                                                        cGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATc     <  1:3048274/91‑1 (MQ=255)
                                                                                         gACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATca    >  2:1824399/1‑91 (MQ=255)
                                                                                         gACATATCGGATAAATGATCAACATCATGCTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATca    <  1:3170446/91‑1 (MQ=255)
                                                                                          aCATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATcac   <  1:211318/91‑1 (MQ=255)
                                                                                          aCATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATcac   <  1:2874553/91‑1 (MQ=255)
                                                                                           |                                                                                          
GGCATCTTTATCCACATTTCTATCAGTATCTTTGGCAATTTGAGTATCCATTTCTTTTAAATCATATTTATCGAAGTGATTATTATCCGACATATCGGATAAATGATCAACATCATGTTGCACATCTTTATCAATATTATTATCAGCATCTTTAACCACGTCAGACTCTGGTTCTATCACA  >  USA300TCH1516_ALE/1845202‑1845382

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: