Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,582,474 G→A 8.6% 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 8.6% 671.0 / 37.7 220H107H (CAC→CATclpP1ATP‑dependent Clp protease proteolytic subunit 1
Reads supporting (aligned to +/- strand):  ref base G (105/96);  new base A (9/10);  total (114/106)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 228 total aligned reads displayed.

TTTTTCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGC  >  USA300TCH1516_ALE/1582384‑1582564
                                                                                          |                                                                                          
  tttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTtgtg                                                                                          >  2:209432/1‑91 (MQ=255)
   ttCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGc                                                                                         <  2:2258426/91‑1 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                      <  2:1812572/91‑1 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATa                                                                                      <  2:67784/91‑1 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCaca                                                                                      <  2:3438592/91‑3 (MQ=255)
      aaGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCaca                                                                                      <  2:2775216/91‑3 (MQ=255)
       aGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaa                                                                                     <  1:443725/91‑1 (MQ=255)
       aGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaa                                                                                     <  1:2032687/91‑1 (MQ=255)
       aGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaa                                                                                     >  1:2102193/1‑91 (MQ=255)
       aGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATaa                                                                                     <  1:2183874/91‑1 (MQ=255)
            aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATa                                                                                >  1:466717/1‑91 (MQ=255)
            aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                                <  1:925806/91‑1 (MQ=255)
            aTCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATa                                                                                <  1:1501750/91‑1 (MQ=255)
              ccGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgt                                                                              >  1:3861424/1‑91 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTGATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                             <  2:3057876/91‑1 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                             >  1:2702569/1‑91 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAgtg                                                                             <  1:407785/91‑1 (MQ=255)
               cGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgtt                                                                             <  1:2045486/91‑2 (MQ=255)
                gCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTATGCACAAAAATAgttt                                                                            <  1:2845348/91‑3 (MQ=39)
                 cTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTc                                                                           >  2:2717679/1‑91 (MQ=255)
                   gTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcac                                                                         <  1:1035234/91‑1 (MQ=255)
                      tttcttAATTCTTCTGCATTACCTTCAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcaccac                                                                      <  2:2661243/91‑1 (MQ=255)
                      tttcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcaccac                                                                      <  2:2073452/91‑1 (MQ=255)
                        tcttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTcac                                                                         <  2:696220/86‑1 (MQ=255)
                          ttAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCAt                                                                  <  1:733007/91‑1 (MQ=255)
                           tAATTCTGCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATa                                                                 <  1:1952110/91‑1 (MQ=255)
                             aTTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGc                                                               >  1:1301365/1‑91 (MQ=255)
                              ttcttcTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCg                                                              >  2:3704763/1‑91 (MQ=255)
                                 ttctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATa                                                           >  2:1021296/1‑91 (MQ=255)
                                   ctGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAac                                                         >  2:1632557/1‑91 (MQ=255)
                                     gCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAacac                                                       >  2:1996576/1‑91 (MQ=255)
                                         tACCTACAGTCATAACCCATGAATTATGAACCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGc                                                   <  2:1847837/91‑1 (MQ=255)
                                          aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                  >  1:2303175/1‑91 (MQ=255)
                                          aCCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCa                                                  <  2:254688/91‑1 (MQ=255)
                                           ccTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCaa                                                 >  1:3412876/1‑91 (MQ=255)
                                            cTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAAt                                                <  2:1394420/91‑1 (MQ=255)
                                             tACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAAtt                                               >  2:73242/1‑91 (MQ=255)
                                               caatCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                             <  2:811956/88‑1 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                             <  2:565758/91‑1 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                             <  1:83775/91‑1 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                             >  2:446114/1‑91 (MQ=255)
                                               cAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGa                                             <  1:742403/91‑1 (MQ=255)
                                                 gTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATg                                           <  2:3837169/91‑1 (MQ=255)
                                                    aTAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGcct                                        <  2:498058/91‑2 (MQ=255)
                                                    aTAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCg                                        <  1:789223/91‑1 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGcctt                                       <  1:1265412/91‑3 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGcctt                                       <  1:2218990/91‑3 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGc                                       >  2:3569673/1‑91 (MQ=255)
                                                     tAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGc                                       >  2:127501/1‑91 (MQ=255)
                                                       aCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGccttag                                     <  2:1774780/91‑5 (MQ=255)
                                                       aCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGccttag                                     <  2:3775440/91‑5 (MQ=255)
                                                       aCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTa                                     <  2:2907529/91‑1 (MQ=255)
                                                       aCCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTa                                     <  2:934157/91‑1 (MQ=255)
                                                        cccATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTaa                                    >  1:3673167/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                 >  1:2366722/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                 >  1:2357639/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                 >  2:2513437/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCACAGCGATAACACTAGCAATTGATGCCGCTAAGGc                                 >  1:3450170/1‑91 (MQ=255)
                                                           aTGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACCCATAGCGATAACACTAGCAATTGATGCCGCTAAggg                                 >  2:3848815/1‑90 (MQ=255)
                                                             gAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCAt                               <  2:2605993/91‑1 (MQ=255)
                                                              aaTTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATc                              >  1:3171492/1‑91 (MQ=255)
                                                               aTTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCg                             >  1:580898/1‑91 (MQ=255)
                                                                 tATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGAc                           >  1:3260354/1‑91 (MQ=255)
                                                                  aTGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACa                          <  1:2706505/91‑1 (MQ=255)
                                                                  aTGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACa                          >  1:2440762/1‑91 (MQ=255)
                                                                  aTGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACa                          >  2:3212022/1‑91 (MQ=255)
                                                                  aTGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACa                          >  1:3766804/1‑91 (MQ=255)
                                                                   tGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACAt                         <  2:3067854/91‑1 (MQ=255)
                                                                   tGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACAt                         >  1:1848275/1‑91 (MQ=255)
                                                                   tGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATGGATGCCGCTAAGGCATCGACAt                         >  2:2950499/1‑91 (MQ=255)
                                                                    gAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATa                        >  2:2339000/1‑91 (MQ=255)
                                                                    gAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATa                        >  1:2931962/1‑91 (MQ=255)
                                                                    gAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATa                        >  2:3300233/1‑91 (MQ=255)
                                                                    gAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATa                        >  2:1950321/1‑91 (MQ=255)
                                                                         aTTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGatat                   <  2:710675/91‑1 (MQ=255)
                                                                           tAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTa                 >  1:2098224/1‑91 (MQ=255)
                                                                              aaaaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAAtt              <  1:3337696/91‑1 (MQ=255)
                                                                                aaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAAtttt            <  1:3479691/91‑1 (MQ=255)
                                                                                aaCTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAAtttt            >  1:3828529/1‑91 (MQ=255)
                                                                                  cTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTg          >  2:1495098/1‑91 (MQ=255)
                                                                                    aTTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCGATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCa        >  1:2969336/1‑91 (MQ=255)
                                                                                    aTTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCa        >  2:2371125/1‑91 (MQ=255)
                                                                                     tttttATGCACAAAAATAGTGTCACCACTCATAGCGAAAACACTAGCAATTTATGCCGCTAAGGCGTCGACATAGATATTAATTTCTGCAg       >  1:3615713/1‑91 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      >  1:2244596/1‑91 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      <  2:3668510/91‑1 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      >  2:1490030/1‑91 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      <  2:3439644/91‑1 (MQ=255)
                                                                                      ttttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      >  2:2039988/1‑91 (MQ=255)
                                                                                      ttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      >  1:2205042/1‑91 (MQ=255)
                                                                                      ttttATGCACAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAgg      >  1:1109635/1‑91 (MQ=255)
                                                                                       tttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGa     >  2:1400007/1‑91 (MQ=255)
                                                                                        ttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt    <  2:3260354/91‑1 (MQ=255)
                                                                                        ttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt    <  2:3171492/91‑1 (MQ=255)
                                                                                        ttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt    <  1:3854874/91‑1 (MQ=255)
                                                                                        ttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt    <  1:1810443/91‑1 (MQ=255)
                                                                                        ttGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGAt    <  1:1765517/91‑1 (MQ=255)
                                                                                         tgtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATg   <  2:31098/91‑1 (MQ=255)
                                                                                          gtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc  <  2:2638808/91‑1 (MQ=255)
                                                                                          gtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc  <  1:2077307/91‑1 (MQ=255)
                                                                                          gtgCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGc  <  2:1766677/91‑1 (MQ=255)
                                                                                          |                                                                                          
TTTTTCAAGTAAATCCGCTGTCTTTCTTAATTCTTCTGCATTACCTACAGTCATAACCCATGAATTATGAATCATTAAAAAACTATTTTTGTGCATAAAAATAGTGTCACCACTCATAGCGATAACACTAGCAATTGATGCCGCTAAGGCATCGACATAGATATTAATTTTTGCAGGATGC  >  USA300TCH1516_ALE/1582384‑1582564

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: