Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,776,622 A→G 5.1% intergenic (‑99/+52) 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,6220AG 5.1% 932.1 / 25.6 292intergenic (‑99/+52)clpX/tigATP‑dependent Clp protease ATP‑binding subunit ClpX/Trigger factor
Reads supporting (aligned to +/- strand):  ref base A (141/134);  new base G (8/7);  total (151/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 293 total aligned reads displayed.

CACCCCTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAT  >  USA300TCH1516_ALE/1776532‑1776711
                                                                                          |                                                                                          
cacCCCTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATa                                                                                            <  2:960405/91‑1 (MQ=255)
 acCCCTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATAc                                                                                           >  2:2418319/1‑91 (MQ=255)
  ccccTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACa                                                                                          >  1:1587703/1‑91 (MQ=255)
  ccccTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACa                                                                                          >  1:1986122/1‑91 (MQ=255)
  ccccTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACa                                                                                          <  1:1926103/91‑1 (MQ=255)
     cTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATaa                                                                                       >  2:2768393/1‑91 (MQ=255)
      tATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAg                                                                                      <  2:3148466/91‑1 (MQ=255)
        ttCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCt                                                                                   <  2:849846/91‑1 (MQ=255)
         tCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTa                                                                                  >  1:2463618/1‑91 (MQ=255)
         tCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTa                                                                                  <  1:2164506/91‑1 (MQ=255)
          cGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTAt                                                                                 >  1:330110/1‑91 (MQ=255)
          cGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTAt                                                                                 <  1:29850/91‑1 (MQ=255)
          cGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTAt                                                                                 >  1:2083712/1‑91 (MQ=255)
            ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGc                                                                               >  1:1139290/1‑91 (MQ=255)
            ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGc                                                                               >  2:1370428/1‑91 (MQ=255)
            ttAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGc                                                                               >  2:3344315/1‑91 (MQ=255)
             tAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTACGCa                                                                              >  2:1572431/1‑91 (MQ=255)
              aaaGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAACGAATAGTAGCAGATACATAAGCTATGCAt                                                                             >  2:1755248/1‑91 (MQ=255)
                  acgataCTAGATTGGATGTTACAATGCAACATATTAACCTACAAACTTTTTGCTTAAAGAATAGGAGCTGATGCATAGCTTATGTATCTGc                                                                         <  2:1832784/88‑1 (MQ=255)
                  aCTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGc                                                                         <  1:2581062/91‑1 (MQ=255)
                  aCTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGc                                                                         <  1:2919630/91‑1 (MQ=255)
                  aCTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGc                                                                         <  1:660758/91‑1 (MQ=255)
                  aCTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGc                                                                         <  2:2612023/91‑1 (MQ=255)
                    tataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTa                                                                       <  1:1382676/91‑1 (MQ=255)
                    tataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTa                                                                       <  2:814374/91‑1 (MQ=255)
                    tataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTa                                                                       >  1:513991/1‑91 (MQ=255)
                    tataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTa                                                                       <  1:3362212/91‑1 (MQ=255)
                     ataCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTAc                                                                      >  2:3258115/1‑91 (MQ=255)
                      taCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACt                                                                     >  1:2786340/1‑91 (MQ=255)
                       aCTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATGCATAGCTTATGTATCTGCTACta                                                                    <  2:1434262/91‑2 (MQ=37)
                        cTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCa                                                                   >  1:1790874/1‑91 (MQ=255)
                          aGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAtt                                                                 <  1:981800/91‑1 (MQ=255)
                            aTTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAtttt                                                               <  2:3957165/91‑1 (MQ=255)
                             ttGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAttttt                                                              <  2:338145/91‑1 (MQ=255)
                             ttGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCAttttt                                                              <  1:425176/91‑1 (MQ=255)
                               ggATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGt                                                            <  1:3184899/91‑1 (MQ=255)
                                 aTGTTACAATGCAACATATTAACATACAAACTTCTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTaa                                                          <  1:3501227/91‑1 (MQ=255)
                                  tGTTACAATGCAACATAGTAACATACAAACTTTTTGCTGAAAGAATAGTAGCAGATAGATAAGCTATGCAGCTGCTACTCATTTTTGTAat                                                         <  2:1253270/91‑1 (MQ=255)
                                   gTTACAACGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATATTTGTAatt                                                        <  1:1972309/91‑1 (MQ=255)
                                     tACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAattta                                                      >  2:1491442/1‑91 (MQ=255)
                                     gaCAACGAAACATATTAACATACAAACATTTGGCTTCAATAATAGTTGCAGATAGATAAGCTATGCATCTGCTACTCATTTTTGTAattta                                                      <  2:3640054/90‑1 (MQ=255)
                                        aaTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATATTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaata                                                   <  2:3478653/91‑1 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                  <  2:2666789/91‑1 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                  >  2:2929541/1‑91 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                  >  1:1469368/1‑91 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                  >  1:1002615/1‑91 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatat                                                  <  2:678187/91‑1 (MQ=255)
                                         aTGCAACATATTAACATACAAACTTTGTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTACTACTCATTTTTGTAatttaatat                                                  >  1:3621159/1‑91 (MQ=255)
                                          tGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatt                                                 >  1:872117/1‑91 (MQ=255)
                                             aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGTGGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                              <  2:395616/91‑1 (MQ=255)
                                             aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                              >  2:859523/1‑91 (MQ=255)
                                             aaCATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAatttaatatttaa                                              <  1:928308/91‑1 (MQ=255)
                                               cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                            <  2:1448649/91‑1 (MQ=255)
                                               cATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATg                                            >  1:684726/1‑91 (MQ=255)
                                                 tatTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGaa                                          >  2:3074330/1‑91 (MQ=255)
                                                 tatTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGaa                                          >  1:1054741/1‑91 (MQ=255)
                                                  atTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAg                                         >  2:3712620/1‑91 (MQ=255)
                                                   ttAACATACAAACTTTTTGCTTAAAGAATAGTGGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGa                                        <  2:2825067/91‑1 (MQ=255)
                                                     aaCATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGAtt                                      <  2:527107/91‑1 (MQ=255)
                                                      aCATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTa                                     >  1:3694395/1‑91 (MQ=255)
                                                      aCATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTa                                     >  1:3882870/1‑91 (MQ=255)
                                                       cATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTaa                                    <  2:2450652/91‑1 (MQ=255)
                                                        aTACAAACTTGTTGCTTAAAGAATAGTAGCTGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAt                                   <  2:3412008/91‑1 (MQ=255)
                                                          aCAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtct                                 <  2:1656046/91‑1 (MQ=255)
                                                          aCAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtct                                 <  1:1598270/91‑1 (MQ=255)
                                                          aCAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtct                                 <  1:3032043/91‑1 (MQ=255)
                                                          aCAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtct                                 <  1:2848458/91‑1 (MQ=255)
                                                            aaaCTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtcttc                               <  1:3652810/91‑1 (MQ=255)
                                                             aaCTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAAtcttct                              >  2:1928103/1‑91 (MQ=255)
                                                               cTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCTTTTTTGTAATTTAATATTTATTGAAGATTAATCTTCttt                            >  2:2578047/1‑91 (MQ=255)
                                                               cTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCttt                            <  1:3602819/91‑1 (MQ=255)
                                                               cTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATGTTTGTAATGTCATATTTAATGAAGATTAATCTTCttt                            >  1:1669257/1‑91 (MQ=255)
                                                                ttttttCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTa                           >  2:2306412/1‑91 (MQ=255)
                                                                tttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTa                           <  1:3461789/91‑1 (MQ=255)
                                                                 ttttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAg                          <  2:3753245/91‑1 (MQ=255)
                                                                  tttGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGt                         >  1:3164193/1‑91 (MQ=255)
                                                                   ttGCTTAAAGAATAGTAGCAGATACATAAGCTACGCAACTGCTACTCATTTTTGTAATTAAATATTTAATGAAGATTAATCTTCTTTAGtt                        >  1:3151991/1‑91 (MQ=255)
                                                                     gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                      >  1:2227703/1‑91 (MQ=255)
                                                                     gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                      <  1:2335401/91‑1 (MQ=255)
                                                                     gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                      >  2:721719/1‑91 (MQ=255)
                                                                     gCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTcc                      <  1:73644/91‑1 (MQ=255)
                                                                      cTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCt                     >  2:392606/1‑91 (MQ=255)
                                                                           aGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAAc                <  2:229082/91‑1 (MQ=255)
                                                                           aGAATAGTAGCAGATACATAAGCGATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAAc                <  2:573488/91‑1 (MQ=255)
                                                                             atgagtagCAGATGCATAGCTTATGTATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGa              >  1:913112/4‑91 (MQ=21)
                                                                              tgagtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATCTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGaa             >  2:16767/3‑91 (MQ=255)
                                                                              aTAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGaa             <  1:2645570/91‑1 (MQ=255)
                                                                                agtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACt           <  2:472543/91‑1 (MQ=255)
                                                                                agtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACt           >  1:1093629/1‑91 (MQ=255)
                                                                                 gtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACtt          >  2:2548299/1‑91 (MQ=255)
                                                                                 gtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACtt          >  2:648460/1‑91 (MQ=255)
                                                                                 gtagCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACtt          <  1:2080032/91‑1 (MQ=255)
                                                                                   agCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGACCTTTg        >  2:1187859/1‑91 (MQ=255)
                                                                                   agCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTg        >  1:386716/1‑91 (MQ=255)
                                                                                     cAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCg      >  1:268298/1‑91 (MQ=255)
                                                                                     cAGATACATAAGCTATGCATCTGCTACGCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCGTTAGTTCCTTCAACGAACTTTGCg      <  2:1320498/91‑1 (MQ=255)
                                                                                      aGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGt     >  1:1311750/1‑91 (MQ=255)
                                                                                       gATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGtt    <  1:2122379/91‑1 (MQ=255)
                                                                                        aTACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTa   <  2:1404816/91‑1 (MQ=255)
                                                                                        aTACATAAGCTATGCAGCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTa   <  2:2820008/91‑1 (MQ=255)
                                                                                          |                                                                                          
CACCCCTATTCGTTAAAGACTATACTAGATTGGATGTTACAATGCAACATATTAACATACAAACTTTTTGCTTAAAGAATAGTAGCAGATACATAAGCTATGCATCTGCTACTCATTTTTGTAATTTAATATTTAATGAAGATTAATCTTCTTTAGTTCCTTCAACGAACTTTGCGTTAT  >  USA300TCH1516_ALE/1776532‑1776711

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