BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,513,975	C→T	5.7%	R54Q (CGA→CAA)	USA300HOU_RS12940 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,513,975	0	C	T	5.7%	518.6 / 14.4	157	R54Q (CGA→CAA)	USA300HOU_RS12940	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (67/81); new base T (5/4); total (72/85)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 157 total aligned reads displayed.

AACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCATAA  >  USA300TCH1516_ALE/2513886‑2514064

 aCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGa                                                                                         <  1:3900289/91‑1 (MQ=255)

 aCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGa                                                                                         <  1:2122613/91‑1 (MQ=255)

 aCCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGa                                                                                         >  1:2667010/1‑91 (MQ=255)

  ccAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGat                                                                                        <  1:365894/91‑1 (MQ=255)

   cAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGatt                                                                                       <  1:181013/91‑1 (MQ=255)

    aaCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGatta                                                                                      <  1:2325129/91‑1 (MQ=255)

     aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                     >  1:116924/1‑91 (MQ=255)

     aCATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattat                                                                                     <  1:3460877/91‑1 (MQ=255)

      cATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattatt                                                                                    <  1:203735/91‑1 (MQ=255)

      cATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattatt                                                                                    >  1:3662388/1‑91 (MQ=255)

      cATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGattatt                                                                                    <  1:3136214/91‑1 (MQ=255)

       aTTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttt                                                                                   <  1:3970543/91‑1 (MQ=255)

       aTTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttt                                                                                   <  1:4019219/91‑1 (MQ=255)

       aTTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttt                                                                                   <  1:522470/91‑1 (MQ=255)

        ttGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGATTAtttt                                                                                  <  1:2290490/91‑1 (MQ=255)

        ttGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAtttt                                                                                  <  1:1679537/91‑1 (MQ=255)

         tGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTAttttt                                                                                 >  1:3710728/1‑91 (MQ=255)

          ggCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTtta                                                                                <  1:4047895/91‑1 (MQ=255)

           gCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCGATTTAAATTGATTAttttaag                                                                               <  1:846525/91‑4 (MQ=38)

           gCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttat                                                                               <  1:3050634/91‑1 (MQ=255)

           gCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttat                                                                               <  1:2573241/91‑1 (MQ=255)

           gCTATTACCATACGCTTGCCTGTTTTCATAACTGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttat                                                                               <  1:3553772/91‑1 (MQ=255)

             tATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTttatta                                                                             <  1:3744518/91‑1 (MQ=255)

              aTTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTAc                                                                            >  1:2557051/1‑91 (MQ=255)

                  ccATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGa                                                                        <  1:3708245/91‑1 (MQ=255)

                     tACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAAc                                                                     <  1:1453206/91‑1 (MQ=255)

                      aCGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACa                                                                    <  1:2954245/91‑1 (MQ=255)

                       cGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACAt                                                                   <  1:1547164/91‑1 (MQ=255)

                       cGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACAt                                                                   >  1:655063/1‑91 (MQ=255)

                         cTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAg                                                                 <  1:333642/91‑1 (MQ=255)

                          tttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGt                                                                >  1:4115664/1‑91 (MQ=255)

                          tttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGt                                                                <  1:3412914/91‑1 (MQ=255)

                           ttCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAAGAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTg                                                               <  1:4092757/91‑1 (MQ=255)

                            tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                              >  1:3731420/1‑91 (MQ=255)

                            tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                              <  1:724732/91‑1 (MQ=255)

                            tCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGc                                                              <  1:2494838/91‑1 (MQ=255)

                             ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGtaga                                                             >  1:2367046/1‑88 (MQ=255)

                             ccTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGcc                                                             >  1:306526/1‑91 (MQ=255)

                               tGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCAt                                                           <  1:707462/91‑1 (MQ=255)

                                gTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCATAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCAtt                                                          >  1:3502471/1‑91 (MQ=255)

                                gTTTTCACAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCAtt                                                          <  1:2068297/91‑1 (MQ=255)

                                  tttCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCAtttt                                                        >  1:2742827/1‑91 (MQ=255)

                                  tttCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCAtttt                                                        >  1:288598/1‑91 (MQ=255)

                                       tAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTgg                                                   >  1:3443647/1‑91 (MQ=255)

                                         aCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgag                                                 >  1:2639445/1‑91 (MQ=255)

                                         aCAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAAGTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgag                                                 >  1:3601760/1‑91 (MQ=255)

                                          cAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGgaga                                                <  1:3238802/91‑1 (MQ=255)

                                            gTTCCAGCTAATCCACGACCTTTTCTTAAAATAATAAATTTAAATCGATTATTTTTATTACCTGAAAAATAGTGCCATTTTATTGGAgatg                                              >  1:1237850/1‑91 (MQ=255)

                                             ttCCAGCTAATCCGCGACCTTTTCTTAAAATAATCAATTTAAAGCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                             <  1:4144585/91‑1 (MQ=255)

                                             ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                             <  1:2481703/91‑1 (MQ=255)

                                             ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                             <  1:2895690/91‑1 (MQ=255)

                                             ttCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatga                                             <  1:1341737/91‑1 (MQ=255)

                                              tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                            >  1:2589863/1‑91 (MQ=255)

                                              tCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAgatgat                                            <  1:4044493/91‑1 (MQ=255)

                                               ccAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGatgatg                                           >  1:1591783/1‑91 (MQ=255)

                                                cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                          >  1:3069940/1‑91 (MQ=255)

                                                cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                          >  1:1935534/1‑91 (MQ=255)

                                                cAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATgg                                          >  1:3639342/1‑91 (MQ=255)

                                                 aGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGt                                         <  1:633458/91‑1 (MQ=255)

                                                 aGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGt                                         >  1:1104879/1‑91 (MQ=255)

                                                  gCTAATCCACGTCCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGtt                                        >  1:1556231/1‑91 (MQ=255)

                                                  gCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGtt                                        >  1:1153475/1‑91 (MQ=255)

                                                  gCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGtt                                        >  1:269790/1‑91 (MQ=255)

                                                  gCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGta                                        >  1:1506808/1‑90 (MQ=255)

                                                   cTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGttt                                       >  1:1035637/1‑91 (MQ=255)

                                                    tAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTg                                      <  1:422413/91‑1 (MQ=255)

                                                    tAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTg                                      <  1:3039529/91‑1 (MQ=255)

                                                    tAATCCACGACCTTTACTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAATCATAGTGCCATTTTATTGGAGATGATGGTTTg                                      >  1:3536537/1‑91 (MQ=255)

                                                      aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                                    >  1:1270273/1‑91 (MQ=255)

                                                      aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                                    <  1:3811348/91‑1 (MQ=255)

                                                      aTCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGtt                                    <  1:1863914/91‑1 (MQ=255)

                                                         cACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGa                                 >  1:3389343/1‑91 (MQ=255)

                                                           cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                               >  1:929638/1‑91 (MQ=255)

                                                           cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                               <  1:1404548/91‑1 (MQ=255)

                                                           cGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGAtt                               >  1:1797091/1‑91 (MQ=255)

                                                            gACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTc                              >  1:3949543/1‑91 (MQ=255)

                                                              ccTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCAt                            <  1:5132/91‑1 (MQ=255)

                                                                ttttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAg                          <  1:3869776/91‑1 (MQ=255)

                                                                 tttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGa                         <  1:486196/91‑1 (MQ=255)

                                                                  ttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaa                        >  1:3462060/1‑91 (MQ=255)

                                                                  ttCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaa                        <  1:4129607/91‑1 (MQ=255)

                                                                   tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                       <  1:202168/91‑1 (MQ=255)

                                                                   tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                       >  1:1741512/1‑91 (MQ=255)

                                                                   tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                       <  1:3628703/91‑1 (MQ=255)

                                                                   tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                       >  1:3072607/1‑91 (MQ=255)

                                                                   tCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGaaa                       <  1:2842014/91‑1 (MQ=255)

                                                                       aaaaTAATCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGa                   >  1:81587/1‑91 (MQ=255)

                                                                          ataatCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg                >  1:2831372/1‑91 (MQ=255)

                                                                          ataatCGATTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTg                >  1:1586518/1‑91 (MQ=255)

                                                                           taatCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGc               <  1:1725335/91‑1 (MQ=255)

                                                                              tCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGc            <  1:2940990/91‑1 (MQ=255)

                                                                               cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa           >  1:224128/1‑91 (MQ=255)

                                                                               cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa           >  1:1200524/1‑91 (MQ=255)

                                                                               cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa           >  1:1200446/1‑91 (MQ=255)

                                                                               cAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGAAGGTTTGTTAGATTCATAGAAAGCGATTGCCGCa           >  1:3397586/1‑91 (MQ=255)

                                                                                gaTTTAAATTGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaa          >  1:1407108/2‑91 (MQ=255)

                                                                                aaTTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCaa          >  1:1934851/1‑91 (MQ=255)

                                                                                      aaTCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCAt    <  1:3943652/91‑1 (MQ=255)

                                                                                      aaTCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCAt    >  1:2309690/1‑91 (MQ=255)

                                                                                       aTCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCATa   <  1:2910764/91‑1 (MQ=255)

AACCAACATTGGCTATTACCATACGCTTTCCTGTTTTCATAACAGTTCCAGCTAATCCACGACCTTTTCTTAAAATAATCAATTTAAATCGATTATTTTTATTACCTGAAACATAGTGCCATTTTATTGGAGATGATGGTTTGTTAGATTCATAGAAAGCGATTGCCGCAAAATCATAA  >  USA300TCH1516_ALE/2513886‑2514064

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