Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 598,501 G→A 13.3% intergenic (+50/+232) tuf → / ← yxeP_2 Elongation factor Tu/putative hydrolase YxeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE598,5010GA13.3% 302.3 / 28.3 105intergenic (+50/+232)tuf/yxeP_2Elongation factor Tu/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base G (46/45);  new base A (5/9);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGAT  >  USA300TCH1516_ALE/598411‑598591
                                                                                          |                                                                                           
gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                             >  1:851052/1‑91 (MQ=255)
  aCTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTc                                                                                           <  1:1820540/91‑1 (MQ=255)
   cTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCt                                                                                          >  1:2603588/1‑91 (MQ=255)
     gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCttt                                                                                       <  1:2160371/91‑1 (MQ=255)
     gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                                       <  1:2572770/91‑3 (MQ=25)
      tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                                       <  1:24767/91‑3 (MQ=25)
      tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                                       <  1:700312/91‑3 (MQ=25)
      tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                                       <  1:700386/91‑3 (MQ=25)
        ggATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtttttt                                                                                    >  1:444288/1‑91 (MQ=255)
        ggATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAGGGATCttcttt                                                                                    <  1:1703377/91‑1 (MQ=37)
        agATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtttttt                                                                                    >  1:2048316/2‑91 (MQ=255)
         gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCttttttt                                                                                   <  1:3242525/91‑1 (MQ=255)
         gATCAGGCGTTGTTACTGAAATCATTAAATAATGTCTAATTTCTTAGATTTTATATAAAGAGAAGATCCCTCAATCGAGGGGTCtttttta                                                                                   <  1:3685827/91‑2 (MQ=255)
           tCAGGCGTTGTTACTGAAATCATTATATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtttttttta                                                                                 <  1:3848303/91‑3 (MQ=255)
            cAGGCGTTGTTACTGAAATCATTAAAGAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTAta                                                                               <  1:2547536/91‑2 (MQ=18)
             aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtg                                                                               >  1:14042/1‑91 (MQ=255)
             aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtg                                                                               <  1:3880914/91‑1 (MQ=255)
             aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatat                                                                               <  1:3070988/91‑4 (MQ=18)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:2221904/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:1664921/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:2966348/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:2221903/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:213168/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:3810069/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                              <  1:1493879/91‑5 (MQ=17)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                              <  1:1493719/91‑5 (MQ=17)
               gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                             <  1:2527328/91‑1 (MQ=255)
                 gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                           <  1:2687307/91‑1 (MQ=255)
                 gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                           <  1:1923416/91‑1 (MQ=255)
                 gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                           <  1:2394658/91‑1 (MQ=255)
                 gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                           <  1:2687369/91‑1 (MQ=255)
                  ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                                          >  1:3707638/1‑91 (MQ=255)
                  ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                                          >  1:2948596/1‑91 (MQ=255)
                     ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                                       <  1:3380831/91‑1 (MQ=255)
                     ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                                       >  1:4127264/1‑91 (MQ=255)
                     ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                                       <  1:696432/91‑1 (MQ=255)
                        cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTg                                                                    >  1:2161660/1‑91 (MQ=255)
                           aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                                 >  1:878740/1‑91 (MQ=255)
                           aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                                 >  1:878581/1‑91 (MQ=255)
                             aTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATg                                                               <  1:3139103/91‑1 (MQ=255)
                               cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGc                                                             <  1:1257629/91‑1 (MQ=255)
                                  tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                                          >  1:1321887/1‑91 (MQ=255)
                                  tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                                          <  1:2733160/91‑1 (MQ=255)
                                  tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                                          <  1:4111331/91‑1 (MQ=255)
                                      taatttctaatttcTTAGATTTTATATAAAAAGAAGATCTCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                                      <  1:2973122/91‑1 (MQ=255)
                                      taatttctaatttcTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                                      >  1:1978576/1‑91 (MQ=255)
                                       aatttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAt                                                     <  1:702355/91‑1 (MQ=255)
                                           tctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTag                                                 <  1:151099/91‑1 (MQ=255)
                                              aatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaag                                              <  1:493467/91‑1 (MQ=255)
                                              aatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaag                                              <  1:493907/91‑1 (MQ=255)
                                                tttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAgaagaa                                            <  1:139600/91‑1 (MQ=255)
                                                   ctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaa                                         >  1:2550930/1‑91 (MQ=255)
                                                   ctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaa                                         >  1:335575/1‑91 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      <  1:4053617/91‑1 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      <  1:2673797/91‑1 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      >  1:1976701/1‑91 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      >  1:2425300/1‑91 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      >  1:1095114/1‑91 (MQ=255)
                                                        aTTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAAtt                                    <  1:4032462/91‑1 (MQ=255)
                                                        aTTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAAtt                                    >  1:3264751/1‑91 (MQ=255)
                                                          tttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                                  <  1:2223108/91‑1 (MQ=255)
                                                          tttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                                  <  1:1057732/91‑1 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  1:4112016/1‑91 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  1:1809949/1‑91 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               <  1:375750/91‑1 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  1:1053856/1‑91 (MQ=255)
                                                                taAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAg                            >  1:1043604/1‑91 (MQ=255)
                                                                 aaaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                           >  1:3427985/1‑91 (MQ=255)
                                                                  aaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                          >  1:3657156/1‑91 (MQ=255)
                                                                  aaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                          >  1:3656818/1‑91 (MQ=255)
                                                                   aaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgact                         >  1:3157166/1‑91 (MQ=255)
                                                                     agaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                       <  1:2234235/91‑1 (MQ=255)
                                                                     aaaagaCCCCTCGATTGAGGGATCTTCTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactg                        >  1:2478825/3‑91 (MQ=25)
                                                                      gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac                      >  1:430875/1‑91 (MQ=255)
                                                                      aaagaCCCCTCGATTGAGGGATCTTCTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                       >  1:1787072/2‑91 (MQ=25)
                                                                        agaccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATAGATGAAAGACTGACta                    >  1:3256999/5‑91 (MQ=38)
                                                                        agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    >  1:4120704/1‑91 (MQ=255)
                                                                         gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa                   >  1:3671001/1‑91 (MQ=255)
                                                                         gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa                   >  1:2469059/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:505748/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  <  1:2928760/91‑1 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:2738234/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:3585481/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTCAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:2737849/1‑91 (MQ=255)
                                                                           tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata                 >  1:3426965/1‑91 (MQ=255)
                                                                            cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa                >  1:458035/1‑91 (MQ=255)
                                                                            cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa                <  1:1152823/91‑1 (MQ=255)
                                                                             ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa               <  1:1439350/91‑1 (MQ=255)
                                                                              cTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaaa              <  1:3251672/91‑1 (MQ=255)
                                                                               tCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc             >  1:634866/1‑91 (MQ=255)
                                                                                   tCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTtat         <  1:1942909/91‑1 (MQ=255)
                                                                                   tCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTtat         <  1:1941826/91‑1 (MQ=255)
                                                                                    cGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTtata        >  1:3717684/1‑91 (MQ=255)
                                                                                    cGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTtata        >  1:3717810/1‑91 (MQ=255)
                                                                                     gAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAGTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATaa       <  1:3290043/91‑1 (MQ=255)
                                                                                       ggggTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACt     >  1:1652407/1‑91 (MQ=255)
                                                                                       ggggTCTTGTTTTAATGTGTAAATTTTGTAATGGCGATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACt     >  1:1236004/1‑91 (MQ=255)
                                                                                          gTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGat  >  1:426159/1‑91 (MQ=255)
                                                                                          gTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGat  <  1:510849/91‑1 (MQ=255)
                                                                                          gTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGat  <  1:1033566/91‑1 (MQ=255)
                                                                                          |                                                                                           
GTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGAT  >  USA300TCH1516_ALE/598411‑598591

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: