Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,168 G→A 10.4% intergenic (+75/‑56) USA300HOU_RS04665 → / → pepA_1 NADH dehydrogenase‑like protein/Cytosol aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE937,1680GA10.4% 720.0 / 47.1 245intergenic (+75/‑56)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (115/103);  new base A (23/3);  total (138/107)
Fisher's exact test for biased strand distribution p-value = 5.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

AAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGA  >  USA300TCH1516_ALE/937078‑937258
                                                                                          |                                                                                          
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           <  2:3957015/91‑1 (MQ=255)
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           <  2:3033572/91‑1 (MQ=255)
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           >  1:2186674/1‑91 (MQ=255)
    ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                                        >  2:1064816/1‑91 (MQ=255)
    ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                                        >  1:2119079/1‑91 (MQ=255)
    ggTAAATTCTAGGCATTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTCCAGTGACACTGAAGGACttt                                                                                        >  2:1175530/1‑91 (MQ=39)
     gTAAAGTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAGAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTa                                                                                       <  2:2643814/91‑1 (MQ=255)
          ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                                  <  2:1869609/91‑1 (MQ=255)
          ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                                  <  1:1677493/91‑1 (MQ=255)
           tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                                 <  2:1033758/91‑1 (MQ=255)
             tAGGCTTTAAATAAAGTTTCAGCTAAGCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCGCTGAAGGGCTTTATTGTTGTa                                                                               <  2:1044917/91‑1 (MQ=255)
              aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATATAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                              <  2:3632641/91‑1 (MQ=255)
              aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                              <  2:1638898/91‑1 (MQ=255)
               ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                             >  1:413625/1‑91 (MQ=255)
               ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                             <  1:1537495/91‑1 (MQ=255)
               ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                             <  1:1730218/91‑1 (MQ=255)
                 cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                           >  1:1737014/1‑91 (MQ=255)
                 cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                           <  2:1490072/91‑1 (MQ=255)
                 cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                           >  1:1154746/1‑91 (MQ=255)
                  tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                          <  1:1196371/91‑1 (MQ=255)
                  tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                          >  1:2901081/1‑91 (MQ=255)
                     aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCCTTTCTAAAGTCCTTCAGTGTCACTGACGGGCTCTATTGTTGTACAATTTGa                                                                       >  1:1229529/1‑91 (MQ=255)
                      aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                                      >  1:2097395/1‑91 (MQ=255)
                      aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                                      <  2:2794875/91‑1 (MQ=255)
                        taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGt                                                                    <  1:2296246/91‑1 (MQ=255)
                         aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                                   <  2:269203/91‑1 (MQ=255)
                           aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                                 <  1:1091313/91‑1 (MQ=255)
                           aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                                 >  2:3921189/1‑91 (MQ=255)
                           aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                                 >  1:2558156/1‑91 (MQ=255)
                           aGTTTCAGCTAAACACATAAAATTAAATTATACATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTAATTGTTGTACAGTTTGAAGTGAg                                                                 >  2:3024596/1‑91 (MQ=255)
                              ttCAGCTAAACCCATAAAATTAAATGATGCATTTATAGAGTCCGTCAGTGTCACTGAAGGGCGTTATTGTTGTACAGTTTGAAGTGAGTaa                                                              <  2:179115/91‑1 (MQ=255)
                              ttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTTTTGTTGTACAATTTGAAGTGAGTaa                                                              >  2:1722832/1‑91 (MQ=255)
                               tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGATGTGAGTAAc                                                             >  1:3545041/1‑91 (MQ=255)
                                cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                            >  1:2685649/1‑91 (MQ=255)
                                 aGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACat                                                           <  1:1361866/91‑1 (MQ=255)
                                  gCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACata                                                          >  1:2155079/1‑91 (MQ=255)
                                  gCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACata                                                          <  1:2473776/91‑1 (MQ=255)
                                    tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                                        <  2:2784800/91‑1 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTCACATATCaa                                                      >  1:3788938/1‑91 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                                      >  1:3451289/1‑91 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                                      <  2:151126/91‑1 (MQ=255)
                                       acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                                     >  2:503436/1‑91 (MQ=255)
                                       acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                                     <  1:2263980/91‑1 (MQ=255)
                                       acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                                     <  2:297277/91‑1 (MQ=255)
                                        cacgTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGCAACATATCAAAt                                                    <  1:3056110/91‑1 (MQ=255)
                                        cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                                    >  1:838738/1‑91 (MQ=255)
                                          caTAAAATTGAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                                  <  1:3790028/91‑1 (MQ=255)
                                          caTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                                  <  1:1139211/91‑1 (MQ=255)
                                            tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                                >  2:245619/1‑91 (MQ=255)
                                             aaaaTTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAag                                               <  2:1307330/91‑1 (MQ=255)
                                              aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                              >  1:3864574/1‑91 (MQ=255)
                                                attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                            <  1:3926914/91‑1 (MQ=255)
                                                attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATAGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                            <  1:3419759/91‑1 (MQ=255)
                                                  taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                                          >  2:3764023/1‑91 (MQ=255)
                                                   aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGCGTAACATATCAAATCTAAGGAGGTa                                         >  2:1841221/1‑91 (MQ=255)
                                                   aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                                         <  1:1856815/91‑1 (MQ=255)
                                                   aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                                         >  1:1326830/1‑91 (MQ=255)
                                                   aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                                         <  1:3041173/91‑1 (MQ=255)
                                                    aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                                        >  1:505831/1‑91 (MQ=255)
                                                     attattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACa                                       >  1:841006/1‑91 (MQ=255)
                                                     attattCATTTATAAAGTCCTAAAGTGCCACTGAAGGGTTTTATTGTTGTACAATTTGAAGTGAGTAACATATGAAATCTAGGGAGGAAcg                                       >  2:3524273/1‑90 (MQ=255)
                                                      ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                                      >  2:566438/1‑91 (MQ=255)
                                                        attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                                    >  1:2823086/1‑91 (MQ=255)
                                                         ttCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGATACATTAt                                   <  2:330033/91‑1 (MQ=255)
                                                           cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa                                 <  2:238878/91‑1 (MQ=255)
                                                           cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa                                 <  1:1258590/91‑1 (MQ=255)
                                                            aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                                <  2:1195544/91‑1 (MQ=255)
                                                              ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGCAGTGAGTAACATATCAACTCTAAGGAGGTACATCATGAAtt                              >  2:3487160/1‑91 (MQ=255)
                                                              ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtt                              <  1:3317091/91‑1 (MQ=255)
                                                               tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                             <  1:1571243/91‑1 (MQ=255)
                                                                 taAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta                           >  1:142448/1‑91 (MQ=38)
                                                                   aaGTCCTTTATTTTCACTTAAGGTCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCACATCTAAGGAGGTACATTATGAATTttaaa                         <  2:3435252/91‑1 (MQ=255)
                                                                    agcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat                        >  1:1891829/4‑91 (MQ=38)
                                                                    aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat                        <  1:2398576/91‑1 (MQ=255)
                                                                     gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt                       <  1:403989/91‑1 (MQ=255)
                                                                      tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta                      >  2:1955424/1‑91 (MQ=255)
                                                                      tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACACATCAAATCTCAGGAGGTACATTATGAATTttaaatta                      >  2:1390596/1‑91 (MQ=255)
                                                                      cccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta                      >  2:2485417/2‑91 (MQ=255)
                                                                       ccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa                     >  2:1926323/1‑91 (MQ=255)
                                                                          tCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa                  >  2:2184249/1‑91 (MQ=255)
                                                                           cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                 >  2:739819/1‑91 (MQ=255)
                                                                           cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                 <  2:23255/91‑1 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  2:2206666/1‑91 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  1:2448692/1‑91 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  1:264808/1‑91 (MQ=255)
                                                                            aGTGACACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  1:621578/1‑91 (MQ=255)
                                                                            aGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  1:1133678/1‑91 (MQ=255)
                                                                              tgacACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacac              >  2:1260766/4‑91 (MQ=255)
                                                                               gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATGAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca             <  2:3056926/91‑1 (MQ=255)
                                                                                acACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGCGGTACATTATGAATTTTAAATTAAATAacacac            >  1:3234605/2‑91 (MQ=255)
                                                                                 cACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCANATCTCAGGAGGTACATTATNAATTTTAAATTAAATAACACACt           >  2:1485034/1‑91 (MQ=255)
                                                                                 cACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACt           >  2:461115/1‑91 (MQ=255)
                                                                                  aCTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTa          <  1:3351962/91‑1 (MQ=255)
                                                                                  aCTGAAGGACTTTACTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTa          >  2:3354234/1‑91 (MQ=255)
                                                                                       aGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCaa     <  2:2951292/91‑1 (MQ=255)
                                                                                        gggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAAc    <  1:3576760/91‑1 (MQ=255)
                                                                                        gggCTTTATTGGTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAAc    <  1:3590868/91‑1 (MQ=255)
                                                                                         ggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACg   >  1:3727269/1‑91 (MQ=255)
                                                                                         ggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACg   <  1:1920898/91‑1 (MQ=255)
                                                                                          gCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGCGGTACATTATGAATTTTACATTAAATAACGCACTAAGCAACGa  >  1:3952365/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGA  >  USA300TCH1516_ALE/937078‑937258

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