Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,158 T→A 13.8% intergenic (+65/‑66) 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,1580TA13.8% 589.6 / 65.7 215intergenic (+65/‑66)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (86/99);  new base A (13/17);  total (99/116)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248
                                                                                          |                                                                                          
tAGCATACAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                            <  2:2520380/91‑1 (MQ=255)
 aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                           <  2:2755144/91‑1 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGacac                                                                                         <  1:283555/91‑1 (MQ=255)
    aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACt                                                                                        <  2:1978974/91‑1 (MQ=255)
     ttCAAAAAAGGTAAATACTAGGCGTTAAAGAAAGTTTCAGCTAAACACATAAAATTAAGTTATTCATTTATAAAGCCCTTCAGTGACACTg                                                                                       <  2:3814955/91‑1 (MQ=255)
     tgcAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAAGTAAATGATTCATTTATAAAGTCCTTCAGTTTCACTg                                                                                       <  2:3645226/89‑1 (MQ=255)
     gtCAAAAAAGGTAAAATCTAAGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTg                                                                                       <  2:1611255/90‑1 (MQ=255)
      tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                                      >  2:3042337/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  2:3865276/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  1:997975/1‑91 (MQ=255)
        aaaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAg                                                                                    <  2:2783536/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:3533696/91‑1 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 <  1:1706233/91‑1 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 <  1:573091/91‑1 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 >  1:1334519/1‑91 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 <  2:2974808/91‑1 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 >  1:792778/1‑91 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAggac                                                                                 <  1:2027717/91‑3 (MQ=38)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCGAAACACATAAAATTAAATTATTCATTTATAAAGCCCGTCAGTGACACTGAAggac                                                                                 <  2:695599/91‑3 (MQ=38)
           aaaGGGAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAggac                                                                                 <  2:3139430/91‑3 (MQ=38)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              >  1:700305/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAAAACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:1003859/1‑91 (MQ=255)
                  aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttg                                                                          <  2:2434904/91‑1 (MQ=255)
                   aTTCTAGGTTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                         <  2:3371292/91‑1 (MQ=255)
                   aTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                         >  1:131515/1‑91 (MQ=255)
                    ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                        <  2:350646/91‑1 (MQ=255)
                     tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGAGCTTTAttgttg                                                                       <  1:2994656/91‑1 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    >  1:3552657/1‑91 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    >  2:3737420/1‑91 (MQ=255)
                         ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                   <  1:2603737/91‑1 (MQ=255)
                         ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                   <  1:2885542/91‑1 (MQ=255)
                          gCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACaa                                                                  >  1:1487424/1‑91 (MQ=255)
                           cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATGCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                 <  2:1752490/91‑1 (MQ=255)
                            tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                >  1:1125130/1‑91 (MQ=255)
                            tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                                <  2:191063/91‑1 (MQ=255)
                              taaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTg                                                              <  2:1110702/91‑1 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             <  2:855223/91‑1 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTAGAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             <  1:3079121/91‑1 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            >  1:63981/1‑91 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            <  2:2173775/91‑1 (MQ=255)
                                 ataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCTTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAg                                                           <  1:3356972/91‑1 (MQ=255)
                                      gTTTCAGCTAAACACATAAAATTAGATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTGATTGTTGTACGATTTGAAGTTAGt                                                      <  2:1378751/91‑1 (MQ=255)
                                          cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                  >  1:454992/1‑91 (MQ=255)
                                          cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                  >  1:1752442/1‑91 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  1:1702718/91‑1 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               >  2:1027917/1‑91 (MQ=255)
                                              tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                              >  2:896022/1‑91 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            <  2:871244/91‑1 (MQ=255)
                                                  cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                          <  2:793477/91‑1 (MQ=255)
                                                  cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                          <  2:587641/91‑1 (MQ=255)
                                                   acaTAAAATTAAATTATTCATTTATACAGTCCTTCAGTGTGACCGAAGGGCCTTATTGTTGTACAATTTGCAGTGAGTAACCTATCAGATc                                         >  2:805770/1‑91 (MQ=255)
                                                   acaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATc                                         >  2:805929/1‑91 (MQ=255)
                                                     aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                                       <  1:194601/91‑1 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    <  1:2185438/91‑1 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    <  1:2185773/91‑1 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    >  1:2921410/1‑91 (MQ=255)
                                                         aattaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagga                                   <  2:2850386/91‑1 (MQ=255)
                                                         aattaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagga                                   >  2:1701101/1‑91 (MQ=255)
                                                            taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGGAGTGAGTAACATATCAAATCTAAGGAGGt                                >  2:1964266/1‑91 (MQ=255)
                                                            taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                                <  2:3300120/91‑1 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              >  1:3427786/1‑91 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              >  1:2050150/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            >  1:2695495/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            <  2:3594861/91‑1 (MQ=255)
                                                                 tattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAtt                           <  1:2311418/91‑1 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          >  1:3301303/1‑91 (MQ=255)
                                                                   ttCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTAt                         <  2:2618045/91‑1 (MQ=255)
                                                                    tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                        <  2:1182518/91‑1 (MQ=255)
                                                                    tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                        >  1:3188573/1‑91 (MQ=255)
                                                                     cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa                       <  1:3165318/91‑1 (MQ=255)
                                                                       tttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAt                     <  2:1022224/91‑1 (MQ=255)
                                                                       tttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAt                     >  2:2682141/1‑91 (MQ=255)
                                                                        ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtt                    <  2:3578668/91‑1 (MQ=255)
                                                                         tataAGGTCCTTCAGTGTCACTGAAGGGCTGTATTGTTGGACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   <  2:461710/91‑1 (MQ=255)
                                                                         tataAAGTCCTTCAGTGTCACTGAAGGTCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATTAAttt                   <  2:115766/91‑1 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCTAATCTAAGGAGGTACATTATGAATTttaa                >  1:3116829/1‑91 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                <  1:738948/91‑1 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                <  2:370835/91‑1 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                <  2:3066785/91‑1 (MQ=255)
                                                                             aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa               >  2:2312143/1‑91 (MQ=255)
                                                                              agcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTATGGAGGTACATTATGAATTttaaat              >  2:3606360/4‑91 (MQ=38)
                                                                              agcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              >  2:1428744/4‑91 (MQ=38)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              <  1:382460/91‑1 (MQ=255)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTACGAATTttaaat              <  2:1593524/91‑1 (MQ=255)
                                                                               gcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt             >  2:2936031/3‑91 (MQ=255)
                                                                               gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt             >  2:3413738/1‑91 (MQ=255)
                                                                               gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt             >  2:2570754/1‑91 (MQ=255)
                                                                                 ccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa           >  1:1018752/1‑91 (MQ=255)
                                                                                 ccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa           >  1:3705130/1‑91 (MQ=255)
                                                                                 ccTTCAGTGACACTGAAGGACTTTATCGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa           >  2:2211247/1‑91 (MQ=39)
                                                                                   ttCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  2:3774361/1‑91 (MQ=255)
                                                                                   ttCAGTGACACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  1:1277993/1‑91 (MQ=255)
                                                                                    tCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa        >  2:2700187/1‑91 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       <  2:2382244/91‑1 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       <  2:1329466/91‑1 (MQ=255)
                                                                                      aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      >  1:2404019/1‑91 (MQ=255)
                                                                                      aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      <  1:2896259/91‑1 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     >  1:1793829/1‑91 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     >  2:710727/1‑91 (MQ=255)
                                                                                         gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca   <  1:2510179/91‑1 (MQ=255)
                                                                                          |                                                                                          
TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248

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: