Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,158 T→A 16.1% 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,1580TA16.1% 677.0 / 109.0 261intergenic (+65/‑66)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (109/109);  new base A (22/20);  total (132/129)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 276 total aligned reads displayed.

TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248
                                                                                          |                                                                                          
tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                            >  1:1263901/1‑91 (MQ=255)
tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTAATTTATAAAGTCCTTCAgtgt                                                                                            >  2:473133/1‑91 (MQ=255)
 aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATGAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                           >  2:3902961/1‑91 (MQ=255)
 aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                           <  2:3459499/91‑1 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGacac                                                                                         <  2:2429805/91‑1 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGacac                                                                                         <  1:2416374/91‑1 (MQ=255)
    aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                                        <  2:2658237/91‑1 (MQ=255)
    aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                                        <  2:1467843/91‑1 (MQ=255)
      tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                                      >  1:271272/1‑91 (MQ=255)
       cagaaaaGGTAAATTCTATGCTTTAAATAAAGTTTCAGCTAAACACATAGAATTAAATTATGCATTTATAAAGCCCTTCAGTGACACTGaa                                                                                     <  2:1814933/88‑1 (MQ=38)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  2:2464383/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGaa                                                                                     <  1:1670324/91‑1 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGaa                                                                                     <  1:3652043/91‑1 (MQ=255)
        aaaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                                    >  2:1918051/1‑91 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:1508803/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:2242570/91‑1 (MQ=255)
          aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAggg                                                                                  <  1:2579496/91‑1 (MQ=255)
          aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                                  <  1:3291582/91‑2 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 <  2:3033680/91‑1 (MQ=255)
             aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACtt                                                                               <  1:155035/91‑1 (MQ=38)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              >  1:2119079/1‑91 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              >  2:1743678/1‑91 (MQ=255)
               gTAAAGTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAGAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTa                                                                             <  2:2643814/91‑1 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           <  1:3592571/91‑1 (MQ=255)
                   aTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                         <  1:380271/91‑1 (MQ=255)
                     tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                       <  2:313532/91‑1 (MQ=255)
                     tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                       <  2:1034111/91‑1 (MQ=255)
                       tAGGCTTTAAATAAAGTTTCAGCTAAGCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCGCTGAAGGGCTTTATTGTTGTa                                                                     <  2:1044917/91‑1 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    >  2:1501427/1‑91 (MQ=255)
                          gCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACaa                                                                  >  1:1723948/1‑91 (MQ=255)
                           cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                 >  1:1736900/1‑91 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCCTTTCTAAAGTCCTTCAGTGTCACTGACGGGCTCTATTGTTGTACAATTTGa                                                             >  1:1229529/1‑91 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             <  2:335382/91‑1 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             >  1:1871164/1‑91 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             <  1:3545268/91‑1 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            >  1:2097395/1‑91 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            >  1:2555464/1‑91 (MQ=255)
                                  taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGt                                                          <  2:3571774/91‑1 (MQ=255)
                                  taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGt                                                          <  1:2296246/91‑1 (MQ=255)
                                   aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                         >  2:2365833/1‑91 (MQ=255)
                                     aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                       <  1:1091313/91‑1 (MQ=255)
                                     aGTTTCAGCTAAACACATAAAATTAAATTATACATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTAATTGTTGTACAGTTTGAAGTGAg                                                       >  2:3024799/1‑91 (MQ=255)
                                     aGTTTAAGATAAATACATAAAATTAAATTATTAATTTATAAAGTCCTTCAGTGTTAATGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                       >  2:364832/1‑91 (MQ=255)
                                       tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                                     >  2:1867945/1‑91 (MQ=255)
                                        ttCAGCTAAACCCATAAAATTAAATGATGCATTTATAGAGTCCGTCAGTGTCACTGAAGGGCGTTATTGTTGTACAGTTTGAAGTGAGTaa                                                    <  2:179115/91‑1 (MQ=255)
                                         tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                   <  2:2161350/91‑1 (MQ=255)
                                          cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                  >  1:3473312/1‑91 (MQ=255)
                                            gCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACata                                                <  1:2473669/91‑1 (MQ=255)
                                              tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                              <  2:2784800/91‑1 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTGACATATCaa                                            >  1:3788939/1‑91 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            >  1:1199489/1‑91 (MQ=255)
                                                 acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                           <  2:297277/91‑1 (MQ=255)
                                                 acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                           <  1:2263980/91‑1 (MQ=255)
                                                  cacaTAAAATTGAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGCTGTACAATTTGAAGTGAGTAACATATCAAAt                                          <  2:3024343/91‑1 (MQ=255)
                                                  cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                          >  1:838738/1‑91 (MQ=255)
                                                  cacaTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGGAACATATCAAAt                                          <  1:3056224/91‑1 (MQ=255)
                                                    caTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                        <  1:1139211/91‑1 (MQ=255)
                                                     aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                                       >  1:412395/1‑91 (MQ=255)
                                                      tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                      >  2:245861/1‑91 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    >  1:3864574/1‑91 (MQ=255)
                                                         aattaaattaTTCATTTATAAAGTCCTTCAGAGTCACTGAAGGGCTTTATTGTTGTACACTTGGAAGTGAGTAACATATCAAATCTAagga                                   >  2:675429/1‑91 (MQ=255)
                                                          attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                  <  1:3927034/91‑1 (MQ=255)
                                                          attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATAGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                  <  1:3419759/91‑1 (MQ=255)
                                                            taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                                <  1:2662812/91‑1 (MQ=255)
                                                            taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                                >  1:1802090/1‑91 (MQ=255)
                                                             aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                               >  1:1326830/1‑91 (MQ=255)
                                                             aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                               <  1:3041173/91‑1 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              >  2:810161/1‑91 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              >  1:505831/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTCCAATTTGAAGTGAGTAACATATCAAATCTCAGGAGGTACATTa                          >  1:3961720/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACACTTTGAAGTGAGTAACATATCAAATCTAAGGAGGTCCATTa                          >  2:1679005/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          >  1:2172369/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          <  2:525529/91‑1 (MQ=255)
                                                                      aacaataAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                      >  2:1213507/5‑91 (MQ=25)
                                                                      aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTCATGTTGTACAATCTGAAGTGAGTAACATATAAAATCTAAGGCGGTACAGTATGaa                      >  2:1143619/1‑91 (MQ=255)
                                                                      aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAATTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                      <  1:3343707/91‑1 (MQ=255)
                                                                      aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                      <  2:1195055/91‑1 (MQ=255)
                                                                           taAAGTCCTTCCGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTAGTAATTtta                 >  2:3500625/1‑91 (MQ=255)
                                                                           taAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta                 >  1:142448/1‑91 (MQ=38)
                                                                           taAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta                 >  1:1370591/1‑91 (MQ=38)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                >  1:1899729/1‑91 (MQ=255)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              <  1:2398576/91‑1 (MQ=255)
                                                                               gcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTCAGGAGGTACATTATGAATTttaaatt             >  2:1283986/3‑91 (MQ=38)
                                                                               gcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt             >  1:886565/3‑91 (MQ=38)
                                                                                tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta            >  2:1955424/1‑91 (MQ=255)
                                                                                 ccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAATGAGGTACATTATGAATTttaaattaa           >  2:3419624/1‑91 (MQ=255)
                                                                                  cTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaaa          <  2:471323/91‑1 (MQ=255)
                                                                                   ttCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  1:1178714/1‑91 (MQ=255)
                                                                                   ttCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  2:1603411/1‑91 (MQ=255)
                                                                                     cAGTGTCGCTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATAGCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       <  2:1345151/91‑1 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATat       >  2:3932044/1‑90 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       >  1:276750/1‑91 (MQ=255)
                                                                                      aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      >  1:264808/1‑91 (MQ=255)
                                                                                      aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      <  2:2713986/91‑1 (MQ=255)
                                                                                      aGTGACACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      >  1:621578/1‑91 (MQ=255)
                                                                                      aGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      >  1:1133678/1‑91 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     >  1:2880423/1‑91 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     >  1:2941015/1‑91 (MQ=255)
                                                                                         gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca   <  2:2849664/91‑1 (MQ=255)
                                                                                          tCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacacac  >  2:2230113/1‑91 (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: