Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,758 Δ1 bp 13.1% intergenic (+25/‑126) btr → / → yxeP_1 HTH‑type transcriptional activator Btr/putative hydrolase YxeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE115,7580A.13.1% 658.4 / 89.6 213intergenic (+25/‑126)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base A (89/95);  new base . (22/6);  total (111/102)
Fisher's exact test for biased strand distribution p-value = 3.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 222 total aligned reads displayed.

TTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAT  >  USA300TCH1516_ALE/115668‑115848
                                                                                          |                                                                                           
ttaaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGa                                                                                             <  1:1044155/91‑1 (MQ=255)
  aaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGata                                                                                           >  1:2799608/1‑91 (MQ=255)
  aaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGata                                                                                           >  1:2184353/1‑91 (MQ=255)
   aacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGatat                                                                                          >  1:28660/1‑91 (MQ=255)
    acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                         <  2:3119664/91‑1 (MQ=255)
     cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                       <  2:1789536/91‑1 (MQ=255)
      attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                      >  2:1750198/1‑91 (MQ=255)
       ttaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACAt                                                                                     <  1:1926522/91‑1 (MQ=255)
        taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                    >  2:3308352/1‑91 (MQ=255)
        taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                    >  1:1565641/1‑91 (MQ=255)
            aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                                <  2:2678589/91‑1 (MQ=255)
             tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtg                                                                               <  2:684201/91‑1 (MQ=255)
              gATCACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                              <  1:3127089/91‑1 (MQ=255)
               aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                             >  2:1630767/1‑91 (MQ=255)
               aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTtgtg                                                                             <  1:1082539/91‑1 (MQ=37)
                 tACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAt                                                                           >  1:2083565/1‑91 (MQ=255)
                  aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtt                                                                          <  1:1205158/91‑1 (MQ=255)
                   cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                         >  2:2972106/1‑91 (MQ=255)
                   cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                         >  1:1338398/1‑91 (MQ=255)
                   cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAttt                                                                         <  1:1762629/91‑1 (MQ=25)
                    gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                        >  1:717251/1‑91 (MQ=255)
                    gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                        >  2:482781/1‑91 (MQ=255)
                    gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                        >  1:164839/1‑91 (MQ=255)
                    gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTATGTATATCAAAATTTGTGAtttg                                                                        <  1:2266492/91‑2 (MQ=25)
                      aTGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttttt                                                                      <  2:414086/91‑1 (MQ=255)
                       tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                     >  2:3519451/1‑91 (MQ=255)
                        gggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtat                                                                    >  2:3835650/1‑91 (MQ=255)
                        gggATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtat                                                                    <  2:2943428/91‑1 (MQ=255)
                           aTTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATAtt                                                                 <  2:1344044/91‑1 (MQ=255)
                             tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                               >  1:2406759/1‑91 (MQ=255)
                              cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                              >  1:285198/1‑91 (MQ=255)
                               aCTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaaa                                                             <  2:3782907/91‑1 (MQ=255)
                                cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                            >  1:2973576/1‑91 (MQ=255)
                                cTTAATTGAAGTAAAACGTCACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                            >  2:3393416/1‑91 (MQ=255)
                                  tAATTGAAGTAATACGTTACCCAAGTTCATAAAATGCTCACAAATCACATATTGTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                          <  2:2053375/91‑1 (MQ=255)
                                  tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                          <  1:2017433/91‑1 (MQ=255)
                                   aaTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                         >  1:752812/1‑91 (MQ=255)
                                     ttGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaa                                                       >  2:3333822/1‑91 (MQ=255)
                                       gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                     <  1:1061508/91‑1 (MQ=255)
                                       gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                     <  2:3635145/91‑1 (MQ=255)
                                        aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                    >  2:1878296/1‑91 (MQ=255)
                                        aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                    <  2:1698860/91‑1 (MQ=255)
                                        aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                    >  1:641914/1‑91 (MQ=255)
                                        aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                    >  1:2938174/1‑91 (MQ=255)
                                          gTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGc                                                  >  1:2382474/1‑91 (MQ=255)
                                            aaaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaa                                                >  1:76822/1‑91 (MQ=255)
                                             aaaCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaa                                               >  1:845764/1‑91 (MQ=255)
                                               aCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaa                                             <  2:394274/91‑1 (MQ=255)
                                                cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                            >  2:2156333/1‑91 (MQ=255)
                                                cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                            <  2:792544/91‑1 (MQ=255)
                                                cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                            >  1:3644793/1‑91 (MQ=255)
                                                 gTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                           >  2:2595678/1‑91 (MQ=255)
                                                 gTTACCCAAGTTCATAAAATGATCACAAATCACAAATGCGGATATACATAATTTGGGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                           >  2:1556423/1‑91 (MQ=255)
                                                  ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                          <  1:3375875/91‑1 (MQ=255)
                                                  ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATAGACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                          <  1:3881515/91‑1 (MQ=255)
                                                   tACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTa                                         <  2:1641924/91‑1 (MQ=255)
                                                    aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                        >  1:763347/1‑91 (MQ=255)
                                                      ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                      >  1:813746/1‑91 (MQ=255)
                                                          gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                  <  1:482781/91‑1 (MQ=255)
                                                           ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                 >  2:2009034/1‑91 (MQ=255)
                                                           ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                 >  1:665185/1‑91 (MQ=255)
                                                              aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                              >  1:2667391/1‑91 (MQ=255)
                                                              aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                              >  1:2667986/1‑91 (MQ=255)
                                                              aTAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGAAAAAAATTAATGGTTAACAt                              <  2:514312/91‑1 (MQ=255)
                                                               tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                             >  1:3770566/1‑91 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           <  2:540903/91‑1 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           <  2:28661/91‑1 (MQ=255)
                                                                 aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                           >  1:3433707/1‑91 (MQ=255)
                                                                  aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAAAGAAAAAAATTAATGGTTAACATctct                          <  2:2372992/91‑1 (MQ=255)
                                                                   aTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtg                         >  1:252784/1‑91 (MQ=255)
                                                                    tGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTGAAAGCTAAAATAGAAAAAAATTAATGGTTAACATCTCtgt                        <  2:3188135/91‑1 (MQ=255)
                                                                      atcacaaatcacaaATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttg                      <  1:2255277/91‑1 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:1867202/1‑91 (MQ=255)
                                                                       tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:1279180/1‑91 (MQ=255)
                                                                       taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:1966800/5‑91 (MQ=21)
                                                                       taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  2:2889556/5‑91 (MQ=21)
                                                                       taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                     >  1:1843626/5‑91 (MQ=21)
                                                                        cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                    >  2:622880/1‑91 (MQ=255)
                                                                        aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                    >  2:3717809/4‑91 (MQ=21)
                                                                        aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                    >  1:3800803/4‑91 (MQ=21)
                                                                        aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                    >  2:132486/4‑91 (MQ=21)
                                                                         acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                   <  1:2509622/91‑1 (MQ=255)
                                                                         acaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                   >  1:1728174/1‑91 (MQ=255)
                                                                         aaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                   >  2:2064120/3‑91 (MQ=21)
                                                                          caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                  <  2:1132608/91‑1 (MQ=255)
                                                                          caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                  <  1:1319543/91‑1 (MQ=255)
                                                                          caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                  >  2:2366776/1‑91 (MQ=255)
                                                                             atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa               <  2:1427893/91‑1 (MQ=255)
                                                                             atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa               <  2:1427710/91‑1 (MQ=255)
                                                                              tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat              <  2:2645429/91‑1 (MQ=255)
                                                                              tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat              >  1:1751342/1‑91 (MQ=37)
                                                                                acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat            <  1:3665898/91‑1 (MQ=255)
                                                                                 caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTGAAATTGCATAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata           <  1:3091079/91‑1 (MQ=255)
                                                                                 caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata           <  1:3497972/91‑1 (MQ=255)
                                                                                 caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata           <  2:1518161/91‑1 (MQ=255)
                                                                                  aaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataa          >  2:1405498/1‑91 (MQ=255)
                                                                                      tttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaataa      <  1:85758/91‑1 (MQ=255)
                                                                                        tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg    >  1:1188506/1‑91 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa   >  1:916536/1‑91 (MQ=255)
                                                                                         gATATACATAATTTGTGATTTTTTATATGCAAAGCTAAAATTGCAAAAAAGTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa   <  2:1648539/91‑1 (MQ=255)
                                                                                          |                                                                                           
TTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAT  >  USA300TCH1516_ALE/115668‑115848

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: