BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,538,744	A→T	11.6%	intergenic (‑258/+69)	gpmA_2 ← / ← fieF	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase/Ferrous‑iron efflux pump FieF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,538,744	0	A	T	11.6%	805.3 / 83.1	272	intergenic (‑258/+69)	gpmA_2/fieF	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase/Ferrous‑iron efflux pump FieF
Reads supporting (aligned to +/- strand): ref base A (90/150); new base T (2/30); total (92/180)
Fisher's exact test for biased strand distribution p-value = 2.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

Only 100 of 281 total aligned reads displayed.

TCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCTGGG  >  USA300TCH1516_ALE/2538654‑2538834

 cATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTaaa                                                                                           <  1:2446030/91‑1 (MQ=255)

   tAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAAcc                                                                                         <  1:2651377/91‑1 (MQ=255)

    aaaTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCt                                                                                        <  1:3698514/91‑1 (MQ=255)

     taTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTa                                                                                       >  1:2903789/2‑91 (MQ=255)

     aaTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTa                                                                                       <  2:2596203/91‑1 (MQ=255)

      aTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaa                                                                                      <  1:2646634/91‑1 (MQ=255)

       tGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaa                                                                                     <  1:3496516/91‑1 (MQ=255)

       tGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaa                                                                                     <  2:1392699/91‑1 (MQ=255)

       tGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaaa                                                                                     >  2:2752697/1‑91 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaaa                                                                                    <  1:2252254/91‑1 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaaa                                                                                    <  1:2252037/91‑1 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaaa                                                                                    <  1:2026730/91‑1 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaaaa                                                                                    >  2:2416224/1‑91 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaaaa                                                                                    <  1:3639536/91‑1 (MQ=255)

           tCAATAGCCACCTAATAACAAATCATATTATACAACTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATa                                                                                 <  2:829031/91‑1 (MQ=255)

           tCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTAAAAATa                                                                                 <  1:3913646/91‑1 (MQ=255)

            cAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTAAAAATat                                                                                <  1:2748674/91‑2 (MQ=255)

             aaTGACAACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTACAAATAGa                                                                               >  2:3500235/1‑91 (MQ=255)

             aaTAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGa                                                                               <  1:1796121/91‑1 (MQ=255)

               tAACCACCTAATAACAAATCATATTATACACATTTGTTCGCTATTTTTATAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATg                                                                             <  2:3258892/91‑1 (MQ=255)

                aaCCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCAAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGt                                                                            <  2:1287687/91‑1 (MQ=255)

                  ccaccTAATAACAAATCATATTATACTCCTTTGTCCGCTATATTTCTAAGGTCTAAAAAATATTTTTAGGTAAACCAAAAAATAGATGtaa                                                                          <  2:332952/91‑1 (MQ=255)

                   caccTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAATATATTTTTAGGTAAACCTAAAAATAGATGtaat                                                                         >  2:1534675/1‑91 (MQ=255)

                        aataaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATaaaaa                                                                    <  1:1392777/91‑1 (MQ=255)

                           aaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGc                                                                 <  1:1164362/91‑1 (MQ=255)

                             cAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcct                                                               >  1:152448/1‑91 (MQ=255)

                              aaaTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctc                                                              >  1:2062776/1‑91 (MQ=255)

                              aaaTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctc                                                              <  1:3176762/91‑1 (MQ=255)

                               aaTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctcc                                                             >  2:1813309/1‑91 (MQ=255)

                               aaTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctcc                                                             <  2:3150935/91‑1 (MQ=255)

                                aTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctcct                                                            <  1:2840022/91‑1 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCctcctc                                                           >  2:2162463/1‑91 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCctcctc                                                           >  1:1929336/1‑91 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGCTAAACCTAAAAATAGATGTAATAAAAACGCttcctc                                                           >  2:3906695/1‑91 (MQ=255)

                                  cATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCa                                                          <  2:2344461/91‑1 (MQ=255)

                                  cATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCa                                                          <  1:959907/91‑1 (MQ=255)

                                    tattatACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGa                                                        >  2:3401012/1‑91 (MQ=255)

                                       tataCACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGatat                                                     >  1:3858310/1‑91 (MQ=255)

                                        ataCACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATAtt                                                    >  1:1061444/1‑91 (MQ=255)

                                           cacCTTAGTGCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATATAATAAAAACGCCTCCTCAGATATTtat                                                 <  2:489283/91‑1 (MQ=255)

                                            acCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTtata                                                <  1:252625/91‑1 (MQ=255)

                                             ccTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCGAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTtatat                                               <  1:254173/91‑1 (MQ=255)

                                              cgttGTTCGCTATGTTTCTAAGGGTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTtatata                                              <  2:1648679/89‑1 (MQ=255)

                                               tttGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTatatat                                             <  1:1077109/91‑1 (MQ=255)

                                                ttGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATc                                            <  1:2077041/91‑1 (MQ=255)

                                                 tGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCt                                           >  1:3414792/1‑91 (MQ=255)

                                                  gTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTa                                          <  2:2452093/91‑1 (MQ=255)

                                                  gTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTa                                          <  1:1607458/91‑1 (MQ=255)

                                                  gTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTa                                          >  2:2957821/1‑91 (MQ=255)

                                                     cGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGa                                       >  1:2155498/1‑91 (MQ=255)

                                                      gCTATATTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGaa                                      <  2:1763074/91‑1 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                     >  1:2596301/1‑91 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                     <  1:570280/91‑1 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                     <  1:2294609/91‑1 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATATATATATCTATGAAg                                     >  2:1976370/1‑91 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGACATTCACACATCTATGAAg                                     >  2:3336490/1‑91 (MQ=255)

                                                        tATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGa                                    >  1:2687873/1‑91 (MQ=255)

                                                           ttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGt                                 <  1:2201122/91‑1 (MQ=255)

                                                           ttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGt                                 <  1:3118716/91‑1 (MQ=255)

                                                            tttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGtt                                >  1:1386473/1‑91 (MQ=255)

                                                               cTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaa                             >  1:2699491/1‑91 (MQ=255)

                                                               cTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaa                             >  2:2679720/1‑91 (MQ=255)

                                                                 aaGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAAt                           <  2:3761587/91‑1 (MQ=255)

                                                                  aGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATa                          <  1:3455756/91‑1 (MQ=255)

                                                                   ggTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATAc                         <  1:140441/91‑1 (MQ=255)

                                                                    gTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACa                        <  2:769806/91‑1 (MQ=255)

                                                                    gTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACa                        <  1:544047/91‑1 (MQ=255)

                                                                     tttAGAAAATATTTTTAGGTAAGCGTAAGAATAGATGTAAGAAAAACGCCACCTCAGATATTTATATATCTATGAAGACGTTTAAATACAt                       <  2:3664001/91‑1 (MQ=255)

                                                                     tttAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACAt                       <  1:2485802/91‑1 (MQ=255)

                                                                     tttAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATAGTTATATATCTATGAAGACGTTTAAATACAt                       >  2:1297893/1‑91 (MQ=255)

                                                                       tAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATta                     <  2:3211139/91‑1 (MQ=255)

                                                                       tAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATta                     <  2:2463888/91‑1 (MQ=255)

                                                                        aaaaaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAGAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATtat                    <  2:1134922/91‑1 (MQ=255)

                                                                        aaaaaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATtat                    >  1:1716002/1‑91 (MQ=255)

                                                                           aaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGa                 <  1:454635/91‑1 (MQ=255)

                                                                           aaaTATTTTTAGGTAAACCTAAAAATAGAAGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGa                 <  2:358610/91‑1 (MQ=255)

                                                                             atatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATg               >  2:597216/1‑91 (MQ=255)

                                                                             atatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATg               >  2:2475761/1‑91 (MQ=255)

                                                                              tatTTTTAGGTTTACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATgg              >  1:1228798/1‑91 (MQ=255)

                                                                              tatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATgg              >  2:943463/1‑91 (MQ=255)

                                                                              tatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATgg              >  1:163990/1‑91 (MQ=255)

                                                                              tatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATgg              <  1:3224094/91‑1 (MQ=255)

                                                                               atTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGt             <  2:1054267/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGCTATTCATATATCTATGAAGACGTTTAAATACATTATAGATGGTc            >  2:1218025/1‑91 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGTAGACGTTTAAATACATTATAGATGGTc            <  1:145518/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTc            <  2:2047148/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTc            <  1:1597638/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCGAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTGATATATCTATGAAGACGTTTAAATACATGATAGATGGTc            <  2:724657/91‑1 (MQ=255)

                                                                                 ttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCt           <  1:3843565/91‑1 (MQ=255)

                                                                                 ttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCt           <  1:2284431/91‑1 (MQ=255)

                                                                                  tttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTg          <  2:8987/91‑1 (MQ=255)

                                                                                  tttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTg          <  1:1354332/91‑1 (MQ=255)

                                                                                     aGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGtt       <  1:754737/91‑1 (MQ=255)

                                                                                     aGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGtt       <  2:162647/91‑1 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc      >  1:461296/1‑91 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc      <  2:2017722/91‑1 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc      >  1:3789020/1‑91 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCGCAGATATTTATAGATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc      <  2:3130924/91‑1 (MQ=255)

                                                                                          aaCCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCTggg  <  2:963792/91‑1 (MQ=255)

                                                                                          aaCCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCTggg  >  1:2360205/1‑91 (MQ=255)

TCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCTGGG  >  USA300TCH1516_ALE/2538654‑2538834

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: `‑`