Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 85,835 T→A 5.5% intergenic (+854/‑166) nikE_1 → / → copB Nickel import ATP‑binding protein NikE/Copper‑exporting P‑type ATPase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE85,8350TA 5.5% 785.7 / 23.2 234intergenic (+854/‑166)nikE_1/copBNickel import ATP‑binding protein NikE/Copper‑exporting P‑type ATPase B
Reads supporting (aligned to +/- strand):  ref base T (110/111);  new base A (9/4);  total (119/115)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 240 total aligned reads displayed.

ATATTAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCGAAA  >  USA300TCH1516_ALE/85745‑85925
                                                                                          |                                                                                          
atatTAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACAt                                                                                            <  2:2792778/91‑1 (MQ=255)
  atTAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTg                                                                                          >  1:1380237/1‑91 (MQ=255)
   ttAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGt                                                                                         >  2:35195/1‑91 (MQ=255)
   ttAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGt                                                                                         <  1:2772730/91‑1 (MQ=255)
   ttAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGt                                                                                         >  2:3218022/1‑91 (MQ=255)
   ttAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGt                                                                                         <  1:3732552/91‑1 (MQ=255)
    tAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGtt                                                                                        >  2:1162235/1‑91 (MQ=255)
    tAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGtt                                                                                        >  1:3765487/1‑91 (MQ=255)
     aaTATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGttt                                                                                       <  2:746962/91‑1 (MQ=255)
      atatatGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGtttt                                                                                      >  2:3655562/1‑91 (MQ=255)
       tatatGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTa                                                                                     >  2:2581011/1‑91 (MQ=255)
       tatatGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTa                                                                                     <  1:1873686/91‑1 (MQ=255)
        atatGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTGGATCTTTTAATTAAAACATTGTTTTaa                                                                                    >  1:2697725/1‑91 (MQ=255)
             aCTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATc                                                                               <  1:3412937/91‑1 (MQ=255)
              cTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCg                                                                              <  1:697820/91‑1 (MQ=255)
              cTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCg                                                                              >  1:2151040/1‑91 (MQ=255)
               ttGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGa                                                                             >  2:2137615/1‑91 (MQ=255)
                tgtgTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGAt                                                                            <  2:3490838/91‑1 (MQ=255)
                   gtTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGAc                                                                         >  2:1028885/1‑91 (MQ=255)
                   gtTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGAc                                                                         >  1:3294759/1‑91 (MQ=255)
                   gtTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGAc                                                                         >  2:2885092/1‑91 (MQ=255)
                   gtTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTCTTAAATATCGATGAc                                                                         >  1:2189370/1‑91 (MQ=255)
                    ttGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACa                                                                        <  1:1374349/91‑1 (MQ=255)
                      ggTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAg                                                                      <  1:3814803/91‑1 (MQ=255)
                      ggTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAg                                                                      >  1:1023992/1‑91 (MQ=255)
                      ggTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAg                                                                      >  2:1372966/1‑91 (MQ=255)
                        tAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATTGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGAGAAGGt                                                                    >  1:3614695/1‑91 (MQ=255)
                        tAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGt                                                                    >  2:640684/1‑91 (MQ=255)
                          aTGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAat                                                                                  <  1:122067/75‑1 (MQ=255)
                          aTGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCt                                                                  <  2:1210797/91‑1 (MQ=255)
                           tGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTa                                                                 >  2:1556033/1‑91 (MQ=255)
                               aaCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGt                                                             <  1:870118/91‑1 (MQ=255)
                               aaCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGt                                                             <  2:1926110/91‑1 (MQ=255)
                               aaCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGt                                                             <  2:3148602/91‑1 (MQ=255)
                                aCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTa                                                            <  2:2697098/91‑1 (MQ=255)
                                aCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTa                                                            <  1:3398271/91‑1 (MQ=255)
                                aCCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTa                                                            >  1:580703/1‑91 (MQ=255)
                                aCCTTGGTTAGAGTGTATTGTCATGTGAAAGTTAAGATTAGATCTTTTAATTAAAACATTGTTGTAAATATCGATGACAAGGTCTAATGTa                                                            <  2:3420662/91‑1 (MQ=255)
                                  cTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAgg                                                          >  1:3387978/1‑91 (MQ=255)
                                  cTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAgg                                                          <  1:1529388/91‑1 (MQ=255)
                                  cTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAgg                                                          <  2:2065496/91‑1 (MQ=255)
                                  cTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAgg                                                          <  1:2339293/91‑1 (MQ=255)
                                   ttGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACACGGTCTAACGCAggc                                                         >  2:3795945/1‑90 (MQ=255)
                                   ttGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGa                                                         <  2:1900017/91‑1 (MQ=255)
                                   ttGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGa                                                         <  1:3439701/91‑1 (MQ=255)
                                    tGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGAc                                                        >  2:3250289/1‑91 (MQ=255)
                                     ggTTAGAGTGTAATGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACg                                                       >  1:2990672/1‑91 (MQ=255)
                                     ggTTAGAGAGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACg                                                       >  1:2058334/1‑91 (MQ=255)
                                      gTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACgt                                                      >  2:495212/1‑91 (MQ=255)
                                      gTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACgt                                                      >  2:2629886/1‑91 (MQ=255)
                                        tAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACgtgt                                                    <  1:1491808/91‑1 (MQ=255)
                                           agTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGa                                                 >  1:1109848/1‑91 (MQ=255)
                                            gtgtATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGat                                                >  1:2250359/1‑91 (MQ=255)
                                            gtgtATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTGTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGat                                                <  2:903367/91‑1 (MQ=255)
                                              gtATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGatat                                              <  1:3188277/91‑1 (MQ=255)
                                              gtATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGatat                                              >  2:577900/1‑91 (MQ=255)
                                               tATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGAtata                                             >  2:2814081/1‑91 (MQ=255)
                                                aTTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATaa                                            >  2:126458/1‑91 (MQ=255)
                                                aTTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATaa                                            <  1:1463715/91‑1 (MQ=255)
                                                   gTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGt                                         >  2:3877315/1‑91 (MQ=255)
                                                   gTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGt                                         <  1:904228/91‑1 (MQ=255)
                                                   gTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGt                                         <  1:638825/91‑1 (MQ=255)
                                                    tCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGtt                                        >  1:1061780/1‑91 (MQ=255)
                                                    tCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGtt                                        >  2:332052/1‑91 (MQ=255)
                                                    tCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGtt                                        <  2:1469399/91‑1 (MQ=255)
                                                    tCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGtt                                        <  2:1442210/91‑1 (MQ=255)
                                                     cATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTc                                       >  2:1151330/1‑91 (MQ=255)
                                                      aTGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCg                                      >  1:2417345/1‑91 (MQ=255)
                                                      aTGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCg                                      <  1:2599408/91‑1 (MQ=255)
                                                        gtgATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGa                                    >  1:776561/1‑91 (MQ=255)
                                                          gATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAAt                                  <  1:332052/91‑1 (MQ=255)
                                                           aTAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAAtt                                 >  2:407711/1‑91 (MQ=255)
                                                           aTAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAAtt                                 >  1:2824242/1‑91 (MQ=255)
                                                            tAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATta                                >  2:2316663/1‑91 (MQ=255)
                                                             aGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAg                                                  >  1:2320322/1‑72 (MQ=255)
                                                                  aGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATaaa                          >  1:1550357/1‑91 (MQ=255)
                                                                   gATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAAt                         >  2:1512632/1‑91 (MQ=255)
                                                                   gATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAAt                         <  2:3765487/91‑1 (MQ=255)
                                                                    aTTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATc                        <  1:3024702/91‑1 (MQ=255)
                                                                     ttAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATcc                       <  2:1112000/91‑1 (MQ=255)
                                                                       aGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCAt                     >  2:562588/1‑91 (MQ=255)
                                                                       aGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCAt                     >  2:132289/1‑91 (MQ=255)
                                                                        gATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATa                    >  2:1586025/1‑91 (MQ=255)
                                                                             tttAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAgg               <  2:2417345/91‑1 (MQ=255)
                                                                              ttaattaaAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGga              >  1:1889637/1‑91 (MQ=255)
                                                                                aattaaAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGgaga            <  2:1062219/91‑1 (MQ=255)
                                                                                   taaAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATaa         <  1:3463062/91‑1 (MQ=255)
                                                                                   taaAACAATGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATaa         >  2:1492823/1‑91 (MQ=255)
                                                                                    aaaaCAATGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATaaa        >  2:1145439/1‑91 (MQ=255)
                                                                                     aaaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAAt       <  2:1747049/91‑1 (MQ=255)
                                                                                     aaaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAAt       <  2:841806/91‑1 (MQ=255)
                                                                                     aaaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAAt       <  2:842768/91‑1 (MQ=255)
                                                                                     aaaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAAt       <  2:2782894/91‑1 (MQ=255)
                                                                                     aaaCAATGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAAt       >  1:2128718/1‑91 (MQ=255)
                                                                                      aaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATc      >  1:3540011/1‑91 (MQ=255)
                                                                                      aaCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATc      <  2:1380237/91‑1 (MQ=255)
                                                                                       aCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCg     <  1:403565/91‑1 (MQ=255)
                                                                                       aCATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCg     <  2:1097632/91‑1 (MQ=255)
                                                                                          ttGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCGaaa  >  1:2516134/1‑91 (MQ=255)
                                                                                          ttGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCGaaa  >  2:2378849/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATATTAATATATGACTTGTGTTGGTAATGATAACCTTGGTTAGAGTGTATTGTCATGTGATAGTTAAGATTAGATCTTTTAATTAAAACATTGTTTTAAATATCGATGACAAGGTCTAATGTAGGACGTGTAGATATAAGGTTCGGAATTATATAAATCCATAAAAGGAGATAAATCGAAA  >  USA300TCH1516_ALE/85745‑85925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: