Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,561,495 Δ1 bp 5.7% coding (546/807 nt) cpdA ← 3',5'‑cyclic adenosine monophosphate phosphodiesterase CpdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,561,4950A. 5.7% 641.4 / 28.9 193coding (546/807 nt)cpdA3',5'‑cyclic adenosine monophosphate phosphodiesterase CpdA
Reads supporting (aligned to +/- strand):  ref base A (96/86);  new base . (2/9);  total (98/95)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 208 total aligned reads displayed.

AAAGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAATTTTGATCT  >  USA300TCH1516_ALE/2561406‑2561585
                                                                                         |                                                                                           
aaaGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCacaa                                                                                            <  2:942950/91‑1 (MQ=255)
aaaGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCacaa                                                                                            <  1:2384307/91‑1 (MQ=255)
 aaGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCacaac                                                                                           <  2:807277/91‑1 (MQ=255)
  aGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCacaaca                                                                                          <  2:1065173/91‑1 (MQ=255)
   ggATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTcacacatg                                                                                         <  1:3356845/91‑5 (MQ=255)
    gATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACAtg                                                                                        >  1:3212560/1‑91 (MQ=255)
     aTTGAAATCACTTGTCCCAATAAATGCATTATAAGAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgtt                                                                                     <  2:1608287/91‑2 (MQ=255)
     aTTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACAtgt                                                                                      <  2:3295987/91‑1 (MQ=255)
     aTTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACAtgt                                                                                      >  1:2790324/1‑91 (MQ=255)
     aTTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgtt                                                                                     <  1:3639292/91‑2 (MQ=255)
     aTTGAAATCACTTGTCCCAATAAATGCATTATAAAAAGCGAATCTACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgtt                                                                                     <  1:3945343/91‑2 (MQ=255)
      ttGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgttg                                                                                    <  2:2871283/91‑3 (MQ=255)
      ttGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgttg                                                                                    <  2:322745/91‑3 (MQ=255)
      ttGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACACAtgttg                                                                                    <  1:2578139/91‑3 (MQ=255)
      ttGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGCCACAACAtgtg                                                                                     >  2:1992814/1‑91 (MQ=255)
        gAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACAtgtgtg                                                                                   >  1:796629/1‑91 (MQ=255)
         aaaTCACTTGTCCCAATAAATGCATTATAGAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGa                                                                                  <  1:1444193/91‑1 (MQ=255)
         aaaTCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGa                                                                                  >  2:1302708/1‑91 (MQ=255)
         aaaTCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGa                                                                                  >  2:1274744/1‑91 (MQ=255)
         aaaTCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGa                                                                                  >  1:497940/1‑91 (MQ=255)
          aaTCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGAc                                                                                 >  1:226435/1‑91 (MQ=255)
            tCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTa                                                                               <  1:2789221/91‑1 (MQ=255)
             cACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTaa                                                                              <  1:2308321/91‑1 (MQ=255)
              aCTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTaaa                                                                             <  2:2490604/91‑1 (MQ=255)
              aCTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTaaa                                                                             <  2:3070526/91‑1 (MQ=255)
               cTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTaaaa                                                                            <  2:3385320/91‑1 (MQ=255)
               cTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTaaaa                                                                            <  1:821066/91‑1 (MQ=255)
                   tCCCAATAAATGCATTATAAAAATTGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCAc                                                                        <  2:3679613/91‑1 (MQ=255)
                   tCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCAc                                                                        <  1:2720378/91‑1 (MQ=255)
                     ccAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAGTCAGGGTGCGTCACAACATGTGTGACTAAAATCACTc                                                                      <  1:3346163/91‑1 (MQ=255)
                     ccAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTc                                                                      >  1:1690357/1‑91 (MQ=255)
                      cAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCg                                                                     <  2:2015741/91‑1 (MQ=255)
                      cAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCg                                                                     <  1:439584/91‑1 (MQ=255)
                      cAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCg                                                                     >  2:1557313/1‑91 (MQ=255)
                         tAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCt                                                                  <  2:1703851/91‑1 (MQ=255)
                            aTGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTAtt                                                               <  1:3252597/91‑1 (MQ=255)
                            aTGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTAtt                                                               <  1:1619461/91‑1 (MQ=255)
                             tGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTAttt                                                              >  1:2509434/1‑91 (MQ=255)
                             tGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTAttt                                                              <  2:2832083/91‑1 (MQ=255)
                              gCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTc                                                             >  1:3726638/1‑91 (MQ=255)
                               cATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTcc                                                            >  2:2969865/1‑91 (MQ=255)
                                aTTATAAGAAGCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCt                                                           <  2:711303/91‑1 (MQ=255)
                                  tataAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTAc                                                         >  2:3613769/1‑91 (MQ=255)
                                    taAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACtt                                                       >  2:3551604/1‑91 (MQ=255)
                                     aaaaaTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTc                                                      >  2:2045615/1‑91 (MQ=255)
                                     aaaaaTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTc                                                      >  1:1186603/1‑91 (MQ=255)
                                     aaaaaTCGAATATACGATGCGGCATAGCAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTc                                                      >  2:84416/1‑91 (MQ=255)
                                      aaaaTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCt                                                     >  2:1898777/1‑91 (MQ=255)
                                        aaTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTaa                                                   <  2:1440678/91‑1 (MQ=255)
                                        aaTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCATACTTCTaa                                                   >  1:3777270/1‑91 (MQ=255)
                                         aTCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAat                                                  >  2:409906/1‑91 (MQ=255)
                                            gAATATACGATGCGGCATAGGAACAATGAAATCTGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAatata                                               <  1:3235403/91‑1 (MQ=255)
                                             aaTATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAatatat                                              >  1:2064537/1‑91 (MQ=255)
                                              atatACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtc                                             >  1:3446398/1‑91 (MQ=255)
                                              atatACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtc                                             <  2:1098419/91‑1 (MQ=255)
                                               tataCGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtct                                            >  2:3137231/1‑91 (MQ=255)
                                                ataCGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtctt                                           <  2:103270/91‑1 (MQ=255)
                                                 taCGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtcttt                                          >  2:2990379/1‑91 (MQ=255)
                                                   cGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACGTGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtctttct                                        <  2:3472937/91‑1 (MQ=255)
                                                    gATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATAtctttctt                                       >  1:1077627/1‑91 (MQ=255)
                                                     aTGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTttcttc                                      <  2:1797957/91‑1 (MQ=255)
                                                      tGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCa                                     <  1:2620402/91‑1 (MQ=255)
                                                      tGCCGCATAGGAACAATGAAATCAGAGTGCGTCACAACATCTGTGACTAAAATCACTCCTCTGTTTCCTACTTCTAATATATCTTTCTTCa                                     <  2:642400/91‑1 (MQ=255)
                                                          gCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTg                                 <  1:3833368/91‑1 (MQ=255)
                                                           cATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGt                                >  1:1126957/1‑91 (MQ=255)
                                                           cATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGt                                >  1:309056/1‑91 (MQ=255)
                                                            aTAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGtt                               <  2:196634/91‑1 (MQ=255)
                                                             tAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTc                              >  2:889933/1‑91 (MQ=255)
                                                              agtaaCAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCa                             <  2:2885127/88‑1 (MQ=255)
                                                              aGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCa                             >  2:2048660/1‑91 (MQ=255)
                                                               ggAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAg                            <  2:78788/91‑1 (MQ=255)
                                                               ggAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAg                            <  1:339167/91‑1 (MQ=255)
                                                                gAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGc                           <  1:922811/91‑1 (MQ=255)
                                                                 aaCAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGcc                          <  2:1200134/91‑1 (MQ=255)
                                                                 aaCAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGcc                          <  2:2658377/91‑1 (MQ=255)
                                                                 aaCAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGcc                          <  1:215807/91‑1 (MQ=255)
                                                                       gAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTAtt                    >  2:1378169/1‑91 (MQ=255)
                                                                         aaTCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTAtttt                  >  2:488631/1‑91 (MQ=255)
                                                                          aTCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTa                 >  1:183729/1‑91 (MQ=255)
                                                                           tCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAg                >  1:1583490/1‑91 (MQ=255)
                                                                           tCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATCTTAg                <  1:3127968/91‑1 (MQ=255)
                                                                            cAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGa               <  1:1849824/91‑1 (MQ=255)
                                                                            cAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGa               <  1:1189791/91‑1 (MQ=255)
                                                                            cAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGa               <  1:2711165/91‑1 (MQ=255)
                                                                            cAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATAGATCTTTCTTCACTTGTTCAGCCGCTATTTTAGa               >  2:2131740/1‑91 (MQ=255)
                                                                              gggTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATa             >  2:3339923/1‑91 (MQ=255)
                                                                              gggTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATa             <  1:3395205/91‑1 (MQ=255)
                                                                               ggTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATaa            >  1:1423739/1‑91 (MQ=255)
                                                                                tttagTCACACATGTTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACCTGTTCAGCCGCTATTTCAGATaaa           >  2:3847176/5‑91 (MQ=255)
                                                                                gTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATaaa           <  1:719927/91‑1 (MQ=255)
                                                                                gTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATaaa           <  1:887165/91‑1 (MQ=255)
                                                                                gTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATaaa           <  2:892020/91‑1 (MQ=255)
                                                                                  tagTCACACATGTTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAtt         >  1:2000792/3‑91 (MQ=255)
                                                                                  gCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAtt         >  1:1812474/1‑91 (MQ=255)
                                                                                  gCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAtt         <  1:1316401/91‑1 (MQ=255)
                                                                                   cGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAttt        >  1:622505/1‑91 (MQ=255)
                                                                                   cGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAttt        <  2:2930974/91‑1 (MQ=255)
                                                                                    gTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAAtttt       >  1:2920559/1‑91 (MQ=255)
                                                                                    gTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGAGAAAtttt       >  1:1955519/1‑91 (MQ=255)
                                                                                        caacaTGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAATTTTGATc   <  1:659999/91‑1 (MQ=255)
                                                                                         |                                                                                           
AAAGGATTGAAATCACTTGTCCCAATAAATGCATTATAAAAATCGAATATACGATGCGGCATAGGAACAATGAAATCAGGGTGCGTCACAACATGTGTGACTAAAATCACTCGTCTATTTCCTACTTCTAATATATCTTTCTTCACTTGTTCAGCCGCTATTTTAGATAAATTTTGATCT  >  USA300TCH1516_ALE/2561406‑2561585

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: