Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,168 G→A 6.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,1680GA 6.4% 630.8 / 18.1 200intergenic (+75/‑56)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (89/97);  new base A (11/2);  total (100/100)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

AAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGA  >  USA300TCH1516_ALE/937078‑937258
                                                                                          |                                                                                          
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           <  1:573091/91‑1 (MQ=255)
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           <  2:2974808/91‑1 (MQ=255)
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                           >  1:1334519/1‑91 (MQ=255)
 aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAACTCCTTCAGTGTCACTGAAGGGc                                                                                           >  2:3119746/1‑91 (MQ=255)
   aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                                         >  2:47710/1‑91 (MQ=255)
   aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                                         <  2:3746484/91‑1 (MQ=255)
   aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACtt                                                                                         <  2:2637443/91‑1 (MQ=38)
       aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAAAACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                                     >  1:1003859/1‑91 (MQ=255)
         aTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                                   >  1:131515/1‑91 (MQ=255)
           tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                                 <  1:2604360/91‑1 (MQ=255)
           tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGAGCTTTAttgttg                                                                                 <  1:2994656/91‑1 (MQ=255)
            cTAGGCTTTAAATAAAGTTTCAGCTAAACACATTAAATTAAATTATTCATTTATAAAGTGCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                                                <  1:3450204/91‑1 (MQ=255)
              aGGCTTTCAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                              <  2:3452143/91‑1 (MQ=255)
              aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                              >  2:3737420/1‑91 (MQ=255)
                gCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACaa                                                                            >  1:1487424/1‑91 (MQ=255)
                 cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                           <  2:360844/91‑1 (MQ=255)
                 cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                           <  1:2134025/91‑1 (MQ=255)
                  tttAAATAAAGTTTCAGCTAAACGCATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                          <  2:3550725/91‑1 (MQ=255)
                  tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                          >  1:3747468/1‑91 (MQ=255)
                    taaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTg                                                                        <  2:1110702/91‑1 (MQ=255)
                     aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTAGAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                                       <  1:3079121/91‑1 (MQ=255)
                       ataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCTTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAg                                                                     <  1:3356972/91‑1 (MQ=255)
                        taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGt                                                                    >  2:738948/1‑91 (MQ=255)
                        taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGt                                                                    <  1:1690319/91‑1 (MQ=255)
                         aaaGTTTCTGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                                   >  2:3327986/1‑91 (MQ=255)
                         aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                                   >  1:1113841/1‑91 (MQ=255)
                            gTTTCAGCTAAACACATAAAATTAGATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTGATTGTTGTACGATTTGAAGTTAGt                                                                <  2:1378751/91‑1 (MQ=255)
                               tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                             >  2:2188104/1‑91 (MQ=255)
                                   cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                                         <  1:1702719/91‑1 (MQ=255)
                                   cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                                         >  2:1027917/1‑91 (MQ=255)
                                   cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                                         >  2:39011/1‑91 (MQ=255)
                                     aaaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCa                                                       >  1:783908/1‑91 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                                      <  2:871174/91‑1 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                                      <  2:871397/91‑1 (MQ=255)
                                      aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGACACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                                      >  2:1164803/1‑91 (MQ=255)
                                       acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                                     <  1:2866460/91‑1 (MQ=255)
                                        cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                                    <  2:793477/91‑1 (MQ=255)
                                         acaTAAAATTAAATTATTCATTTATACAGTCCTTCAGTGTGACCGAAGGGCCTTATTGTTGTACAATTTGCAGTGAGTAACCTATCAGATc                                                   >  2:805770/1‑91 (MQ=255)
                                         acaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATc                                                   <  1:1814126/91‑1 (MQ=255)
                                         acaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATc                                                   >  2:805929/1‑91 (MQ=255)
                                         acaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATc                                                   <  1:1142362/91‑1 (MQ=255)
                                          caTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                                  >  2:2415371/1‑91 (MQ=255)
                                           aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                                                 <  1:194601/91‑1 (MQ=255)
                                              aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                              <  1:3684763/91‑1 (MQ=255)
                                              aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                              <  1:2185773/91‑1 (MQ=255)
                                              aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                              >  1:2921959/1‑91 (MQ=255)
                                               aattaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagga                                             <  2:2850386/91‑1 (MQ=255)
                                                attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTCaggag                                            >  2:3587966/1‑91 (MQ=255)
                                                attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                            <  2:1719756/91‑1 (MQ=255)
                                                 ttaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggagg                                           >  2:919067/1‑91 (MQ=255)
                                                  taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGGAGTGAGTAACATATCAAATCTAAGGAGGt                                          >  2:1964266/1‑91 (MQ=255)
                                                    aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                                        >  1:2050150/1‑91 (MQ=255)
                                                    aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                                        >  1:3427665/1‑91 (MQ=255)
                                                     attattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACa                                       <  2:2720044/91‑1 (MQ=255)
                                                      ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                                      >  1:2695495/1‑91 (MQ=255)
                                                       tattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGGTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAtt                                     <  2:3056167/91‑1 (MQ=255)
                                                       tattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAtt                                     <  1:2311418/91‑1 (MQ=255)
                                                          tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                                  <  1:2586885/91‑1 (MQ=255)
                                                          tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                                  >  1:3188573/1‑91 (MQ=255)
                                                          tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                                  <  2:2613838/91‑1 (MQ=255)
                                                           cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa                                 <  1:3165318/91‑1 (MQ=255)
                                                             tttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAt                               <  2:1022258/91‑1 (MQ=255)
                                                              ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtt                              <  2:2474239/91‑1 (MQ=255)
                                                                 taAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta                           >  1:2275537/1‑91 (MQ=255)
                                                                 taAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta                           >  1:243360/1‑91 (MQ=255)
                                                                  aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                          <  2:370835/91‑1 (MQ=255)
                                                                  aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                          <  1:738948/91‑1 (MQ=255)
                                                                  aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                          >  2:1905655/1‑91 (MQ=255)
                                                                   aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa                         >  2:2312143/1‑91 (MQ=255)
                                                                    aGTGCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat                        <  1:3134842/91‑1 (MQ=255)
                                                                    aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat                        <  1:641392/91‑1 (MQ=255)
                                                                    aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat                        <  1:382460/91‑1 (MQ=255)
                                                                    aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTACGAATTttaaat                        <  2:1593524/91‑1 (MQ=255)
                                                                     gcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt                       >  2:2936031/3‑91 (MQ=255)
                                                                     gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt                       >  2:2570754/1‑91 (MQ=255)
                                                                     gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTCGAAGTGCGTAACATATCAAATCTAAAAAGGTACATTATGAATTttaaatt                       >  2:3603897/1‑91 (MQ=255)
                                                                      cccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATAACAAATCTAAGGAGGTACATTATGAATTttaaatta                      >  1:1719609/2‑91 (MQ=38)
                                                                         ttCAGTGACACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat                   >  2:2111789/1‑91 (MQ=255)
                                                                          tCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa                  >  2:2700101/1‑91 (MQ=255)
                                                                           cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                 >  1:1762237/1‑91 (MQ=255)
                                                                           cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                 <  2:2382244/91‑1 (MQ=255)
                                                                           cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                 <  2:1329466/91‑1 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                <  1:2896259/91‑1 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  1:2404019/1‑91 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                <  2:702952/91‑1 (MQ=255)
                                                                            aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac                >  2:841869/1‑91 (MQ=255)
                                                                              tgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacac              <  1:13083/91‑1 (MQ=255)
                                                                               gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAGCATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca             <  2:3412314/91‑1 (MQ=255)
                                                                               gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca             >  1:519986/1‑91 (MQ=255)
                                                                               gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca             >  1:808501/1‑91 (MQ=255)
                                                                                 cACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACt           <  1:712737/91‑1 (MQ=255)
                                                                                    tGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAg        <  2:1762558/91‑1 (MQ=255)
                                                                                    tGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAg        >  2:2812712/1‑91 (MQ=255)
                                                                                     gAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGc       >  1:1010970/1‑91 (MQ=255)
                                                                                        gggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAAc    >  1:3270576/1‑91 (MQ=255)
                                                                                         ggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACg   <  2:3734310/91‑1 (MQ=255)
                                                                                         ggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACg   >  1:107180/1‑91 (MQ=255)
                                                                                         ggCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACg   >  2:1211025/1‑91 (MQ=255)
                                                                                          gCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTaa         >  2:2291521/1‑84 (MQ=255)
                                                                                          gCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGa  <  2:1515255/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGA  >  USA300TCH1516_ALE/937078‑937258

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: