Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 598,501 G→A 8.9% 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,5010GA 8.9% 728.6 / 41.3 236intergenic (+50/+232)tuf/yxeP_2Elongation factor Tu/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base G (117/97);  new base A (8/13);  total (126/110)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 251 total aligned reads displayed.

GTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGAT  >  USA300TCH1516_ALE/598411‑598591
                                                                                          |                                                                                           
gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                             <  1:1118163/91‑1 (MQ=255)
gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                             <  1:1280133/91‑1 (MQ=255)
gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                             <  2:221409/91‑1 (MQ=255)
   cTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCt                                                                                          <  2:1509605/91‑1 (MQ=255)
    tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtt                                                                                         <  2:665289/91‑1 (MQ=255)
     gTAGGCTCAGGCGTTGTTACTGAACTCAATAAATAATTTCTAATTTCATAGATTTTATATAAAAAGAAGATACCTCAATCGAGGGTTCttt                                                                                       >  2:3451199/1‑91 (MQ=255)
     gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAGGGATCttc                                                                                        <  1:627314/91‑2 (MQ=38)
      tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCtttt                                                                                      >  1:3648582/1‑91 (MQ=255)
         gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCttttttt                                                                                   <  1:2260547/91‑1 (MQ=255)
         gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCttttttt                                                                                   >  1:2919038/1‑91 (MQ=255)
         gATCAGGCGTTGTTACTGAAATCATTAAAGAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCttctttt                                                                                   <  1:3173383/91‑1 (MQ=25)
           tCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTaa                                                                                 >  2:323413/1‑91 (MQ=255)
            cAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAt                                                                                >  1:1353313/1‑91 (MQ=255)
            cAGGCGTTGTTACTGAAATAATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTata                                                                                <  2:1122329/91‑3 (MQ=18)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              >  1:3887238/1‑91 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              <  1:1517337/91‑1 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgt                                                                              >  2:2354113/1‑91 (MQ=255)
              ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAGGGATCTTCTTTTTatat                                                                               <  1:2420037/91‑4 (MQ=21)
               gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                             >  1:212023/1‑91 (MQ=255)
               gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAtgtg                                                                             >  2:2741262/1‑91 (MQ=255)
               gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                              <  1:873807/91‑5 (MQ=18)
               gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATCTTCTTTTTatata                                                                              <  1:2281300/91‑5 (MQ=18)
                  ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaa                                                                          <  2:2685199/91‑1 (MQ=255)
                   tgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTaaa                                                                         <  2:267798/91‑1 (MQ=255)
                      tACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAttt                                                                      <  2:3405993/91‑1 (MQ=255)
                       aCTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAAtttt                                                                     >  2:3294694/1‑91 (MQ=255)
                         tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGt                                                                   <  1:2978291/91‑1 (MQ=255)
                         tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGt                                                                   <  1:1679280/91‑1 (MQ=255)
                             aTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATg                                                               <  2:293974/91‑1 (MQ=255)
                              tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATgg                                                              >  2:1610347/1‑91 (MQ=255)
                              tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATgg                                                              >  2:1610268/1‑91 (MQ=255)
                               cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGc                                                             <  2:1353357/91‑1 (MQ=255)
                               cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGc                                                             <  2:3722693/91‑1 (MQ=255)
                                aTTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAAAGTGTAAATTTTGTAATGGCt                                                            >  2:2498247/1‑91 (MQ=255)
                                 ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTa                                                           >  1:323258/1‑91 (MQ=255)
                                 ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTa                                                           >  2:3630308/1‑91 (MQ=255)
                                    aataatTTCTCATTTCTTAGATTTTATATAATAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTc                                                        <  2:2605280/91‑1 (MQ=255)
                                     ataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCg                                                       >  1:1694852/1‑91 (MQ=255)
                                      taatttctaatttcTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                                      >  1:839042/1‑91 (MQ=255)
                                       aatttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCgtt                                                     >  2:2624294/1‑89 (MQ=255)
                                       aatttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAt                                                     <  2:3560436/91‑1 (MQ=255)
                                       aatttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAt                                                     <  1:766928/91‑1 (MQ=255)
                                         tttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGAttt                                                   <  1:483929/91‑1 (MQ=255)
                                            ctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTaga                                                >  2:838982/1‑91 (MQ=255)
                                             taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                               >  2:401654/1‑91 (MQ=255)
                                              aatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaag                                              <  2:1410951/91‑1 (MQ=255)
                                                tttctTAGATTTTATATAAAAAGAAGGTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTGAATTTTGTAGTGGCTATTCGATTTAgaagaa                                            <  2:1464062/91‑1 (MQ=255)
                                                 ttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAAc                                           >  1:1603783/1‑91 (MQ=255)
                                                  tctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACa                                          >  1:3625353/1‑91 (MQ=255)
                                                    ttAGATTTTATATAAAAAGAAGATCCCTCAATTGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                        >  2:2868330/1‑91 (MQ=255)
                                                    ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                        <  1:417971/91‑1 (MQ=255)
                                                     tAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaata                                       >  2:646625/1‑91 (MQ=255)
                                                      aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                                      <  2:3076676/91‑1 (MQ=255)
                                                       gcttttATATAAACAGAAGCCCCCTCAACCGAGGGGTCCTTTTTTAATGTGTAAATTTTGTAAGGGCTATTCGATTTAGACGAACaataat                                     >  2:3651771/3‑91 (MQ=255)
                                                         ttttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtg                                   <  1:3006960/91‑1 (MQ=255)
                                                          tttATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtga                                  <  2:1298545/91‑1 (MQ=255)
                                                            tatatgAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                                <  2:512565/91‑1 (MQ=255)
                                                            tatataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                                <  1:369813/91‑1 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  2:78372/1‑91 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  2:1223538/1‑91 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  1:1870711/1‑91 (MQ=255)
                                                             atataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                               >  2:1210148/1‑91 (MQ=255)
                                                              tataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaa                              >  2:2038095/1‑91 (MQ=255)
                                                               ataAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaaa                             >  1:1018340/1‑91 (MQ=255)
                                                                 aaaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                           <  2:1663539/91‑1 (MQ=255)
                                                                 aaaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                           <  2:2484053/91‑1 (MQ=255)
                                                                  aaaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                          <  1:2589366/91‑1 (MQ=255)
                                                                   aaaGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgact                         >  2:3502222/1‑91 (MQ=255)
                                                                    aagaagATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactg                        <  1:9973/91‑1 (MQ=255)
                                                                     agaagaTCCCTCAATCGAGGGGTCTTTTTTTACTGTGTCCATTTTGTACTGGCTATTCGCTTTCGAAGAACACTAATTGACGAACgactga                       >  2:456300/1‑91 (MQ=255)
                                                                     agaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                       <  1:2702459/91‑1 (MQ=255)
                                                                      gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac                      >  2:743794/1‑91 (MQ=255)
                                                                      gaagaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactgac                      <  2:1847744/91‑1 (MQ=255)
                                                                        agaccccTCGATTGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    >  1:714759/5‑91 (MQ=255)
                                                                        agaccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    >  1:3811275/5‑91 (MQ=38)
                                                                        agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    >  1:1816542/1‑91 (MQ=255)
                                                                        agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    >  1:139557/1‑91 (MQ=255)
                                                                        agaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta                    <  2:3095922/91‑1 (MQ=255)
                                                                         gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa                   >  1:1188679/1‑91 (MQ=255)
                                                                         gaTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGACTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa                   >  2:3314615/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:3734842/1‑91 (MQ=255)
                                                                          aTCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat                  >  1:3524402/1‑91 (MQ=255)
                                                                           tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata                 >  1:3385766/1‑91 (MQ=255)
                                                                           tCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATCCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata                 >  2:2895903/1‑91 (MQ=255)
                                                                           ccccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata                 >  2:1733769/2‑91 (MQ=255)
                                                                            cccTCGATTGAGGGATCTTTTTTTAATGTGTAAATTTCGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa                >  1:2155395/1‑91 (MQ=255)
                                                                            cccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa                >  2:2080624/1‑91 (MQ=255)
                                                                             ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa               >  2:2267042/1‑91 (MQ=255)
                                                                             ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa               <  2:392818/91‑1 (MQ=255)
                                                                             ccTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa               >  1:279827/1‑91 (MQ=255)
                                                                                cAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACt            >  1:3346782/1‑91 (MQ=255)
                                                                                 aaTCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACtt           <  1:2056537/91‑1 (MQ=255)
                                                                                  aTCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTta          >  2:3012935/1‑91 (MQ=255)
                                                                                     gAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATaa       <  2:2431538/91‑1 (MQ=255)
                                                                                     gAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATaa       <  2:2724516/91‑1 (MQ=255)
                                                                                     gAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATaa       >  1:2461925/1‑91 (MQ=255)
                                                                                     gAGGGGTCTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAAc      <  2:1667447/91‑1 (MQ=255)
                                                                                       ggggTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACt     >  1:579092/1‑91 (MQ=255)
                                                                                        gggTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTg    >  1:3593707/1‑91 (MQ=255)
                                                                                          gTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGat  <  1:1427579/91‑1 (MQ=255)
                                                                                          |                                                                                           
GTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTTTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTATAACTGAT  >  USA300TCH1516_ALE/598411‑598591

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: