Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 9,719 G→A 21.4% intergenic (+24/+63) 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,7190GA21.4% 220.1 / 55.5 103intergenic (+24/+63)gyrA/nnrDDNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (37/44);  new base A (15/7);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 103 total aligned reads displayed.

AGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTA  >  USA300TCH1516_ALE/9632‑9808
                                                                                       |                                                                                         
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:2680181/91‑1 (MQ=255)
    aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                                    <  1:3673558/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:3109075/1‑91 (MQ=255)
              tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtttt                                                                          >  1:158151/1‑91 (MQ=255)
              tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAtttt                                                                          >  1:3691692/1‑91 (MQ=255)
               cATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTAttttt                                                                         >  1:3703005/1‑91 (MQ=255)
                 tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                       >  1:3489397/1‑91 (MQ=255)
                 tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                       <  1:2972108/91‑1 (MQ=255)
                 tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTAt                                                                       <  1:3928548/91‑1 (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:30969/91‑1 (MQ=255)
                      cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAg                                                                  <  1:3893458/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:3594482/1‑91 (MQ=255)
                              agaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAAt                                                          >  1:1155287/1‑91 (MQ=255)
                               gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATg                                                         <  1:2555294/91‑1 (MQ=255)
                               gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATg                                                         <  1:2555889/91‑1 (MQ=255)
                                aagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGc                                                        <  1:1744051/91‑1 (MQ=255)
                                aagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGc                                                        <  1:3949599/91‑1 (MQ=255)
                                 agaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCt                                                       >  1:1117795/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:1526495/91‑1 (MQ=255)
                                             atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                           <  1:1527959/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:553960/1‑91 (MQ=255)
                                                 tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                                       >  1:3903292/1‑91 (MQ=255)
                                                  caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                      >  1:1688895/1‑91 (MQ=255)
                                                  caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                      <  1:278057/91‑1 (MQ=255)
                                                     aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                   >  1:3815369/1‑91 (MQ=255)
                                                     aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                   >  1:3372383/1‑91 (MQ=255)
                                                     aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                   >  1:816957/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:2054344/1‑91 (MQ=255)
                                                                         aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt               >  1:2052701/1‑91 (MQ=255)
                                                                         aGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt               >  1:2759346/1‑91 (MQ=255)
                                                                         aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt               >  1:1145817/1‑91 (MQ=37)
                                                                         aGACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt               >  1:3833062/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:3521938/1‑91 (MQ=255)
                                                                          gACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg              >  1:2280421/1‑91 (MQ=255)
                                                                          gACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg              >  1:3521989/1‑91 (MQ=255)
                                                                          gACTCCCCTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg              >  1:3824505/1‑91 (MQ=37)
                                                                             tccccTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt           >  1:1859187/3‑91 (MQ=38)
                                                                             tccccTACATATAGGGAAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt           >  1:3408162/3‑91 (MQ=38)
                                                                                ccTATATGTAGGGGAGTCTTATTTTTATGCTATAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAt        <  1:3521153/91‑1 (MQ=255)
                                                                                      tGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTa  >  1:912616/1‑91 (MQ=255)
                                                                                       |                                                                                         
AGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTATTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTA  >  USA300TCH1516_ALE/9632‑9808

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: