Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,582,479 T→C 15.1% M106V (ATG→GTG)  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,4790TC15.1% 289.3 / 37.1 112M106V (ATG→GTG) clpP1ATP‑dependent Clp protease proteolytic subunit 1
Reads supporting (aligned to +/- strand):  ref base T (46/49);  new base C (10/7);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGCATTT  >  USA300TCH1516_ALE/1582389‑1582568
                                                                                          |                                                                                         
cAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCAt                                                                                           >  1:1058257/1‑91 (MQ=255)
 aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                          >  1:322434/1‑91 (MQ=255)
 aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                          >  1:323009/1‑91 (MQ=255)
     aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                      >  1:2604601/1‑91 (MQ=255)
     aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                      >  1:612988/1‑91 (MQ=255)
     aaaTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaaaaa                                                                                      <  1:3264298/91‑1 (MQ=255)
       aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATa                                                                                    >  1:581508/1‑91 (MQ=255)
       aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                                    <  1:2977477/91‑1 (MQ=255)
       aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                                    <  1:3125643/91‑1 (MQ=255)
        tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAg                                                                                   <  1:1147393/91‑1 (MQ=255)
        tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAg                                                                                   >  1:2933570/1‑91 (MQ=255)
        tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAg                                                                                   <  1:2524638/91‑1 (MQ=255)
        tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAg                                                                                   <  1:1603838/91‑1 (MQ=255)
        tCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAg                                                                                   <  1:1603483/91‑1 (MQ=255)
          cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                                 <  1:4067251/91‑1 (MQ=255)
          cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                                 >  1:3313293/1‑91 (MQ=255)
          cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                                 <  1:4065116/91‑1 (MQ=255)
          cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgtt                                                                                 <  1:1222141/91‑2 (MQ=255)
          cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgtt                                                                                 <  1:1221628/91‑2 (MQ=255)
            cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                               >  1:3051723/1‑91 (MQ=255)
            cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                               >  1:2793700/1‑91 (MQ=255)
             tGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTca                                                                              <  1:1271906/91‑1 (MQ=255)
              gTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcac                                                                             >  1:353354/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)
                    cttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCa                                                                       <  1:923371/91‑1 (MQ=255)
                       aaTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCATTCATAg                                                                    <  1:3226948/91‑1 (MQ=255)
                       aaTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAg                                                                    <  1:122090/91‑1 (MQ=255)
                        aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                                   <  1:1240675/91‑1 (MQ=255)
                        aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                                   <  1:1900374/91‑1 (MQ=255)
                        aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                                   <  1:2231554/91‑1 (MQ=255)
                        aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                                   <  1:2231595/91‑1 (MQ=255)
                          tcttctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGa                                                                 >  1:3483931/1‑91 (MQ=255)
                          tcttctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGa                                                                 >  1:1567852/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)
                                     aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                      <  1:3094060/91‑1 (MQ=255)
                                      ccTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCaa                                                     <  1:2058502/91‑1 (MQ=255)
                                      ccTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCaa                                                     >  1:489965/1‑91 (MQ=255)
                                         aCAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTg                                                  >  1:2277563/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:2788088/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)
                                                     cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                      <  1:1808214/91‑1 (MQ=255)
                                                     cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                      >  1:527860/1‑91 (MQ=255)
                                                     cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                      >  1:3613706/1‑91 (MQ=255)
                                                     cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                      <  1:1807907/91‑1 (MQ=255)
                                                     cATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAgg                                      >  1:530187/1‑91 (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)
                                                                                tttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGTAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAg           >  1:3301980/1‑91 (MQ=255)
                                                                                tttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAg           >  1:2189105/1‑91 (MQ=255)
                                                                                 ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg          >  1:2928843/1‑91 (MQ=255)
                                                                                 ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg          >  1:2930334/1‑91 (MQ=255)
                                                                                  tttatgCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGa         >  1:1130792/5‑91 (MQ=255)
                                                                                  tttGTGCATAAAAATAGTGTCACCACTCATAGGGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATCTTAATTTTTGCAGGa         <  1:3541761/91‑1 (MQ=255)
                                                                                  tttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGa         <  1:3894615/91‑1 (MQ=255)
                                                                                   ttatgCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATCTTAATTTTTGCAGGAt        >  1:273818/4‑91 (MQ=255)
                                                                                   ttatgCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt        >  1:829560/4‑91 (MQ=255)
                                                                                    tgtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGCTATTAATTTTTGCAGGATg       >  1:266799/1‑91 (MQ=255)
                                                                                    tgtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATg       <  1:3355463/91‑1 (MQ=255)
                                                                                     atgCACAAAAATAGTGTCACCACTCATAGCGATAACAGTAGCAATTCATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc      >  1:3302890/2‑91 (MQ=255)
                                                                                     atgCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc      >  1:1449567/2‑91 (MQ=255)
                                                                                     atgCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTCATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc      >  1:3302200/2‑91 (MQ=255)
                                                                                       gCACAAAAATAGCGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGCAt    >  1:3196000/1‑91 (MQ=255)
                                                                                        cATAAAAATCGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGCAtt   <  1:3972513/91‑1 (MQ=255)
                                                                                         aTAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGCAttt  >  1:1698177/1‑91 (MQ=255)
                                                                                          |                                                                                         
CAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGCATTT  >  USA300TCH1516_ALE/1582389‑1582568

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: