Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,582,474 G→A 9.3% H107H (CAC→CAT clpP1 ← ATP‑dependent Clp protease proteolytic subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,582,4740GA 9.3% 330.7 / 18.4 107H107H (CAC→CATclpP1ATP‑dependent Clp protease proteolytic subunit 1
Reads supporting (aligned to +/- strand):  ref base G (46/51);  new base A (3/7);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

TTTTTCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATG  >  USA300TCH1516_ALE/1582384‑1582563
                                                                                          |                                                                                         
tttttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTtg                                                                                           <  1:525560/91‑1 (MQ=255)
  tttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTtgtg                                                                                         <  1:1804872/91‑1 (MQ=255)
   ttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGc                                                                                        <  1:3380807/91‑1 (MQ=255)
   ttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGc                                                                                        <  1:3897688/91‑1 (MQ=255)
    tCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCa                                                                                       <  1:3355085/91‑1 (MQ=255)
    tCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCa                                                                                       >  1:3631980/1‑91 (MQ=255)
     cAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCAt                                                                                      >  1:1058257/1‑91 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                     >  1:323009/1‑91 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                     >  1:322434/1‑91 (MQ=255)
       aGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaa                                                                                    <  1:3584278/91‑1 (MQ=255)
          aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                 >  1:612988/1‑91 (MQ=255)
          aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                 <  1:3264298/91‑1 (MQ=255)
          aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                 >  1:2604601/1‑91 (MQ=255)
            aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                               <  1:3125643/91‑1 (MQ=255)
            aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                               <  1:2977477/91‑1 (MQ=255)
             tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAg                                                                              <  1:1147393/91‑1 (MQ=255)
             tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAg                                                                              >  1:2933570/1‑91 (MQ=255)
             tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAg                                                                              <  1:1603483/91‑1 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                            <  1:4067251/91‑1 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                            <  1:4065116/91‑1 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                            >  1:3313738/1‑91 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                            >  1:3313293/1‑91 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgtt                                                                            <  1:1221628/91‑2 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgtt                                                                            <  1:1222141/91‑2 (MQ=255)
                 cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                          >  1:2750086/1‑91 (MQ=255)
                 cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                          >  1:2793700/1‑91 (MQ=255)
                 cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                          >  1:3051723/1‑91 (MQ=255)
                   gTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcac                                                                        >  1:622506/1‑91 (MQ=255)
                    tctttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcacc                                                                       <  1:2634062/91‑1 (MQ=255)
                     ctttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAACAAACTATTTTTGTGCATAAAAATAGTGTcacca                                                                      >  1:3401175/1‑91 (MQ=255)
                     ctttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcacca                                                                      >  1:535361/1‑91 (MQ=255)
                     ctttattAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcacca                                                                      <  1:3458365/91‑1 (MQ=255)
                      tttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcaccac                                                                     <  1:26879/91‑1 (MQ=255)
                      tttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcaccac                                                                     <  1:26404/91‑1 (MQ=255)
                      tttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcaccac                                                                     >  1:969689/1‑91 (MQ=255)
                       ttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACt                                                                    <  1:1223006/91‑1 (MQ=255)
                         cttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCa                                                                  <  1:923371/91‑1 (MQ=255)
                            aaTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCATTCATAg                                                               <  1:3226948/91‑1 (MQ=255)
                            aaTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAg                                                               <  1:122090/91‑1 (MQ=255)
                             aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                              <  1:2231595/91‑1 (MQ=255)
                             aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                              <  1:2231554/91‑1 (MQ=255)
                             aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                              <  1:1240675/91‑1 (MQ=255)
                               tcttctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGa                                                            >  1:1567852/1‑91 (MQ=255)
                               tcttctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGa                                                            >  1:3483931/1‑91 (MQ=255)
                                    tGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAaca                                                       <  1:1267397/91‑1 (MQ=255)
                                          aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                 <  1:2359319/91‑1 (MQ=255)
                                          aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                 <  1:3623887/91‑1 (MQ=255)
                                          aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                 >  1:2711131/1‑91 (MQ=255)
                                           ccTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCaa                                                <  1:2058502/91‑1 (MQ=255)
                                           ccTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCaa                                                >  1:489965/1‑91 (MQ=255)
                                             tACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAAtt                                              <  1:1752869/91‑1 (MQ=255)
                                              aCAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTg                                             >  1:2277563/1‑91 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCGTAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                            >  1:4146150/1‑91 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                            <  1:1303599/91‑1 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                            <  1:273904/91‑1 (MQ=255)
                                                aGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGAt                                           >  1:1185755/1‑91 (MQ=255)
                                                    aTAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGcct                                       <  1:3167169/91‑2 (MQ=255)
                                                    aTAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCg                                       <  1:3267356/91‑1 (MQ=255)
                                                    aTAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCg                                       <  1:1897239/91‑1 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGc                                      >  1:603906/1‑91 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGc                                      >  1:604280/1‑91 (MQ=255)
                                                      aaCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCt                                     >  1:923942/1‑91 (MQ=255)
                                                      aaCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCt                                     >  1:924366/1‑91 (MQ=255)
                                                       aCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGccttag                                    <  1:254907/91‑5 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 <  1:1808214/91‑1 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 >  1:3613706/1‑91 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 <  1:3225197/91‑1 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 >  1:527860/1‑91 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 >  1:530187/1‑91 (MQ=255)
                                                          cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                 <  1:1807907/91‑1 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                >  1:3093713/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                >  1:1561285/1‑91 (MQ=255)
                                                            tGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCa                               <  1:3714307/91‑1 (MQ=255)
                                                               aTTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCg                            >  1:456707/1‑91 (MQ=255)
                                                               aTTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCg                            <  1:214119/91‑1 (MQ=255)
                                                                ttATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGa                           <  1:162474/91‑1 (MQ=255)
                                                                ttATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGa                           <  1:162849/91‑1 (MQ=255)
                                                                  aTGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACa                         >  1:2806495/1‑91 (MQ=255)
                                                                   tGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACAt                        >  1:3138859/1‑91 (MQ=255)
                                                                    gAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATa                       <  1:3454418/91‑1 (MQ=255)
                                                                      aTCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGa                     >  1:1522111/1‑91 (MQ=255)
                                                                       tCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGat                    >  1:368282/1‑91 (MQ=255)
                                                                         aTTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGatat                  <  1:1317218/91‑1 (MQ=255)
                                                                           tAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTa                <  1:1310411/91‑1 (MQ=255)
                                                                           tAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTa                <  1:1187486/91‑1 (MQ=255)
                                                                            aaaaaaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTaa               <  1:1836825/91‑1 (MQ=255)
                                                                            aaaaaaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTaa               <  1:1837111/91‑1 (MQ=255)
                                                                                aaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATCACGCTAGCAATTGATGCCGCTAAGGGATCGACATAGATATTAAtttt           <  1:2340294/91‑1 (MQ=255)
                                                                                aaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGTCGCTAAGGCATCGACATAGATATTAAtttt           >  1:1557456/1‑91 (MQ=255)
                                                                                  cTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTg         >  1:3014744/1‑91 (MQ=255)
                                                                                    aTTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCa       >  1:1549748/1‑91 (MQ=255)
                                                                                     tttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAg      >  1:2189105/1‑91 (MQ=255)
                                                                                     tttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAg      >  1:611471/1‑91 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg     >  1:2928843/1‑91 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg     >  1:2930334/1‑91 (MQ=255)
                                                                                      ttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTATTGCAgg     >  1:2064713/1‑91 (MQ=255)
                                                                                       tttGTGCATAAAAATAGTGTCACCACTCATAGGGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATCTTAATTTTTGCAGGa    <  1:3541761/91‑1 (MQ=255)
                                                                                       tttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGa    <  1:3894615/91‑1 (MQ=255)
                                                                                         tgtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGCTATTAATTTTTGCAGGATg  >  1:266799/1‑91 (MQ=255)
                                                                                         tgtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATg  <  1:3355463/91‑1 (MQ=255)
                                                                                          |                                                                                         
TTTTTCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATG  >  USA300TCH1516_ALE/1582384‑1582563

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: