Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,754,732 A→T 5.7% intergenic (‑139/+125) gloB ← / ← opuD_3 Hydroxyacylglutathione hydrolase/Glycine betaine transporter OpuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,754,7320AT 5.7% 375.0 / 10.0 106intergenic (‑139/+125)gloB/opuD_3Hydroxyacylglutathione hydrolase/Glycine betaine transporter OpuD
Reads supporting (aligned to +/- strand):  ref base A (53/47);  new base T (3/3);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 113 total aligned reads displayed.

TATTTAATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTTTCAC  >  USA300TCH1516_ALE/2754642‑2754822
                                                                                          |                                                                                          
tatttaATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTa                                                                                            >  1:3675186/1‑91 (MQ=255)
tatttaATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTa                                                                                            >  1:1593937/1‑91 (MQ=255)
 atttaatttaATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTaa                                                                                           <  1:793697/91‑1 (MQ=255)
  tttaatttaaTAATCAATGAGAGTCTTTCAAACTTCGATTTAGATTATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTaaa                                                                                          <  1:3299633/91‑1 (MQ=255)
  tttaatttaaTAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTaaa                                                                                          <  1:3300559/91‑1 (MQ=255)
     aatttaatAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAGATTTAAttt                                                                                       <  1:25187/91‑3 (MQ=255)
      atttaatAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCtt                                                                                      <  1:1903095/91‑1 (MQ=255)
         taataaTCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAGATTTAATTTTTTa                                                                                   <  1:2420842/91‑1 (MQ=255)
         taataaTCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAGATTTAATTTTTTa                                                                                   <  1:3104913/91‑1 (MQ=255)
            taatCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACGATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTactg                                                                                >  1:1300312/1‑88 (MQ=255)
            taatCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTa                                                                                <  1:967139/91‑1 (MQ=255)
             aatCAGTGAGAGTCTTTCAAACTTCTATTTAGAATATGCGTCATACAATAAATCAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTac                                                                               <  1:3856990/91‑1 (MQ=255)
               tCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaa                                                                             >  1:2129839/1‑91 (MQ=255)
                cAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaac                                                                            >  1:3614663/1‑91 (MQ=255)
                 aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                           <  1:748322/91‑1 (MQ=255)
                 aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                           <  1:864167/91‑1 (MQ=255)
                 aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                           >  1:4007548/1‑91 (MQ=255)
                 aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACNTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                           <  1:2471125/91‑1 (MQ=255)
                  aTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACAt                                                                          >  1:2379007/1‑91 (MQ=255)
                  aTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACAt                                                                          <  1:1852943/91‑1 (MQ=255)
                   tGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACAtt                                                                         <  1:2304861/91‑1 (MQ=255)
                     agagTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTaa                                                                       <  1:3762475/91‑1 (MQ=255)
                      gagTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAAt                                                                      >  1:1261861/1‑91 (MQ=255)
                       agTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATg                                                                     >  1:159573/1‑91 (MQ=255)
                        gTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGa                                                                    >  1:2029741/1‑91 (MQ=255)
                         tCTTTCAAACTTCGATTTAGAATATGCTTTATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaa                                                                   <  1:3326023/91‑1 (MQ=255)
                          cTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaaa                                                                  >  1:2803161/1‑91 (MQ=255)
                          cTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaaa                                                                  >  1:2584288/1‑91 (MQ=255)
                           tttCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaaaa                                                                 <  1:878823/91‑1 (MQ=255)
                             tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                               <  1:2251241/91‑1 (MQ=255)
                             tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                               <  1:512679/91‑1 (MQ=255)
                             tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                               <  1:1816078/91‑1 (MQ=255)
                              cAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCt                                                              <  1:1412104/91‑1 (MQ=255)
                                aaCTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTaa                                                            <  1:3002206/91‑1 (MQ=255)
                                    tCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGt                                                        >  1:787468/1‑91 (MQ=255)
                                         ttAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATc                                                   >  1:2734664/1‑91 (MQ=255)
                                          tAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATcc                                                  >  1:3058259/1‑91 (MQ=255)
                                          tAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATcc                                                  >  1:3058179/1‑91 (MQ=255)
                                            gAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTa                                                >  1:3475764/1‑91 (MQ=255)
                                            gAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTa                                                >  1:1658349/1‑91 (MQ=255)
                                             aaTATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTAc                                               <  1:1656639/91‑1 (MQ=255)
                                              atatGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACa                                              <  1:1014545/91‑1 (MQ=255)
                                              atatGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACa                                              <  1:479695/91‑1 (MQ=255)
                                               tatGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACaa                                             >  1:264374/1‑91 (MQ=255)
                                                atGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAAt                                            >  1:3268709/1‑91 (MQ=255)
                                                atGCCTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTAGAAt                                            >  1:3268556/1‑91 (MQ=255)
                                                 tGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATg                                           >  1:2350931/1‑91 (MQ=255)
                                                 tGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATg                                           >  1:446579/1‑91 (MQ=255)
                                                  gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                          >  1:2035241/1‑91 (MQ=255)
                                                  gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                          >  1:3473458/1‑91 (MQ=255)
                                                  gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                          >  1:3472324/1‑91 (MQ=255)
                                                  gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                          >  1:2032624/1‑91 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       <  1:3919973/91‑1 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       <  1:2883436/91‑1 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       >  1:2906444/1‑91 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       <  1:90501/91‑1 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       <  1:3893240/91‑1 (MQ=255)
                                                     tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                       <  1:2799094/91‑1 (MQ=255)
                                                      cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                      <  1:1442422/91‑1 (MQ=255)
                                                      cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                      <  1:1443683/91‑1 (MQ=255)
                                                      cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                      <  1:862236/91‑1 (MQ=255)
                                                      cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                      >  1:2464404/1‑91 (MQ=255)
                                                       aTACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctacta                                     >  1:2602812/1‑91 (MQ=255)
                                                         aCAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCTTACAATGCTACTaaa                                   >  1:477691/1‑91 (MQ=255)
                                                           aataaataAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTaaaaa                                 <  1:2106150/91‑1 (MQ=255)
                                                           aagaaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTaaaaa                                 <  1:2106480/88‑1 (MQ=255)
                                                            ataaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTaaaaaa                                >  1:4093372/1‑91 (MQ=255)
                                                            ataaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTaaaaaa                                >  1:2934524/1‑91 (MQ=255)
                                                              aaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAg                              <  1:2228011/91‑1 (MQ=255)
                                                              aaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAg                              <  1:2228001/91‑1 (MQ=255)
                                                               aataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAgg                             <  1:63201/91‑1 (MQ=255)
                                                                 taaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAgggg                           <  1:3177222/91‑1 (MQ=255)
                                                                 taaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAgggg                           <  1:1420686/91‑1 (MQ=255)
                                                                  aaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGa                          >  1:2733709/1‑91 (MQ=255)
                                                                  aaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGa                          <  1:4110872/91‑1 (MQ=255)
                                                                   aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                         >  1:3371042/1‑91 (MQ=255)
                                                                   aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                         <  1:956017/91‑1 (MQ=255)
                                                                   aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                         >  1:1686132/1‑91 (MQ=255)
                                                                   aGTACTTGCAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAACAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                         >  1:1239757/1‑91 (MQ=255)
                                                                    gTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAAt                        >  1:2840920/1‑91 (MQ=255)
                                                                    gTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAAt                        <  1:2797510/91‑1 (MQ=255)
                                                                      aCTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATgg                      >  1:3382268/1‑91 (MQ=255)
                                                                       cTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGa                     >  1:1325754/1‑91 (MQ=255)
                                                                         tGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAAc                   >  1:1628524/1‑91 (MQ=255)
                                                                           gAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAg                 >  1:174111/1‑91 (MQ=255)
                                                                             gtaCTAAAAGATTTAATTTTTTAATGCAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaa               >  1:3287946/3‑91 (MQ=38)
                                                                             gCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaa               >  1:232746/1‑91 (MQ=255)
                                                                              taCTAAAAGATTTAATTTTTTAGTGCAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa              >  1:4030201/2‑91 (MQ=38)
                                                                              taCTAAAAGATTTAATTTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa              >  1:1254927/2‑91 (MQ=39)
                                                                              cACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa              <  1:2335227/91‑1 (MQ=255)
                                                                              cACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa              <  1:4038991/91‑1 (MQ=255)
                                                                               attAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAAt             <  1:831229/89‑1 (MQ=255)
                                                                               aCTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAAt             <  1:3915553/91‑1 (MQ=255)
                                                                                cTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATg            <  1:2654696/91‑1 (MQ=255)
                                                                                cTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATg            >  1:1651854/1‑91 (MQ=255)
                                                                                cTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATg            >  1:1651729/1‑91 (MQ=255)
                                                                                     aaaTTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATAtt       >  1:1333062/1‑91 (MQ=255)
                                                                                      aaTTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATAttt      <  1:1523401/91‑1 (MQ=255)
                                                                                          aaaTCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTTTcac  >  1:2413210/1‑91 (MQ=255)
                                                                                          aaaTCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTTTcac  >  1:1292010/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATTTAATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTTTCAC  >  USA300TCH1516_ALE/2754642‑2754822

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: