Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,165,416 A→T 8.5% intergenic (‑384/‑94) tsaE ← / → ilvD tRNA threonylcarbamoyladenosine biosynthesis protein TsaE/Dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,165,4160AT 8.5% 582.5 / 33.0 185intergenic (‑384/‑94)tsaE/ilvDtRNA threonylcarbamoyladenosine biosynthesis protein TsaE/Dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (87/82);  new base T (5/11);  total (92/93)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TATACTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAATTT  >  USA300TCH1516_ALE/2165326‑2165506
                                                                                          |                                                                                          
tataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataa                                                                                            <  1:1610280/91‑1 (MQ=255)
tataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataa                                                                                            <  1:884526/91‑1 (MQ=255)
tataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataa                                                                                            <  1:2621642/91‑1 (MQ=255)
  taCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAataata                                                                                          >  1:1491769/1‑91 (MQ=255)
   aCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAAtattat                                                                                         <  2:3785622/91‑1 (MQ=255)
   aCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAatat                                                                                         >  2:2478164/1‑91 (MQ=255)
    ctctTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATattatt                                                                                        <  1:3014214/91‑1 (MQ=255)
     tctTTAGGACTCGGACGTTGGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTg                                                                                       <  2:1498459/91‑1 (MQ=255)
      ctTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTgg                                                                                      <  2:2870660/91‑1 (MQ=255)
       tttAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTggg                                                                                     <  2:2577402/91‑1 (MQ=255)
        ttAGGACTCGAACGTTAGTAAATATTTACTAAAGGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGa                                                                                    <  1:3446644/91‑1 (MQ=255)
        ttAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGa                                                                                    <  1:1663731/91‑1 (MQ=255)
        ttAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGa                                                                                    >  2:1397592/1‑91 (MQ=255)
        ttAGGACTCGAACGTTAGTAAATATGTCCGATACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATAgggc                                                                                    <  2:1119417/91‑2 (MQ=255)
         tAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGAc                                                                                   <  1:3219338/91‑1 (MQ=255)
           ggACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGt                                                                                 <  2:3036535/91‑1 (MQ=255)
            gACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtg                                                                                >  2:3849901/1‑90 (MQ=255)
             aCTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGttt                                                                               >  1:1615266/1‑91 (MQ=255)
             aCTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGttt                                                                               >  1:1890891/1‑91 (MQ=255)
              ctctAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGACGttta                                                                              <  2:1590396/91‑2 (MQ=255)
              cTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGACGttta                                                                              <  2:1590562/91‑2 (MQ=255)
              cTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtttt                                                                              >  1:597762/1‑91 (MQ=255)
              cTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtttt                                                                              <  2:3876013/91‑1 (MQ=255)
                cGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtttttt                                                                            <  1:1286297/91‑1 (MQ=255)
                 gAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGACGtttacaa                                                                           <  1:3529540/91‑5 (MQ=255)
                 gAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGttttttt                                                                           >  2:1174299/1‑91 (MQ=255)
                  aaCGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTa                                                                          >  2:1707025/1‑91 (MQ=255)
                   aCGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTAt                                                                         >  1:2508246/1‑91 (MQ=255)
                    cGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAATCGTCCCAATAATATTGGGACGTTTTTTTATg                                                                        <  2:918285/91‑1 (MQ=255)
                     gTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGt                                                                       <  2:2202096/91‑1 (MQ=255)
                     gTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGt                                                                       <  1:804566/91‑1 (MQ=255)
                       tAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGttt                                                                     >  2:3118413/1‑91 (MQ=255)
                       tAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGttt                                                                     >  1:2401277/1‑91 (MQ=255)
                        aGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGtttt                                                                    >  2:197278/1‑91 (MQ=255)
                           aaaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 >  2:915337/1‑91 (MQ=255)
                            aaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCt                                                                <  1:3118413/91‑1 (MQ=255)
                            aaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCt                                                                >  2:2438504/1‑91 (MQ=255)
                                  tACTAAACGCTTGAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAt                                                          >  1:2977987/1‑91 (MQ=255)
                                   aCTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAtt                                                         <  1:1513956/91‑1 (MQ=255)
                                   aCTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAtt                                                         <  2:571643/91‑1 (MQ=255)
                                    cTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTa                                                        >  1:313808/1‑91 (MQ=255)
                                     tAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTAc                                                       >  2:1809850/1‑91 (MQ=255)
                                      aaaCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACt                                                      <  1:181626/91‑1 (MQ=255)
                                      aaaCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACt                                                      <  1:1126263/91‑1 (MQ=255)
                                      aaaCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACt                                                      <  2:2143821/91‑1 (MQ=255)
                                      aaaCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACt                                                      <  1:2418503/91‑1 (MQ=255)
                                        aCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTa                                                    <  2:639426/91‑1 (MQ=255)
                                         cGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATTTTTTATCTTTCAATTACTTAt                                                   >  2:3663217/1‑91 (MQ=255)
                                         cGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAt                                                   <  2:316978/91‑1 (MQ=255)
                                         cGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAt                                                   >  1:1736294/1‑91 (MQ=255)
                                          gCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAtt                                                  >  2:43672/1‑91 (MQ=255)
                                             ttAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAttttt                                               >  2:3293712/1‑91 (MQ=255)
                                             ttAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAttttt                                               <  2:2544447/91‑1 (MQ=255)
                                              tAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTtta                                              >  2:190876/1‑91 (MQ=255)
                                              tAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTtta                                              >  2:2333292/1‑91 (MQ=255)
                                                aGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttatt                                            <  1:2722509/91‑1 (MQ=255)
                                                aGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttatt                                            >  1:2434960/1‑91 (MQ=255)
                                                aGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttatt                                            >  2:373169/1‑91 (MQ=255)
                                                 gTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttatta                                           <  1:2588665/91‑1 (MQ=255)
                                                   ccTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACt                                         <  1:2368207/91‑1 (MQ=255)
                                                    cTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACta                                        <  1:1854911/91‑1 (MQ=255)
                                                       tttCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     >  2:815691/1‑91 (MQ=255)
                                                       tttCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     >  1:2161872/1‑91 (MQ=255)
                                                       tttCTGTTTGAATGGGACATGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     <  2:2860506/91‑1 (MQ=255)
                                                        ttttGTTTGGATGGGGTTTGTAACGTTCCCAATAATTTTGGGTCTTTTTTTTTTTTTTTATCTTTCAATTATTTTTTTTTATTACTAaaaa                                    <  2:2345/88‑1 (MQ=255)
                                                          cTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAAc                                  <  2:965291/91‑1 (MQ=255)
                                                          cTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTGTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAAc                                  <  2:2987140/91‑1 (MQ=255)
                                                          cTGTCTGAATGGGACTTGTAAACGTCCCAATAATATTGGGCCGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAAc                                  >  1:1388726/1‑91 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTCTTACTATAAAACAtt                               >  2:404643/1‑90 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATg                               >  1:1587726/1‑91 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATg                               <  2:2868879/91‑1 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATg                               >  1:3797844/1‑91 (MQ=255)
                                                                  aTGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTaa                          <  1:2100087/91‑1 (MQ=255)
                                                                  aTGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTaa                          <  2:2897957/91‑1 (MQ=255)
                                                                    gggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATc                        <  2:3615393/91‑1 (MQ=255)
                                                                     ggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCa                       <  1:2791201/91‑1 (MQ=255)
                                                                     ggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCa                       >  1:952097/1‑91 (MQ=255)
                                                                      gACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAt                      >  1:1043187/1‑91 (MQ=255)
                                                                      gACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAt                      >  1:3646817/1‑91 (MQ=255)
                                                                      gACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAt                      >  1:3823377/1‑91 (MQ=255)
                                                                       aCTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAtt                     >  2:2513430/1‑91 (MQ=255)
                                                                        cTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTa                    <  2:2977987/91‑1 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTaa                          <  1:2467673/84‑1 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaa                   <  2:858841/91‑1 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaa                   <  2:1615266/91‑1 (MQ=255)
                                                                          tGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaa                  >  1:3416865/1‑91 (MQ=255)
                                                                           gTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaaa                 <  1:3728431/91‑1 (MQ=255)
                                                                           gTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaaa                 >  2:2814050/1‑91 (MQ=255)
                                                                             aaaCGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAAtt               >  2:3509910/1‑91 (MQ=255)
                                                                                cGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTAc            >  1:1077512/1‑91 (MQ=255)
                                                                                  tCCCAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:2230690/1‑91 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          <  2:322042/91‑1 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          <  2:3818879/91‑1 (MQ=255)
                                                                                    ccAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgggg        >  2:3065345/1‑91 (MQ=255)
                                                                                      aataatATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGa      >  2:360681/1‑91 (MQ=255)
                                                                                      aataatATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGa      <  2:1764129/91‑1 (MQ=255)
                                                                                       ataataTTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGaa     <  1:992378/91‑1 (MQ=255)
                                                                                        taataTTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAAt    >  1:2730768/1‑91 (MQ=255)
                                                                                          atatTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTCCTATAAAACATGATTAATCATTAAAATTTACGGGGGAAttt  >  2:226367/1‑91 (MQ=255)
                                                                                          atatTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAAttt  >  2:1295391/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATACTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAATTT  >  USA300TCH1516_ALE/2165326‑2165506

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