Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 979,355 A→C 5.3% intergenic (+25/+34) clpB_1 → / ← ampR Chaperone protein ClpB/HTH‑type transcriptional activator AmpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE979,3550AC 5.3% 740.5 / 23.3 215intergenic (+25/+34)clpB_1/ampRChaperone protein ClpB/HTH‑type transcriptional activator AmpR
Reads supporting (aligned to +/- strand):  ref base A (90/113);  new base C (8/4);  total (98/117)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCAATATTTC  >  USA300TCH1516_ALE/979266‑979445
                                                                                         |                                                                                          
  aTTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATc                                                                                         <  1:1885925/91‑1 (MQ=255)
   ttAAAGTTAATTTAAATTCAGACAATAGCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCg                                                                                        >  1:3516142/1‑91 (MQ=255)
   ttAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCg                                                                                        <  2:3337980/91‑1 (MQ=255)
     aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAACATTGAACCAAGAATGTGATGATTCATCGCa                                                                                      >  2:3879350/1‑91 (MQ=255)
     aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCa                                                                                      <  1:2287260/91‑1 (MQ=255)
     aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCa                                                                                      >  2:1202822/1‑91 (MQ=255)
      aaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCAt                                                                                     >  1:261629/1‑91 (MQ=255)
      aaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCAt                                                                                     >  1:3680397/1‑91 (MQ=255)
       aGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCAtt                                                                                    >  1:1068102/1‑91 (MQ=255)
        gTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                                   >  2:2624815/1‑91 (MQ=255)
        gTTAATGTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                                   <  2:349936/91‑1 (MQ=255)
         ttAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCt                                                                                  <  1:3165372/91‑1 (MQ=255)
         ttAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCt                                                                                  <  1:1868851/91‑1 (MQ=255)
          tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                                 <  2:1875868/91‑1 (MQ=255)
          tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                                 >  2:1625031/1‑91 (MQ=255)
          tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                                 >  1:3457542/1‑91 (MQ=255)
          tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                                 >  1:3458093/1‑91 (MQ=255)
          tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGCGATTAATCATCACATTCtt                                                                                 <  1:2404622/91‑1 (MQ=25)
           aaTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTg                                                                                >  1:2872509/1‑91 (MQ=255)
            aTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTgg                                                                               >  2:132557/1‑91 (MQ=255)
            aTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTgg                                                                               >  2:2650894/1‑91 (MQ=255)
              ttAAATTCAGACAATAACTTAACGTGTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGCGATTAATCATCACATTCTTGGtt                                                                             <  2:3528352/91‑1 (MQ=25)
               tAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGttt                                                                            >  1:629545/1‑91 (MQ=255)
                aaaTTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGCGATTAATCATCACATTCTTgggtgc                                                                           >  1:2053205/1‑87 (MQ=25)
                  aTTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGtttttt                                                                         <  1:1927237/91‑1 (MQ=255)
                    tCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGTATTCTTGGTTTTTTAg                                                                       >  2:2975379/1‑91 (MQ=255)
                    tCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAg                                                                       <  2:2994989/91‑1 (MQ=255)
                     cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                                      <  1:769046/91‑1 (MQ=255)
                     cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                                      <  1:2572430/91‑1 (MQ=255)
                     cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                                      >  1:1815759/1‑91 (MQ=255)
                     cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                                      <  2:1055089/91‑1 (MQ=255)
                      aGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAt                                                                     <  2:3667561/91‑1 (MQ=255)
                       gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAtt                                                                    >  2:2023466/1‑91 (MQ=255)
                       gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATCCTTGGTTTTTTAGAtt                                                                    >  2:816937/1‑91 (MQ=255)
                        aCAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAttt                                                                   <  2:2434405/91‑1 (MQ=255)
                        aCAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAttt                                                                   <  2:383154/91‑1 (MQ=255)
                          aaTAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGa                                                                 <  1:3704486/91‑1 (MQ=255)
                            tAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGaaa                                                               >  2:1983815/1‑91 (MQ=255)
                              aCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAtt                                                             <  2:211840/91‑1 (MQ=255)
                              aCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAtt                                                             >  1:1301944/1‑91 (MQ=255)
                               cTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAttt                                                            >  2:1087373/1‑91 (MQ=255)
                                ttAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATtta                                                           >  2:26198/1‑91 (MQ=255)
                                 tAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATttat                                                          >  1:3203790/1‑91 (MQ=255)
                                 tAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATttat                                                          <  1:1430497/91‑1 (MQ=255)
                                      tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                                     <  2:2454446/91‑1 (MQ=255)
                                      tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                                     <  2:601060/91‑1 (MQ=255)
                                      tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                                     <  2:3203790/91‑1 (MQ=255)
                                        tAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAAtt                                                   <  1:1202822/91‑1 (MQ=255)
                                          atgatgAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATGTATTATAATTGt                                                 <  1:3000004/91‑1 (MQ=255)
                                          aTGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGt                                                 <  2:578542/91‑1 (MQ=255)
                                           tGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGtt                                                >  2:3818851/1‑91 (MQ=255)
                                            gTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTc                                               <  2:3389346/91‑1 (MQ=255)
                                             ttGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCa                                              >  1:3247916/1‑91 (MQ=255)
                                              tGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCAt                                             >  1:1650954/1‑91 (MQ=255)
                                               gAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCAtt                                            <  2:2127858/91‑1 (MQ=255)
                                                aaaaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTg                                           >  1:2066793/1‑91 (MQ=255)
                                                aaaaaaTTCATGAAGAAAATTGACACAAGAATGTGATGATTAAACGCATTTTTGTTTTTTTAGATTTGAAATTTATTATAATTTTTCATTg                                           <  2:3608245/91‑1 (MQ=255)
                                                 aaaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTgg                                          <  1:729804/91‑1 (MQ=255)
                                                  aaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGc                                         >  2:37483/1‑91 (MQ=255)
                                                  aaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGc                                         <  1:941435/91‑1 (MQ=255)
                                                    aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                                       <  1:2440058/91‑1 (MQ=255)
                                                    aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                                       <  1:1533812/91‑1 (MQ=255)
                                                    aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                                       <  1:538570/91‑1 (MQ=255)
                                                    aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATCCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                                       <  1:3814546/91‑1 (MQ=255)
                                                     aTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAat                                      <  1:2023466/91‑1 (MQ=255)
                                                     aTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAat                                      >  1:1591296/1‑91 (MQ=255)
                                                      ttCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAata                                     <  2:1408132/91‑1 (MQ=255)
                                                       tCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAatat                                    <  1:2485544/91‑1 (MQ=255)
                                                        cATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAAtata                                   <  1:2851092/91‑1 (MQ=255)
                                                        cATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAAtata                                   <  1:2296416/91‑1 (MQ=255)
                                                         aTGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATAc                                  >  1:2858715/1‑91 (MQ=255)
                                                         aTGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATAc                                  <  1:1170455/91‑1 (MQ=255)
                                                         aTGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATAc                                  <  2:22017/91‑1 (MQ=255)
                                                         aTGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATAc                                  >  1:1484549/1‑91 (MQ=255)
                                                         aTGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATAc                                  >  1:3301156/1‑91 (MQ=255)
                                                          tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                                 >  1:2918468/1‑91 (MQ=255)
                                                          tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                                 >  1:659673/1‑91 (MQ=255)
                                                          tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                                 <  2:2333799/91‑1 (MQ=255)
                                                           gAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtt                                >  1:2006055/1‑91 (MQ=255)
                                                             aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                              <  1:2775404/91‑1 (MQ=255)
                                                             aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                              <  1:889553/91‑1 (MQ=255)
                                                             aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                              <  1:170902/91‑1 (MQ=255)
                                                                aaaTTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTaa                           <  1:510872/91‑1 (MQ=255)
                                                                       aCCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt                    <  1:1047477/91‑1 (MQ=255)
                                                                       aCCAAGAATGCGATTAATCATCACATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt                    >  2:1910702/1‑91 (MQ=25)
                                                                       aCCAAGAATGCGATTAATCATCACATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt                    >  1:2172801/1‑91 (MQ=25)
                                                                         cAAGAATGCGATTAATCATCACATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaa                  >  1:1251606/1‑91 (MQ=25)
                                                                          aaGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaat                 >  2:1762791/1‑91 (MQ=255)
                                                                          aaGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaat                 <  2:2372287/91‑1 (MQ=255)
                                                                           aGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaata                <  2:3301156/91‑1 (MQ=255)
                                                                             aaTGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTAATATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATaaa              >  2:3477801/1‑91 (MQ=255)
                                                                               tgtgATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAaca            <  2:3053664/91‑1 (MQ=255)
                                                                               tgtgATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAaca            <  2:2754126/91‑1 (MQ=255)
                                                                                 tgatgaTTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACAcc          <  1:442045/91‑1 (MQ=255)
                                                                                   atgatTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCaa        >  2:1584625/1‑91 (MQ=255)
                                                                                   atgatTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCaa        >  2:2859644/1‑91 (MQ=255)
                                                                                   atgatTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCaa        >  1:157455/1‑91 (MQ=255)
                                                                                    tgatTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCAat       <  2:2880479/91‑1 (MQ=255)
                                                                                    tgatTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCAat       <  1:3334093/91‑1 (MQ=255)
                                                                                         aaTCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCAATATTTc  <  2:487525/91‑1 (MQ=255)
                                                                                         |                                                                                          
CGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACCAATATTTC  >  USA300TCH1516_ALE/979266‑979445

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: