Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,754,728 A→G 5.4% intergenic (‑135/+129) 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,7280AG 5.4% 364.0 / 8.1 110intergenic (‑135/+129)gloB/opuD_3Hydroxyacylglutathione hydrolase/Glycine betaine transporter OpuD
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base G (3/3);  total (59/51)
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.74e-01

Only 100 of 113 total aligned reads displayed.

TTTATTTAATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTT  >  USA300TCH1516_ALE/2754640‑2754818
                                                                                        |                                                                                          
tttatttaATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAAt                                                                                          >  1:2913941/1‑91 (MQ=255)
tttatttaATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAAt                                                                                          >  1:3120762/1‑91 (MQ=255)
  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)
      taatttaatAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATc                                                                                    >  1:2645267/1‑91 (MQ=255)
       aatttaatAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAGATTTAAttt                                                                                   <  1:25187/91‑3 (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:4007548/1‑91 (MQ=255)
                   aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                       <  1:748322/91‑1 (MQ=255)
                   aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                       <  1:864167/91‑1 (MQ=255)
                   aaTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACNTGGAGCACTAAAAAATTAAATCTTTTAGTacaaca                                                                       <  1:2471125/91‑1 (MQ=255)
                    aTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACAt                                                                      <  1:1852943/91‑1 (MQ=255)
                    aTGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACAt                                                                      >  1:2379007/1‑91 (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:2584288/1‑91 (MQ=255)
                            cTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaaa                                                              >  1:2803161/1‑91 (MQ=255)
                             tttCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGaaaa                                                             <  1:878823/91‑1 (MQ=255)
                               tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                           <  1:512679/91‑1 (MQ=255)
                               tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                           <  1:2251241/91‑1 (MQ=255)
                               tCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAAcc                                                           <  1:1816078/91‑1 (MQ=255)
                                  aaCTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTaa                                                        <  1:3002206/91‑1 (MQ=255)
                                      tCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGt                                                    >  1:787468/1‑91 (MQ=255)
                                           ttAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATc                                               <  1:1647985/91‑1 (MQ=255)
                                           ttAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATc                                               >  1:2734664/1‑91 (MQ=255)
                                            tAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATcc                                              >  1:3058179/1‑91 (MQ=255)
                                            tAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATcc                                              >  1:3058259/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)
                                                   tGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATg                                       >  1:2350931/1‑91 (MQ=255)
                                                   tGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATg                                       >  1:446579/1‑91 (MQ=255)
                                                    gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                      >  1:3472324/1‑91 (MQ=255)
                                                    gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                      >  1:2032624/1‑91 (MQ=255)
                                                    gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                      >  1:2035241/1‑91 (MQ=255)
                                                    gCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGc                                      >  1:3473458/1‑91 (MQ=255)
                                                       tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                   <  1:3919973/91‑1 (MQ=255)
                                                       tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                   <  1:3893240/91‑1 (MQ=255)
                                                       tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                   <  1:2883436/91‑1 (MQ=255)
                                                       tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                   <  1:90501/91‑1 (MQ=255)
                                                       tCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctac                                   <  1:2799094/91‑1 (MQ=255)
                                                        cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                  <  1:1443683/91‑1 (MQ=255)
                                                        cATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGctact                                  <  1:1442422/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:2934524/1‑91 (MQ=255)
                                                              ataaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTaaaaaa                            >  1:4093372/1‑91 (MQ=255)
                                                                aaataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAg                          <  1:2228011/91‑1 (MQ=255)
                                                                 aataaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAgg                         >  1:1854232/1‑91 (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:2220156/1‑91 (MQ=255)
                                                                    aaGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGa                      <  1:4110872/91‑1 (MQ=255)
                                                                     aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                     <  1:956017/91‑1 (MQ=255)
                                                                     aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                     >  1:3371042/1‑91 (MQ=255)
                                                                     aGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                     >  1:1686132/1‑91 (MQ=255)
                                                                     aGTACTTGCAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAACAGTTCATCCTACAATGCTACTAAAAAAAGGGGaa                     >  1:1239757/1‑91 (MQ=255)
                                                                      gTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAAt                    <  1:2797510/91‑1 (MQ=255)
                                                                      gTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAAt                    >  1:2840920/1‑91 (MQ=255)
                                                                      gTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAAt                    >  1:2841470/1‑91 (MQ=255)
                                                                        aCTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATgg                  >  1:3382318/1‑91 (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)
                                                                               gCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaa           >  1:2390094/1‑91 (MQ=255)
                                                                                taCTAAAAGATTTAATTTTTTAGTGCAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa          >  1:4030201/2‑91 (MQ=38)
                                                                                taCTAAAAGATTTAATTTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa          >  1:1254927/2‑91 (MQ=39)
                                                                                cACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa          <  1:4038991/91‑1 (MQ=255)
                                                                                cACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGaaa          <  1:2335227/91‑1 (MQ=255)
                                                                                 attAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAAt         <  1:831229/89‑1 (MQ=255)
                                                                                 aCTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAAt         <  1:3915553/91‑1 (MQ=255)
                                                                                  cTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATg        >  1:1651854/1‑91 (MQ=255)
                                                                                  cTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATg        <  1:2654696/91‑1 (MQ=255)
                                                                                      aaaaTTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGatat    <  1:1360714/91‑1 (MQ=255)
                                                                                       aaaTTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATAtt   >  1:1333062/1‑91 (MQ=255)
                                                                                        aaTTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATAttt  <  1:1523401/91‑1 (MQ=255)
                                                                                        |                                                                                          
TTTATTTAATTTAATAATCAATGAGAGTCTTTCAAACTTCGATTTAGAATATGCTTCATACAATAAATAAGTACTTGGAGCACTAAAAAATTAAATCTTTTAGTACAACATTAATGAAAACCTAAAAGTTCATCCTACAATGCTACTAAAAAAAGGGGAATGGAACAGAAATGATATTT  >  USA300TCH1516_ALE/2754640‑2754818

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: