Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,776,634 C→T 6.7% intergenic (‑111/+40) clpX ← / ← tig ATP‑dependent Clp protease ATP‑binding subunit ClpX/Trigger factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,776,6340CT 6.7% 387.5 / 12.8 119intergenic (‑111/+40)clpX/tigATP‑dependent Clp protease ATP‑binding subunit ClpX/Trigger factor
Reads supporting (aligned to +/- strand):  ref base C (61/50);  new base T (4/4);  total (65/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 122 total aligned reads displayed.

TTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTTAATAAAT  >  USA300TCH1516_ALE/1776544‑1776723
                                                                                           |                                                                                         
ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGc                                                                                           >  1:2737636/1‑91 (MQ=255)
ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGc                                                                                           >  1:4020540/1‑91 (MQ=255)
ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACAGAAGCTATGc                                                                                           >  1:3727613/1‑91 (MQ=255)
      aCTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGc                                                                                     >  1:1464725/1‑91 (MQ=255)
       cTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATGCATAGCTTATGTATCTGCt                                                                                    <  1:2553973/91‑1 (MQ=37)
        tataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTa                                                                                   <  1:2285280/91‑1 (MQ=255)
         ataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTAc                                                                                  <  1:1901768/91‑1 (MQ=255)
         ataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTAc                                                                                  <  1:1901548/91‑1 (MQ=255)
         ataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTAc                                                                                  >  1:2588286/1‑91 (MQ=255)
         ataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTAc                                                                                  <  1:1901575/91‑1 (MQ=255)
          taCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACt                                                                                 <  1:2567847/91‑1 (MQ=255)
           aCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATGCATAGCTTATGTATCTGCTACta                                                                                <  1:3591788/91‑2 (MQ=25)
           aCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATGCATAGCTTATGTATCTGCTACta                                                                                <  1:1885678/91‑2 (MQ=25)
           aCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATGCATAGCTTATGTATCTGCTACta                                                                                <  1:2756528/91‑2 (MQ=25)
            cTAGATTGGATGTTACAATGCAACATATTAACGTACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCa                                                                               <  1:3522633/91‑1 (MQ=255)
             tAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAt                                                                              >  1:1843190/1‑91 (MQ=255)
             tAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAt                                                                              >  1:1843039/1‑91 (MQ=255)
               gATTGGATGCTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTGCTCAttt                                                                            >  1:2085918/1‑91 (MQ=255)
                  tGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTg                                                                         <  1:1328793/91‑1 (MQ=255)
                    gATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTa                                                                       >  1:1362183/1‑91 (MQ=255)
                     aTGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTaa                                                                      >  1:2377285/1‑91 (MQ=255)
                      tGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAat                                                                     >  1:2989339/1‑91 (MQ=255)
                      tGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAat                                                                     >  1:3677935/1‑91 (MQ=255)
                       gTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatt                                                                    <  1:100428/91‑1 (MQ=255)
                          aCAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaa                                                                 >  1:910268/1‑91 (MQ=255)
                           cAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaat                                                                <  1:2255699/91‑1 (MQ=255)
                           cAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaat                                                                >  1:993225/1‑91 (MQ=255)
                            aaTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaata                                                               <  1:871700/91‑1 (MQ=255)
                            aaTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaata                                                               <  1:3151680/91‑1 (MQ=255)
                             aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                              >  1:3243544/1‑91 (MQ=255)
                              tGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatt                                                             <  1:1201654/91‑1 (MQ=255)
                               gCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatattt                                                            <  1:2394388/91‑1 (MQ=255)
                                cAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatattta                                                           >  1:952955/1‑91 (MQ=255)
                                 aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                                          <  1:2990436/91‑1 (MQ=255)
                                 aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                                          <  1:62329/91‑1 (MQ=255)
                                 aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGGAGCAGATGCATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                                          <  1:64013/91‑1 (MQ=255)
                                  aCATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaacatttaat                                                         >  1:2385275/1‑91 (MQ=255)
                                   cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                                        >  1:1788311/1‑91 (MQ=255)
                                   cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                                        >  1:1227716/1‑91 (MQ=255)
                                   cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                                        >  1:2030444/1‑91 (MQ=255)
                                   cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                                        <  1:2932583/91‑1 (MQ=255)
                                    atatTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGa                                                       <  1:2236830/91‑1 (MQ=255)
                                    atatTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGa                                                       <  1:327735/91‑1 (MQ=255)
                                    atatTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGa                                                       <  1:657408/91‑1 (MQ=255)
                                      atTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAg                                                     >  1:44144/1‑91 (MQ=255)
                                      atTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAg                                                     >  1:43956/1‑91 (MQ=255)
                                         aaCATACAAACTTTTTGCTTAAAGAATAGTAGCCGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGAtt                                                  <  1:4061976/91‑1 (MQ=255)
                                          aCATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTa                                                 <  1:166410/91‑1 (MQ=255)
                                           cATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTaa                                                <  1:3124796/91‑1 (MQ=255)
                                            aTACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAt                                               <  1:3224191/91‑1 (MQ=255)
                                              aCAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtct                                             <  1:2235062/91‑1 (MQ=255)
                                               cAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtctt                                            <  1:171039/91‑1 (MQ=255)
                                               cAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtctt                                            <  1:2098897/91‑1 (MQ=255)
                                               cAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtctt                                            <  1:2098637/91‑1 (MQ=255)
                                               cAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtctt                                            <  1:1547278/91‑1 (MQ=255)
                                                  aCTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATcttctt                                         >  1:1072824/1‑91 (MQ=255)
                                                  aCTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATcttctt                                         >  1:3722534/1‑91 (MQ=255)
                                                   cTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCttt                                        >  1:1424257/1‑91 (MQ=255)
                                                    tttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTa                                       >  1:3032032/1‑91 (MQ=255)
                                                     ttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAg                                      <  1:263010/91‑1 (MQ=255)
                                                     ttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAg                                      <  1:2582117/91‑1 (MQ=255)
                                                     ttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAg                                      >  1:1592902/1‑91 (MQ=255)
                                                      tttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGt                                     <  1:1235450/91‑1 (MQ=255)
                                                        tGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTc                                   >  1:4054412/1‑91 (MQ=255)
                                                         gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                                  >  1:1269078/1‑91 (MQ=255)
                                                         gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                                  >  1:3251788/1‑91 (MQ=255)
                                                         gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                                  >  1:3251789/1‑91 (MQ=255)
                                                         gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                                  >  1:3841549/1‑91 (MQ=255)
                                                          cTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCt                                 <  1:3548826/91‑1 (MQ=255)
                                                          cTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCt                                 >  1:1673375/1‑91 (MQ=255)
                                                               aGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAAc                            >  1:2644547/1‑91 (MQ=255)
                                                                 aaTAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGa                          <  1:1224360/91‑1 (MQ=255)
                                                                 aaTAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGa                          <  1:578780/91‑1 (MQ=255)
                                                                   tagtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAAc                        >  1:1005087/1‑91 (MQ=255)
                                                                    agtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACt                       >  1:1386836/1‑91 (MQ=255)
                                                                      tagCAGATGCATAGCTTATGTATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACttt                     >  1:1012865/1‑91 (MQ=37)
                                                                      tagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACttt                     >  1:3698686/1‑91 (MQ=255)
                                                                      tagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACttt                     >  1:3699411/1‑91 (MQ=255)
                                                                       agCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTg                    >  1:1636278/1‑91 (MQ=255)
                                                                        gCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGc                   <  1:2043950/91‑1 (MQ=255)
                                                                        gCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGc                   >  1:2780262/1‑91 (MQ=255)
                                                                           gATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGtt                >  1:2923469/1‑91 (MQ=255)
                                                                             tACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAt              >  1:241481/1‑91 (MQ=255)
                                                                                aTAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAtctc           <  1:1019713/91‑1 (MQ=255)
                                                                                aTAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAtctc           <  1:1444034/91‑1 (MQ=255)
                                                                                aTAAGCTATGCATCTGCTACTCATTATTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAtctc           <  1:3851248/91‑1 (MQ=255)
                                                                                  aaGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCtt         <  1:1827750/91‑1 (MQ=255)
                                                                                  aaGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCtt         >  1:3564658/1‑91 (MQ=255)
                                                                                  aaGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCtt         >  1:3564504/1‑91 (MQ=255)
                                                                                   aGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTta        >  1:3622491/1‑91 (MQ=255)
                                                                                   aGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTta        >  1:3622492/1‑91 (MQ=255)
                                                                                    gcttATGTATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTta        >  1:1433450/4‑91 (MQ=255)
                                                                                    gcttATGTATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTta        >  1:811765/4‑91 (MQ=255)
                                                                                     cTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTtaat      >  1:1119339/1‑91 (MQ=255)
                                                                                     cTATGCATCTGCAACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTtaat      >  1:3176743/1‑91 (MQ=255)
                                                                                     cTATGCATCTGCAACTCATTTTTGTAATGTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTtaat      >  1:3176720/1‑91 (MQ=255)
                                                                                        aTGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTtaataa    >  1:672224/1‑91 (MQ=255)
                                                                                        aTGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTtaataa    >  1:674294/1‑91 (MQ=255)
                                                                                         tGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTTAATaaa   >  1:3616189/1‑91 (MQ=255)
                                                                                          gCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTTAATAAAt  >  1:2409845/1‑91 (MQ=255)
                                                                                           |                                                                                         
TTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTATCTCTTAATAAAT  >  USA300TCH1516_ALE/1776544‑1776723

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