Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 598,492 A→G 14.5% intergenic (+41/+241) 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,4920AG14.5% 296.2 / 36.5 110intergenic (+41/+241)tuf/yxeP_2Elongation factor Tu/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base A (48/46);  new base G (5/11);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 114 total aligned reads displayed.

AAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAC  >  USA300TCH1516_ALE/598402‑598580
                                                                                          |                                                                                         
aaGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCa                                                                                           >  1:2486779/1‑91 (MQ=255)
 aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCaa                                                                                          >  1:3609374/1‑91 (MQ=255)
 aGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCaa                                                                                          >  1:1983477/1‑91 (MQ=255)
   gTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATc                                                                                        <  1:3940642/91‑1 (MQ=255)
    tGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCg                                                                                       >  1:530065/1‑91 (MQ=255)
      gACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAg                                                                                     <  1:2948555/91‑1 (MQ=255)
      gACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAg                                                                                     >  1:942974/1‑91 (MQ=255)
       aCGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgg                                                                                    <  1:116252/91‑1 (MQ=255)
       aCGTACTATAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAgg                                                                                    <  1:2067866/91‑1 (MQ=255)
        cGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAggg                                                                                   >  1:1172042/1‑91 (MQ=255)
         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:700386/91‑3 (MQ=25)
               tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttct                                                                            <  1:24767/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:3880914/91‑1 (MQ=255)
                      aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtg                                                                    >  1:14042/1‑91 (MQ=255)
                      aGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatat                                                                    <  1:3070988/91‑4 (MQ=18)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                   <  1:2221903/91‑1 (MQ=255)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                   <  1:3810069/91‑1 (MQ=255)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                   <  1:2966348/91‑1 (MQ=255)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                   <  1:213168/91‑1 (MQ=255)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                   <  1:1493719/91‑5 (MQ=17)
                       ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                   <  1:1493879/91‑5 (MQ=17)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                  >  1:4132047/1‑91 (MQ=255)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                  <  1:2527328/91‑1 (MQ=255)
                        gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                  >  1:1977928/1‑91 (MQ=255)
                          gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                <  1:2687369/91‑1 (MQ=255)
                          gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                <  1:2687307/91‑1 (MQ=255)
                          gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                <  1:2394658/91‑1 (MQ=255)
                          gttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTa                                                                <  1:1923416/91‑1 (MQ=255)
                           ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                               >  1:2948596/1‑91 (MQ=255)
                           ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                               >  1:2260643/1‑91 (MQ=255)
                           ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                               <  1:1355860/91‑1 (MQ=255)
                              ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                            >  1:4127264/1‑91 (MQ=255)
                              ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                            <  1:3380831/91‑1 (MQ=255)
                              ttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtt                                                            <  1:696432/91‑1 (MQ=255)
                                 cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTg                                                         >  1:2161660/1‑91 (MQ=255)
                                    aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                      >  1:878581/1‑91 (MQ=255)
                                    aaaTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTaa                                                      >  1:878740/1‑91 (MQ=255)
                                      aTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATg                                                    <  1:3139103/91‑1 (MQ=255)
                                        cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGc                                                  <  1:1257629/91‑1 (MQ=255)
                                           tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                               <  1:4111331/91‑1 (MQ=255)
                                           tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                               <  1:2733160/91‑1 (MQ=255)
                                           tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTAt                                               >  1:3607673/1‑91 (MQ=255)
                                               taatttctaatttcTTAGATTTTATATAAAAAGAAGATCTCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                           <  1:2973122/91‑1 (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:335575/1‑91 (MQ=255)
                                                            ctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaa                              >  1:2550930/1‑91 (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:3264751/1‑91 (MQ=255)
                                                                 aTTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAAtt                         <  1:4032462/91‑1 (MQ=255)
                                                                   tttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                       <  1:1057732/91‑1 (MQ=255)
                                                                   tttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                       <  1:2223108/91‑1 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                    >  1:1809949/1‑91 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                    >  1:1053856/1‑91 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                    >  1:4112016/1‑91 (MQ=255)
                                                                      atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                    <  1:375750/91‑1 (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)
                                                                               gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac           >  1:430875/1‑91 (MQ=255)
                                                                               aaagaCCCCTCGATTGAGGGATCTTCTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga            >  1:1787072/2‑91 (MQ=25)
                                                                                aagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgact          <  1:2903160/91‑1 (MQ=255)
                                                                                 agaccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATAGATGAAAGACTGACta         >  1:3256955/5‑91 (MQ=38)
                                                                                 agaccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATAGATGAAAGACTGACta         >  1:3256999/5‑91 (MQ=38)
                                                                                 agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta         >  1:4120704/1‑91 (MQ=255)
                                                                                  gaccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa        >  1:3788539/4‑91 (MQ=255)
                                                                                  gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa        >  1:2469059/1‑91 (MQ=255)
                                                                                  gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa        >  1:3671001/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat       >  1:505748/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat       >  1:3585481/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat       >  1:2738234/1‑91 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat       <  1:2928760/91‑1 (MQ=255)
                                                                                   aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTCAGAAGAACAATAATTGATGAAAGACTGACtaat       >  1:2737849/1‑91 (MQ=255)
                                                                                    tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata      >  1:3426965/1‑91 (MQ=255)
                                                                                     cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa     <  1:1152823/91‑1 (MQ=255)
                                                                                      ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa    <  1:1439350/91‑1 (MQ=255)
                                                                                      ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa    >  1:923277/1‑91 (MQ=255)
                                                                                        tCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc  >  1:2555942/1‑91 (MQ=255)
                                                                                        tCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc  >  1:634810/1‑91 (MQ=255)
                                                                                        tCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc  >  1:634866/1‑91 (MQ=255)
                                                                                          |                                                                                         
AAGGTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAC  >  USA300TCH1516_ALE/598402‑598580

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: