Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 396,540 G→A A208T (GCA→ACA)  efeM → putative iron uptake system component EfeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE396,5400GA98.0% 357.8 / ‑1.7 101A208T (GCA→ACA) efeMputative iron uptake system component EfeM
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (57/42);  total (58/43)
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.03e-01

Only 100 of 104 total aligned reads displayed.

AATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATA  >  USA300TCH1516_ALE/396453‑396624
                                                                                       |                                                                                    
aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCacac                                                                                   <  1:2926776/91‑1 (MQ=255)
aaTGAAGTTGCAACTTCAAAAAGCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCacac                                                                                   <  1:2619609/91‑1 (MQ=255)
  tGAAGATGCAACTTCTAAAATCACACGTGAAGAAGAAATGTATTCACATACAGGTTTATATGATTTTAAAGCGAACGTTGAAGGCACACaa                                                                                 <  1:3440554/91‑1 (MQ=255)
   gAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaa                                                                                >  1:1594248/1‑91 (MQ=255)
     aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaaaa                                                                              >  1:1232048/1‑91 (MQ=255)
     aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaaaa                                                                              >  1:464257/1‑91 (MQ=255)
      gTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaaaaa                                                                             >  1:2468458/1‑91 (MQ=255)
         gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAttt                                                                          >  1:2772659/1‑91 (MQ=255)
         gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAttt                                                                          <  1:1712112/91‑1 (MQ=255)
          cAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATCCAGATTTATATGATTTCAAAGCGCACGTTGAAGGCACACAAAAAATTTa                                                                         >  1:3672768/1‑91 (MQ=255)
          cAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTa                                                                         >  1:888824/1‑91 (MQ=255)
             cTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                      >  1:2545948/1‑91 (MQ=255)
              ttCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAc                                                                     >  1:289561/1‑91 (MQ=255)
               tCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACt                                                                    <  1:483649/91‑1 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACtt                                                                   <  1:3073262/91‑1 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACtt                                                                   >  1:3229359/1‑91 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACtt                                                                   >  1:1378065/1‑91 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACtt                                                                   <  1:970863/91‑1 (MQ=255)
                  aaaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAt                                                                 >  1:2196453/1‑91 (MQ=255)
                  aaaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAt                                                                 >  1:2196992/1‑91 (MQ=255)
                  aaaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAt                                                                 >  1:2586135/1‑91 (MQ=255)
                   aaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAtt                                                                >  1:1873828/1‑91 (MQ=255)
                    aaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAttt                                                               >  1:2565010/1‑91 (MQ=255)
                    aaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAttt                                                               >  1:727745/1‑91 (MQ=255)
                     aTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTa                                                              <  1:3006900/91‑1 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                             >  1:1232716/1‑91 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                             <  1:455933/91‑1 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                             <  1:456209/91‑1 (MQ=255)
                        acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAc                                                           >  1:2723701/1‑91 (MQ=255)
                        acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAc                                                           >  1:2723716/1‑91 (MQ=255)
                        acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAc                                                           >  1:3301738/1‑91 (MQ=255)
                         caGGTGAAGGAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAcc                                                          <  1:3550278/91‑1 (MQ=255)
                          aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCt                                                         >  1:229610/1‑91 (MQ=255)
                          aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCt                                                         >  1:226475/1‑91 (MQ=255)
                          aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCt                                                         >  1:226491/1‑91 (MQ=255)
                           ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTa                                                        >  1:2290447/1‑91 (MQ=255)
                           ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTa                                                        >  1:3676993/1‑91 (MQ=255)
                           ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTa                                                        >  1:3677008/1‑91 (MQ=255)
                            gTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTAt                                                       >  1:1300464/1‑91 (MQ=255)
                            gTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTAt                                                       >  1:1264259/1‑91 (MQ=255)
                              gaagaagaaCTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTAttt                                                     <  1:3543291/91‑1 (MQ=255)
                               aagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTAtttt                                                    >  1:491806/1‑91 (MQ=255)
                               aagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTActtt                                                    >  1:3334923/1‑91 (MQ=255)
                                 gaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTag                                                  >  1:195299/1‑91 (MQ=255)
                                 gaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTag                                                  >  1:195244/1‑91 (MQ=255)
                                  aagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTaga                                                 <  1:962208/91‑1 (MQ=255)
                                     aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                              >  1:1526853/1‑91 (MQ=255)
                                     aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                              >  1:1172508/1‑91 (MQ=255)
                                       aTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                            <  1:2620080/91‑1 (MQ=255)
                                        tttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaa                                           <  1:2158183/91‑1 (MQ=255)
                                         ttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaaa                                          <  1:1479835/91‑1 (MQ=255)
                                         ttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaaa                                          >  1:3254075/1‑91 (MQ=255)
                                         ttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaaa                                          >  1:3067456/1‑91 (MQ=255)
                                          tATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACCCAAAAAATTTATGACTTATTTAAACCTATCTTAGAGAAAAAAg                                         >  1:3706011/1‑91 (MQ=255)
                                          tATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAg                                         <  1:158229/91‑1 (MQ=255)
                                          tATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATCTTAGAGAAAAAAg                                         >  1:3705947/1‑91 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGGACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTAGTTTAGAGAAAAAAGa                                        <  1:3150388/91‑1 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                        >  1:2770630/1‑91 (MQ=255)
                                             tCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATa                                      <  1:2293481/91‑1 (MQ=255)
                                             tCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATa                                      <  1:1983284/91‑1 (MQ=255)
                                              cgcaTACAGATTTATATGCTTTTAAAGCGAACGTTGAAGGTGCACAAAAAATTTATGATTTATTCAAACCAATTCTAGAGAAAAAAGATaa                                     >  1:1794232/3‑91 (MQ=255)
                                              cacaTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaa                                     <  1:2585277/91‑1 (MQ=255)
                                                  tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaaa                                 <  1:2422849/91‑1 (MQ=255)
                                                   aCAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                <  1:40295/91‑1 (MQ=255)
                                                   aCAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                <  1:3418486/91‑1 (MQ=255)
                                                    cAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATGAAAAAtt                               >  1:1836370/1‑91 (MQ=255)
                                                     aGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                              >  1:1177968/1‑91 (MQ=255)
                                                     aGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                              <  1:1362592/91‑1 (MQ=255)
                                                     aGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                              >  1:2511455/1‑91 (MQ=255)
                                                        tttATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGt                           >  1:229428/1‑91 (MQ=255)
                                                           atatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgat                        <  1:48697/91‑1 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTGAGtgatg                       >  1:1833337/1‑91 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                       <  1:722875/91‑1 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                       >  1:2621941/1‑91 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                       >  1:2621781/1‑91 (MQ=255)
                                                             atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatga                      >  1:3691145/1‑91 (MQ=255)
                                                             atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatga                      >  1:3691165/1‑91 (MQ=255)
                                                               gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGata                    <  1:2316010/91‑1 (MQ=255)
                                                               gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGata                    <  1:2875707/91‑1 (MQ=255)
                                                               gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGata                    <  1:2875679/91‑1 (MQ=255)
                                                                 ttttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATc                  <  1:1138347/91‑1 (MQ=255)
                                                                 ttttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATc                  <  1:3095553/91‑1 (MQ=255)
                                                                 ttttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATc                  <  1:1138731/91‑1 (MQ=255)
                                                                   ttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCa                <  1:1946007/91‑1 (MQ=255)
                                                                   ttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCa                <  1:3580273/91‑1 (MQ=255)
                                                                   ttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCa                <  1:3580974/91‑1 (MQ=255)
                                                                    tAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa               <  1:113468/91‑1 (MQ=255)
                                                                    tAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa               <  1:3423471/91‑1 (MQ=255)
                                                                     aaaGCGAACGTTGAAGGCACACAAATAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa              <  1:1765094/91‑1 (MQ=255)
                                                                      aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAAt             <  1:235413/91‑1 (MQ=255)
                                                                      aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAAt             >  1:1131489/1‑91 (MQ=255)
                                                                        gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGa           >  1:3527082/1‑91 (MQ=255)
                                                                         cGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGaa          >  1:1301738/1‑91 (MQ=255)
                                                                           aaCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt        <  1:3552563/91‑1 (MQ=255)
                                                                           aaCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt        >  1:1340528/1‑91 (MQ=255)
                                                                           aaCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt        >  1:2733337/1‑91 (MQ=255)
                                                                           aaCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt        <  1:2536668/91‑1 (MQ=255)
                                                                            aCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACtt       <  1:216491/91‑1 (MQ=255)
                                                                            aCGTCGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACtt       <  1:1505973/91‑1 (MQ=255)
                                                                              gTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCg     >  1:3377328/1‑91 (MQ=255)
                                                                                       |                                                                                    
AATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATA  >  USA300TCH1516_ALE/396453‑396624

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