Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 9,716 T→C 22.6% intergenic (+21/+66) gyrA → / ← nnrD DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE9,7160TC22.6% 198.0 / 58.1 102intergenic (+21/+66)gyrA/nnrDDNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base C (16/7);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 5.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 103 total aligned reads displayed.

TGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAT  >  USA300TCH1516_ALE/9629‑9802
                                                                                       |                                                                                      
tGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTacata                                                                                     <  1:1617092/91‑2 (MQ=255)
   aGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAg                                                                                  <  1:3096280/91‑1 (MQ=255)
    gTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgg                                                                                 <  1:1721741/91‑1 (MQ=38)
     tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAggg                                                                                >  1:1345945/1‑91 (MQ=255)
      aaGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAggga                                                                               <  1:4109226/91‑2 (MQ=38)
       aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                              <  1:3673558/91‑1 (MQ=255)
       aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                              <  1:2680181/91‑1 (MQ=255)
        gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAg                                                                             >  1:1576157/1‑91 (MQ=255)
         aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                            <  1:2347930/91‑1 (MQ=255)
         aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                            <  1:3614920/91‑1 (MQ=255)
         aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGt                                                                            <  1:1440429/91‑1 (MQ=37)
           aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                          <  1:1463194/91‑1 (MQ=37)
           aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                          <  1:2481503/91‑1 (MQ=37)
           aaGAGTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                          <  1:3457412/91‑1 (MQ=255)
           aaGACTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                          <  1:3457384/91‑1 (MQ=255)
            aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCtt                                                                         >  1:3230026/1‑91 (MQ=255)
            aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCtt                                                                         <  1:1762594/91‑1 (MQ=37)
             gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTa                                                                        >  1:2788708/1‑91 (MQ=255)
               tttCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtt                                                                      <  1:2574277/91‑1 (MQ=255)
               tttCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtt                                                                      <  1:415471/91‑1 (MQ=255)
                 tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtttt                                                                    >  1:158151/1‑91 (MQ=255)
                 tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtttt                                                                    >  1:3691692/1‑91 (MQ=255)
                 tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtttt                                                                    >  1:3109075/1‑91 (MQ=255)
                  cATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAttttt                                                                   >  1:3703005/1‑91 (MQ=255)
                    tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                 <  1:2972108/91‑1 (MQ=255)
                    tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                 <  1:3928548/91‑1 (MQ=255)
                    tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                 >  1:3489397/1‑91 (MQ=255)
                      gATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTAttttttatt                                                               <  1:1161996/91‑5 (MQ=21)
                       aTCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCt                                                              <  1:1071263/91‑1 (MQ=255)
                        tCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTa                                                             >  1:3122640/1‑91 (MQ=255)
                         cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAg                                                            <  1:3893458/91‑1 (MQ=255)
                         cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAg                                                            <  1:30969/91‑1 (MQ=255)
                          gtgtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGa                                                           <  1:3704679/91‑1 (MQ=255)
                            gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGaaa                                                         >  1:1714836/1‑91 (MQ=255)
                              tGAAGAAGATATACAACAATCATCAGATGAAGATTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGt                                                       <  1:3627566/91‑1 (MQ=255)
                              tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGt                                                       <  1:2168064/91‑1 (MQ=255)
                                 agaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAAt                                                    >  1:1155287/1‑91 (MQ=255)
                                 agaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAAt                                                    >  1:3594482/1‑91 (MQ=255)
                                  gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATg                                                   <  1:2555889/91‑1 (MQ=255)
                                  gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATg                                                   <  1:2555294/91‑1 (MQ=255)
                                   aagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGc                                                  <  1:3949599/91‑1 (MQ=255)
                                   aagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGc                                                  <  1:1744051/91‑1 (MQ=255)
                                    agaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCt                                                 >  1:1117795/1‑91 (MQ=255)
                                     gaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCtt                                                >  1:3425507/1‑91 (MQ=255)
                                     gaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCtt                                                <  1:2787432/91‑1 (MQ=255)
                                       tataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTa                                              <  1:672175/91‑1 (MQ=255)
                                        ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTAc                                             <  1:1592669/91‑1 (MQ=255)
                                          acaacaATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACta                                           >  1:2768821/1‑91 (MQ=255)
                                           caacaaTCATCAGATGAAGAATAATCAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACtat                                          <  1:3862185/91‑1 (MQ=255)
                                           caacaaTCATCAGATGAAGAATAATCAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACtat                                          <  1:3862310/91‑1 (MQ=255)
                                           caacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACtat                                          >  1:1853784/1‑91 (MQ=255)
                                              caaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATAtt                                       <  1:778971/91‑1 (MQ=255)
                                               aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTc                                      <  1:2077631/91‑1 (MQ=255)
                                                atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                     <  1:1527959/91‑1 (MQ=255)
                                                atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                     <  1:1526495/91‑1 (MQ=255)
                                                  catcaGATGAAGAATAATAAAAAATAAGACTTCTCTACATGTAGGGCAGTCTTATTTTTCTGCTAGAAAGTAATGCTTTACTATATTCAAt                                   >  1:3310705/1‑91 (MQ=255)
                                                  catcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAAt                                   >  1:2740485/1‑91 (MQ=255)
                                                  catcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAAt                                   <  1:1489927/91‑1 (MQ=255)
                                                   atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                  >  1:1777000/1‑91 (MQ=255)
                                                   atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                  <  1:2425218/91‑1 (MQ=255)
                                                    tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                                 >  1:3903292/1‑91 (MQ=255)
                                                    tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                                 >  1:553960/1‑91 (MQ=255)
                                                     caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                >  1:1688895/1‑91 (MQ=255)
                                                     caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                <  1:278057/91‑1 (MQ=255)
                                                        aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                             >  1:816957/1‑91 (MQ=255)
                                                        aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                             >  1:3372383/1‑91 (MQ=255)
                                                        aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                             >  1:3815369/1‑91 (MQ=255)
                                                         tGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGt                            <  1:3942108/91‑1 (MQ=255)
                                                                taataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACt                     <  1:2679918/91‑1 (MQ=255)
                                                                  ataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTaa                   >  1:3205210/1‑91 (MQ=255)
                                                                    aaaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACt                 <  1:2665341/91‑1 (MQ=255)
                                                                    aaaaaaTAAGACTTCCCTATATGAAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACt                 <  1:460606/91‑1 (MQ=255)
                                                                     aaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACtt                >  1:3030043/1‑91 (MQ=255)
                                                                     aaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACtt                <  1:2183122/91‑1 (MQ=255)
                                                                     aaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATCCAATGATTAGTAATGACTAACtt                >  1:2776475/1‑91 (MQ=255)
                                                                       aaaTAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTc              <  1:2025826/91‑1 (MQ=255)
                                                                       aaaTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTc              >  1:1400882/1‑91 (MQ=25)
                                                                        aaTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGCAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCt             >  1:3055574/1‑91 (MQ=25)
                                                                         aTAAGACTCCCCTACATATAGGGAAGTCTTCTTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa            >  1:1977390/1‑91 (MQ=37)
                                                                         aTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa            >  1:1976431/1‑91 (MQ=25)
                                                                         aTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa            >  1:1977509/1‑91 (MQ=25)
                                                                          tAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa           >  1:1138191/1‑91 (MQ=37)
                                                                          tAAGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa           >  1:2020293/1‑91 (MQ=25)
                                                                           aaGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt          <  1:860219/91‑1 (MQ=255)
                                                                           aaGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt          <  1:106499/91‑1 (MQ=255)
                                                                           aaGACGTCCCTAGATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt          <  1:1644515/91‑1 (MQ=255)
                                                                            aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:2052701/1‑91 (MQ=255)
                                                                            aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:2759346/1‑91 (MQ=255)
                                                                            aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:3833062/1‑91 (MQ=37)
                                                                            aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:1145817/1‑91 (MQ=37)
                                                                            aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:1145650/1‑91 (MQ=37)
                                                                            aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:2910571/1‑91 (MQ=37)
                                                                            aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt         >  1:81226/1‑91 (MQ=37)
                                                                             gACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg        >  1:3521989/1‑91 (MQ=255)
                                                                             gACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg        >  1:3521938/1‑91 (MQ=255)
                                                                             gACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg        >  1:2280421/1‑91 (MQ=255)
                                                                             gACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg        >  1:3824505/1‑91 (MQ=37)
                                                                                tccccTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt     >  1:1859187/3‑91 (MQ=38)
                                                                                   ccTATATGTAGGGGAGTCTTATTTTTATGCTATAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAt  <  1:3521153/91‑1 (MQ=255)
                                                                                   ccTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAt  >  1:206589/1‑91 (MQ=255)
                                                                                       |                                                                                      
TGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAT  >  USA300TCH1516_ALE/9629‑9802

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: