Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 447,278:1 +T 5.9% intergenic (+266/‑20) tst_1 → / → speC_1 Toxic shock syndrome toxin‑1/Exotoxin type C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE447,2781.T 5.9% 837.9 / 43.6 203intergenic (+266/‑20)tst_1/speC_1Toxic shock syndrome toxin‑1/Exotoxin type C
Reads supporting (aligned to +/- strand):  ref base . (99/92);  new base T (4/8);  total (103/100)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 222 total aligned reads displayed.

GTTAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTACAAC  >  USA300TCH1516_ALE/447189‑447368
                                                                                          |                                                                                          
gTTAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTataat                                                                                           >  1:324561/1‑91 (MQ=255)
gTTAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTataat                                                                                           >  1:3414069/1‑91 (MQ=255)
gTTAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTataat                                                                                           <  1:1724697/91‑1 (MQ=255)
 ttAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTataata                                                                                          <  2:3673460/91‑1 (MQ=255)
  tAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAg                                                                                         <  1:2241752/91‑1 (MQ=255)
   aaaaaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGa                                                                                        >  2:1599061/1‑91 (MQ=255)
   aaaaaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGa                                                                                        >  1:1978500/1‑91 (MQ=255)
     aaaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTATTATAATTAGAtt                                                                                      <  2:2901459/91‑1 (MQ=39)
      aaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTg                                                                                     <  1:3621667/91‑1 (MQ=255)
      aaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTg                                                                                     >  1:2323377/1‑91 (MQ=255)
      aaCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTATTATAATTAGATTg                                                                                     <  1:2840325/91‑1 (MQ=38)
       ccAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTATTATAATTAGATTgt                                                                                    <  2:3850839/90‑2 (MQ=38)
       aCAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTATTATAATTAGATTgt                                                                                    <  2:1243884/91‑2 (MQ=38)
        cAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGa                                                                                   <  1:8602/91‑1 (MQ=255)
          ggTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGc                                                                                 >  2:459587/1‑91 (MQ=255)
          ggTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGc                                                                                 >  2:3663596/1‑91 (MQ=255)
          ggTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGc                                                                                 >  1:2143609/1‑91 (MQ=255)
           gTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCa                                                                                <  2:3558658/91‑1 (MQ=255)
            ttAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCAt                                                                               >  2:1908652/1‑91 (MQ=255)
            ttAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCAt                                                                               <  2:3393993/91‑1 (MQ=255)
            ttAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCAt                                                                               <  1:2940799/91‑1 (MQ=255)
            ttAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCAt                                                                               <  1:808299/91‑1 (MQ=255)
             tAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATa                                                                              <  2:1079065/91‑1 (MQ=255)
                    aaTATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATGGGAGCATACAACTat                                                                       <  1:3309768/91‑1 (MQ=255)
                     atatCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatg                                                                      <  1:811410/91‑1 (MQ=255)
                      tatCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatga                                                                     >  2:1914058/1‑91 (MQ=255)
                         ccGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaa                                                                  <  1:267641/91‑1 (MQ=255)
                          cGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaat                                                                 >  2:3773954/1‑91 (MQ=255)
                          cGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaat                                                                 >  2:2522102/1‑91 (MQ=255)
                           gAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatg                                                                <  2:243499/91‑1 (MQ=255)
                           gAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatg                                                                >  2:2323351/1‑91 (MQ=255)
                           gAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatg                                                                <  1:3208077/91‑1 (MQ=255)
                            aaTACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatga                                                               <  1:3556183/91‑1 (MQ=255)
                            aaTACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATCGGAGCATACAACTatgaaaatga                                                               >  2:913059/1‑91 (MQ=255)
                             aTACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACACTCTAATTATAATAGATTGGAGCATACAACTatgaaaatgaa                                                              >  1:3633772/1‑91 (MQ=255)
                             aTACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatgaa                                                              <  2:3864092/91‑1 (MQ=255)
                              tACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTatgaaaatgaaa                                                             <  2:3211472/91‑1 (MQ=255)
                               aCAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTAtgaaaatgaaaa                                                            >  2:2621091/1‑91 (MQ=255)
                               aCAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTAtgaaaatgaaaa                                                            <  2:2198089/91‑1 (MQ=255)
                                  gCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAata                                                         >  1:1834267/1‑91 (MQ=255)
                                    tCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATAtt                                                       >  2:500609/1‑91 (MQ=255)
                                     ccTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTg                                                      <  1:3818329/91‑1 (MQ=255)
                                     ccTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTg                                                      >  2:814503/1‑91 (MQ=255)
                                      cTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGc                                                     <  1:3756909/91‑1 (MQ=255)
                                       tataATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCa                                                    <  2:114533/91‑1 (MQ=255)
                                          aataTGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCaaaa                                                 >  2:2781070/1‑91 (MQ=255)
                                           atatGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCaaaaa                                                >  1:1138580/1‑91 (MQ=255)
                                           atatGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCaaaaa                                                <  1:2516133/91‑1 (MQ=255)
                                            tatGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAAt                                               >  1:468577/1‑91 (MQ=255)
                                               ggggGTATGATTGAAAGTACGGAATAAAACAATCAAATTATAATAGATTGGAGCGTACAACTATGAAAATGAAAAATATTGCAAAAATAAg                                            <  2:168812/91‑1 (MQ=255)
                                               gggTGTATGATTCAAATTACGTAATAAAACACTCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAg                                            >  1:3447202/1‑91 (MQ=255)
                                               gggTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAg                                            >  2:594873/1‑91 (MQ=255)
                                                ggTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGt                                           <  2:741665/91‑1 (MQ=255)
                                                ggTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGt                                           <  2:3076960/91‑1 (MQ=255)
                                                 gtgtATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGtt                                          <  1:761784/91‑1 (MQ=255)
                                                 gtgtATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGtt                                          <  2:1753406/91‑1 (MQ=255)
                                                   gtATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTg                                        >  2:2081546/1‑91 (MQ=255)
                                                    tATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGt                                       <  1:3391256/91‑1 (MQ=255)
                                                    tATGATTCAAATTACGTAATAAAACAATCTAATTATAATACATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTCTGt                                       >  2:3189749/1‑91 (MQ=255)
                                                      tGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGtta                                     >  1:1094579/1‑91 (MQ=255)
                                                       gATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGttat                                    <  1:488901/91‑1 (MQ=255)
                                                       gATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGttat                                    <  1:548464/91‑1 (MQ=255)
                                                         ttCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGttatta                                  <  2:3546854/91‑1 (MQ=255)
                                                          tCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAg                                 >  1:107314/1‑91 (MQ=255)
                                                           cAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAgg                                >  2:1010559/1‑91 (MQ=255)
                                                           cAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAgg                                >  1:1302450/1‑91 (MQ=255)
                                                              aTTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAat                             >  1:195207/1‑91 (MQ=255)
                                                              aTTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAat                             <  2:2495656/91‑1 (MQ=255)
                                                                tACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAatat                           >  1:938566/1‑91 (MQ=255)
                                                                tACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAatat                           >  1:1549972/1‑91 (MQ=255)
                                                                 aCGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAatat                           >  2:2693831/1‑90 (MQ=255)
                                                                  cGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAAACAAGAAAAAATAAGTTTGTTATTAGGAATATTa                         <  2:3804592/91‑1 (MQ=38)
                                                                    taataaAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGc                       <  2:1466707/91‑1 (MQ=255)
                                                                    taataaAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGc                       >  2:2606266/1‑91 (MQ=255)
                                                                      ataaAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCaa                     >  1:2751946/1‑91 (MQ=255)
                                                                       taaAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAAc                    <  1:2079944/91‑1 (MQ=255)
                                                                         aaaCAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAg                  <  1:45377/91‑1 (MQ=255)
                                                                            cAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGgtg               <  2:2090997/91‑1 (MQ=255)
                                                                            cAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGgtg               <  2:1737756/91‑1 (MQ=255)
                                                                             aaTCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGgtgt              <  1:2522102/91‑1 (MQ=255)
                                                                              aTCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTa             >  2:2054009/1‑91 (MQ=255)
                                                                               tCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTaa            <  1:2674366/91‑1 (MQ=255)
                                                                               tCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTaa            >  1:1566076/1‑91 (MQ=255)
                                                                                cTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTaaa           <  1:1319749/91‑1 (MQ=255)
                                                                                cTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTaaa           >  1:2004467/1‑91 (MQ=255)
                                                                                  aaTTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAaca         <  2:1865498/91‑1 (MQ=255)
                                                                                   aTTATAATTAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAaca         >  1:151245/1‑91 (MQ=255)
                                                                                   aTTATAATTAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAaca         >  2:1336215/1‑91 (MQ=255)
                                                                                   aTTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAacac        >  2:3298930/1‑91 (MQ=255)
                                                                                   aTTATAATAGATTGGAGCATACAACTATGAAAATAAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAacac        >  2:1813023/1‑91 (MQ=255)
                                                                                    ttATAATTAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAacac        >  1:2294668/1‑91 (MQ=255)
                                                                                    ttATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACt       <  2:1164101/91‑1 (MQ=255)
                                                                                    ttATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACt       >  2:2230513/1‑91 (MQ=255)
                                                                                     tataATAGATTGGAGCATACAACTATGAAAATGAAAAATCTTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTa      >  1:439812/1‑91 (MQ=255)
                                                                                      ataataGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTAc     <  1:1999437/91‑1 (MQ=255)
                                                                                      ataataGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTAc     >  2:2806985/1‑91 (MQ=255)
                                                                                      ataataGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTAc     >  2:556693/1‑91 (MQ=255)
                                                                                       taataGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTACa    <  1:768437/91‑1 (MQ=255)
                                                                                        aataGATTGGAGCATAGAACTATGAAGATGAAAAAAAATGAAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTACaa   <  2:3493723/91‑1 (MQ=255)
                                                                                        aataGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTACaa   <  2:1730231/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTTAAAAACAGGTTAATGTGAATATCCGAATACAGCTCCTATAATATGGGTGTATGATTCAAATTACGTAATAAAACAATCTAATTATAATAGATTGGAGCATACAACTATGAAAATGAAAAATATTGCAAAAATAAGTTTGTTATTAGGAATATTAGCAACAGGTGTAAACACTACAAC  >  USA300TCH1516_ALE/447189‑447368

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