Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,158 T→A 12.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,1580TA12.1% 439.1 / 41.8 139intergenic (+65/‑66)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (60/62);  new base A (9/8);  total (69/70)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248
                                                                                          |                                                                                          
tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                            >  1:174632/1‑91 (MQ=255)
 aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                           >  1:3835466/1‑91 (MQ=255)
  gCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                          >  1:626136/1‑91 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGacac                                                                                         <  1:922524/91‑1 (MQ=255)
    aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACt                                                                                        <  1:1725624/91‑1 (MQ=255)
    aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGACACt                                                                                        <  1:139102/91‑1 (MQ=255)
     ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTg                                                                                       >  1:182460/1‑91 (MQ=255)
      tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                                      >  1:2817560/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  1:1933731/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATATATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  1:3868814/1‑91 (MQ=255)
        aaaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                                    <  1:1286519/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:2992960/91‑1 (MQ=255)
          aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                                  <  1:4129980/91‑2 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAggac                                                                                 <  1:655312/91‑3 (MQ=255)
             aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                               <  1:2406645/91‑1 (MQ=255)
             aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                               >  1:1759527/1‑91 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              <  1:819151/91‑1 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:1263702/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:1289911/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:554269/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:4030958/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:2778795/1‑91 (MQ=255)
                  aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttg                                                                          >  1:3029695/1‑91 (MQ=255)
                  aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACTTTAtaa                                                                          <  1:4013448/91‑3 (MQ=37)
                     tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                       <  1:1832324/91‑1 (MQ=255)
                      cTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                                      >  1:493963/1‑91 (MQ=255)
                      cTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                                      >  1:493565/1‑91 (MQ=255)
                       tAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTa                                                                     <  1:2464641/91‑1 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    <  1:123352/91‑1 (MQ=255)
                         ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                   >  1:562041/1‑91 (MQ=255)
                           cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                 <  1:2821139/91‑1 (MQ=255)
                           cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                 <  1:2821316/91‑1 (MQ=255)
                             ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                               >  1:3371320/1‑91 (MQ=255)
                             ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                               >  1:670532/1‑91 (MQ=255)
                             ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTGATAAAGTCGTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                               <  1:3316052/91‑1 (MQ=255)
                             ttAAATAAAGTTTCAGCTAAACACATAAAATAAAATTATTCATTGATAAAGTCGTTCGGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                               <  1:3316070/91‑1 (MQ=255)
                               aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                             <  1:4024728/91‑1 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            <  1:3019526/91‑1 (MQ=255)
                                aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            <  1:2010426/91‑1 (MQ=255)
                                   aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                         >  1:589132/1‑91 (MQ=255)
                                       tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                                     <  1:1359371/91‑1 (MQ=255)
                                       tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                                     <  1:3642033/91‑1 (MQ=255)
                                         tCAGCTAAACACATAAACTTAAATTATTCATTTATAAAGTCCTTCAGTGTCCCTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                   <  1:849827/91‑1 (MQ=255)
                                         tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                   <  1:3571394/91‑1 (MQ=255)
                                         tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                   >  1:164501/1‑91 (MQ=255)
                                          cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                  >  1:487809/1‑91 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  1:2618500/91‑1 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  1:3906127/91‑1 (MQ=255)
                                              tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                              <  1:3167538/91‑1 (MQ=255)
                                              tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                              <  1:3168045/91‑1 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            >  1:734800/1‑91 (MQ=255)
                                                 acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                           >  1:1516722/1‑91 (MQ=255)
                                                 acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                           <  1:2727726/91‑1 (MQ=255)
                                                      tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                      <  1:2997095/91‑1 (MQ=255)
                                                      tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                      <  1:403658/91‑1 (MQ=255)
                                                       aaaaTTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAag                                     >  1:2346368/1‑91 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    >  1:3661153/1‑91 (MQ=255)
                                                        aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                    <  1:714133/91‑1 (MQ=255)
                                                         aattaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagga                                   <  1:3528768/91‑1 (MQ=255)
                                                          attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                  <  1:3243393/91‑1 (MQ=255)
                                                          attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                                  <  1:684625/91‑1 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              <  1:120228/91‑1 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              >  1:3922609/1‑91 (MQ=255)
                                                              aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                              <  1:822949/91‑1 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            <  1:717952/91‑1 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            <  1:1396910/91‑1 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            >  1:1823786/1‑91 (MQ=255)
                                                                 tattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAtt                           <  1:3235505/91‑1 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          <  1:1113093/91‑1 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          <  1:2101634/91‑1 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          <  1:2102628/91‑1 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          >  1:1120476/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          >  1:3642780/1‑91 (MQ=255)
                                                                  attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                          <  1:1088817/91‑1 (MQ=255)
                                                                       tttATAAAGTCCTTCAGTGTCAATGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATAAAATCTAAGGAGGTACATTATGAAt                     >  1:1830566/1‑91 (MQ=255)
                                                                        ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtt                    <  1:3006547/91‑1 (MQ=255)
                                                                         tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   <  1:2473091/91‑1 (MQ=255)
                                                                         tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   <  1:3041562/91‑1 (MQ=255)
                                                                         tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   <  1:3041362/91‑1 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                >  1:162476/1‑91 (MQ=255)
                                                                            aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                >  1:3465684/1‑91 (MQ=255)
                                                                            aaaGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa                >  1:2555264/1‑91 (MQ=38)
                                                                             aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa               <  1:273624/91‑1 (MQ=255)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              <  1:2935972/91‑1 (MQ=255)
                                                                               gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt             >  1:3620089/1‑91 (MQ=255)
                                                                                tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta            >  1:3128357/1‑91 (MQ=255)
                                                                                tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta            >  1:3129359/1‑91 (MQ=255)
                                                                                cccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta            >  1:2904171/2‑91 (MQ=255)
                                                                                 ccTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa           <  1:1271881/91‑1 (MQ=255)
                                                                                 ccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaa           >  1:1810448/1‑91 (MQ=255)
                                                                                  cTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaaa          <  1:3162701/91‑1 (MQ=255)
                                                                                   ttCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  1:1374427/1‑91 (MQ=255)
                                                                                   ttCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  1:1867538/1‑91 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       >  1:246832/1‑91 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     <  1:2184649/91‑1 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     <  1:2184603/91‑1 (MQ=255)
                                                                                       gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAaca     <  1:189685/91‑1 (MQ=255)
                                                                                        tgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacac    <  1:312816/91‑1 (MQ=255)
                                                                                         gtCACTGAAGGGCTTTATTGTTGTACAATCTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca   >  1:1497655/1‑91 (MQ=255)
                                                                                          tCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacacac  <  1:2236363/91‑1 (MQ=255)
                                                                                          |                                                                                          
TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248

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