Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,129,392 C→T 8.2% intergenic (+225/+13) hlb_2 → / ← USA300HOU_RS10870 Phospholipase C/putative leukocidin‑like protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,129,3920CT 8.2% 991.6 / 48.1 310intergenic (+225/+13)hlb_2/USA300HOU_RS10870Phospholipase C/putative leukocidin‑like protein 1
Reads supporting (aligned to +/- strand):  ref base C (141/143);  new base T (10/16);  total (151/159)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTGA  >  USA300TCH1516_ALE/2129302‑2129482
                                                                                          |                                                                                          
cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                                            >  2:3367377/1‑91 (MQ=255)
 gcTCAAACAAGATTCATACACTGCACGGCATAATTATTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                                           <  2:3906949/90‑1 (MQ=255)
 acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                                           <  1:305667/91‑1 (MQ=255)
   tCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                                         >  1:960839/1‑91 (MQ=255)
    cAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACt                                                                                        <  1:3496187/91‑1 (MQ=255)
    cAAACAAGATTCATACACTGCACGTCATAATAAAGCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACt                                                                                        <  2:1894211/91‑1 (MQ=18)
       aCAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACtttt                                                                                     <  1:1533558/91‑1 (MQ=18)
       aCAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACtttt                                                                                     >  2:3843796/1‑91 (MQ=255)
        cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTtta                                                                                    <  1:307732/91‑1 (MQ=17)
        cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTtta                                                                                    <  1:1946037/91‑1 (MQ=18)
        cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTtta                                                                                    >  2:1254570/1‑91 (MQ=255)
         aaGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttat                                                                                   >  2:2809317/1‑91 (MQ=255)
          aGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttatt                                                                                  >  2:116503/1‑91 (MQ=255)
            aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttattta                                                                                <  1:105076/91‑3 (MQ=17)
            aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTAc                                                                                >  1:2126433/1‑91 (MQ=255)
            aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTAc                                                                                >  2:3132342/1‑91 (MQ=255)
             ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                                               <  1:450610/91‑1 (MQ=255)
              gcATACGCTGCACGTCGTAATATATCTATCTATTGAAATAGAATTAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTGCtt                                                                              <  2:1708713/90‑1 (MQ=255)
               cATACACTGCCCGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                             <  2:223856/91‑1 (MQ=255)
               cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                             <  2:2080087/91‑1 (MQ=255)
               cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                             >  2:2162688/1‑91 (MQ=255)
               cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                             >  1:2845216/1‑91 (MQ=255)
                 tACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAtt                                                                           >  2:751028/1‑91 (MQ=255)
                 tACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAtt                                                                           >  1:3107725/1‑91 (MQ=255)
                  acacTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAttt                                                                          >  2:2492695/1‑91 (MQ=255)
                  acacTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAttt                                                                          <  2:1112715/91‑1 (MQ=255)
                     cTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCtt                                                                       >  1:234957/1‑91 (MQ=255)
                      tGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCttt                                                                      <  1:93958/91‑1 (MQ=255)
                      tGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCttt                                                                      >  2:3497249/1‑91 (MQ=255)
                      tGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCttt                                                                      >  1:1390071/1‑91 (MQ=255)
                        cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                                                    <  2:655445/91‑1 (MQ=255)
                        cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                                                    <  2:3595697/91‑1 (MQ=255)
                          cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                                                  <  1:80596/91‑1 (MQ=255)
                          cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                                                  >  1:2586797/1‑91 (MQ=255)
                          cGTCATAATAAATATATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                                                  <  1:548778/91‑1 (MQ=255)
                            tCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcatta                                                                >  2:882173/1‑91 (MQ=255)
                               taataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatca                                                             <  2:1394246/91‑1 (MQ=255)
                               taataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatca                                                             >  1:1466094/1‑91 (MQ=255)
                                aataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcat                                                            <  2:2364810/91‑1 (MQ=255)
                                aataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcat                                                            >  1:3872239/1‑91 (MQ=255)
                                aataaACCTATCTATTCAAATATACATAAAAGTTACCTACTACAATCTATGTACCAGGCAACTTTTATTACTTATTTCTTTtcattatcat                                                            >  2:2559474/1‑91 (MQ=255)
                                 ataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatt                                                           >  2:2838577/1‑91 (MQ=255)
                                  taaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatta                                                          <  1:1362140/91‑1 (MQ=255)
                                   aaaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTaa                                                         <  1:489062/91‑1 (MQ=255)
                                   aaaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTaa                                                         <  2:2916428/91‑1 (MQ=255)
                                   aaaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTaa                                                         >  1:1879843/1‑91 (MQ=255)
                                    aaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAg                                                        <  2:3617626/91‑1 (MQ=255)
                                     atctatctATTCTAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACCTTTATTACTTATTTCTTTTCATTATCATTAAGt                                                       <  2:3292519/91‑1 (MQ=255)
                                     atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                                                       <  2:1677655/91‑1 (MQ=255)
                                     atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                                                       <  2:3409610/91‑1 (MQ=255)
                                     atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                                                       <  1:2611178/91‑1 (MQ=255)
                                        tatctatTCAAATATAAATAAAAGTTACCTACTACATTTTCTGTAACCAGCAACAGTTATTAATTATCTTTTTTCATTATCATTAAGTACt                                                    >  2:3675373/1‑91 (MQ=255)
                                        tatctatTCAAATATAAATAAAAGTGACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACt                                                    <  1:3676943/91‑1 (MQ=255)
                                         atctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtt                                                   >  1:1918300/1‑91 (MQ=255)
                                          tctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACttt                                                  >  2:547590/1‑91 (MQ=255)
                                           ctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtttt                                                 <  1:3703108/91‑1 (MQ=255)
                                           ctatTCAAATATAAAAGAACGTTACCTACTACGTGCTATGTAGCAGGCCACTTTTATTACTTATTTCTTTTCATTATCAGTAAGTACtttt                                                 <  2:1100531/91‑1 (MQ=255)
                                            tatTCAAATATAAATAAAAGTTACATACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTa                                                <  2:3685667/91‑1 (MQ=255)
                                             atTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTAc                                               <  2:3486883/91‑1 (MQ=255)
                                                cAAATATAAATAAAAGTTACCTACTACATTCCCTGTAGCAGGCAACTTTTATTACTTATTTCTTTTCGCTATCATTAAGCACTTTTACaaa                                            >  2:3022224/1‑91 (MQ=255)
                                                 aaaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAAc                                           <  1:2073622/91‑1 (MQ=255)
                                                   atatAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACtt                                         >  1:3602443/1‑91 (MQ=255)
                                                    tataAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTc                                        <  1:881884/91‑1 (MQ=255)
                                                     ataaataaAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTca                                       <  2:242538/91‑1 (MQ=255)
                                                      taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                                      <  1:1729576/91‑1 (MQ=255)
                                                      taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                                      <  1:711358/91‑1 (MQ=255)
                                                        aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt                                    >  2:3539717/1‑91 (MQ=255)
                                                         ataaaAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt                                   >  2:1294717/1‑91 (MQ=17)
                                                          taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa                                  <  1:2985310/91‑1 (MQ=255)
                                                          taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa                                  <  1:2874560/91‑1 (MQ=255)
                                                            aaaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtg                                >  1:1813722/1‑91 (MQ=255)
                                                             aaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgt                               <  1:3448263/91‑1 (MQ=255)
                                                             aaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgt                               <  1:3348634/91‑1 (MQ=255)
                                                             aaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgt                               >  1:1648912/1‑91 (MQ=255)
                                                             aaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgt                               >  2:1767609/1‑91 (MQ=255)
                                                                 tACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCt                           <  2:151243/91‑1 (MQ=255)
                                                                 tACCTACTACATTCTATGTAGCAGGCAACTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt                          <  2:1407892/91‑1 (MQ=255)
                                                                  aCCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt                          >  1:2412462/1‑91 (MQ=255)
                                                                  aCCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt                          >  1:2892634/1‑91 (MQ=255)
                                                                   ccTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTc                         <  2:3274362/91‑1 (MQ=255)
                                                                     tactacATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCa                       <  1:2161272/91‑1 (MQ=255)
                                                                      actacATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCaa                      <  1:3828556/91‑1 (MQ=255)
                                                                             agaaTGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACtt               >  1:3193431/4‑91 (MQ=39)
                                                                              cTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTc              >  1:1010619/1‑91 (MQ=255)
                                                                              cTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTc              <  1:821002/91‑1 (MQ=255)
                                                                               tATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCa             >  1:3775928/1‑91 (MQ=255)
                                                                               tATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCa             >  2:455773/1‑91 (MQ=255)
                                                                                aTGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCat            >  1:1802150/1‑91 (MQ=255)
                                                                                 tGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCata           <  2:921633/91‑1 (MQ=255)
                                                                                 tGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCata           >  2:1317997/1‑91 (MQ=255)
                                                                                  gtagtaGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCatat          >  1:3574082/1‑91 (MQ=255)
                                                                                   tAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCAtata         >  1:2219627/1‑91 (MQ=255)
                                                                                     gCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAAt       <  2:2723412/91‑1 (MQ=255)
                                                                                       aGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGc     <  1:2743857/91‑1 (MQ=255)
                                                                                        ggCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCt    >  1:674044/1‑91 (MQ=255)
                                                                                        ggCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCt    <  2:1569375/91‑1 (MQ=255)
                                                                                         gCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTg   <  2:3038923/91‑1 (MQ=255)
                                                                                         gCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTg   <  1:1262409/91‑1 (MQ=255)
                                                                                          cAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTGa  <  1:58901/91‑1 (MQ=255)
                                                                                          cAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTGa  >  1:3271968/1‑91 (MQ=255)
                                                                                          |                                                                                          
CACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTCATATAATGCTGA  >  USA300TCH1516_ALE/2129302‑2129482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: