BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,123,544	T→A	6.2%	N66I (AAT→ATT)	USA300HOU_RS10825 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,123,544	0	T	A	6.2%	510.7 / 18.3	146	N66I (AAT→ATT)	USA300HOU_RS10825	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (76/61); new base A (5/4); total (81/65)
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.86e-01

Only 100 of 162 total aligned reads displayed.

ATGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTAG  >  USA300TCH1516_ALE/2123454‑2123634

aTGTTCACTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAt                                                                                            <  1:2076597/91‑1 (MQ=255)

 tGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAtt                                                                                           >  1:3376009/1‑91 (MQ=255)

 tGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAtt                                                                                           <  1:3885558/91‑1 (MQ=255)

 tGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAtt                                                                                           <  1:2460964/91‑1 (MQ=255)

  gTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAttt                                                                                          >  1:2422414/1‑91 (MQ=255)

  gTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAttt                                                                                          >  1:4099589/1‑91 (MQ=255)

  gTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAttt                                                                                          >  1:639287/1‑91 (MQ=255)

   ttCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAtttt                                                                                         <  1:270663/91‑1 (MQ=255)

   ttCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAtttt                                                                                         <  1:270252/91‑1 (MQ=255)

     cATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCg                                                                                       <  1:1591738/91‑1 (MQ=255)

     cATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAATTTCg                                                                                       <  1:807812/91‑1 (MQ=255)

      aTTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAATTTCGt                                                                                      <  1:608803/91‑1 (MQ=255)

       tttGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTa                                                                                     >  1:3130589/1‑91 (MQ=255)

       tttGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTa                                                                                     <  1:3780848/91‑1 (MQ=255)

       tttGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTa                                                                                     <  1:376612/91‑1 (MQ=255)

           ttATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCa                                                                                 >  1:3192014/1‑91 (MQ=255)

            tATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCAt                                                                                >  1:1216842/1‑91 (MQ=255)

            tATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCAt                                                                                <  1:878161/91‑1 (MQ=255)

                  tcctATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTa                                                                          >  1:528258/1‑91 (MQ=255)

                   cctCTTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaa                                                                         <  1:3735363/91‑1 (MQ=255)

                    ctATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaa                                                                        <  1:971754/91‑1 (MQ=255)

                      aTTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaaaa                                                                      >  1:2436462/1‑91 (MQ=255)

                      aTTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaaaa                                                                      >  1:966115/1‑91 (MQ=255)

                      aTTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTAATTACCTaaaaa                                                                      <  1:2558771/91‑1 (MQ=255)

                        tAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaaaaaa                                                                    <  1:1179161/91‑1 (MQ=255)

                         aaGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAat                                                                   >  1:27458/1‑91 (MQ=255)

                          aGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAata                                                                  >  1:2713821/1‑91 (MQ=255)

                           gATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAatat                                                                 <  1:3825357/91‑1 (MQ=255)

                              gtttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcat                                                              <  1:3965654/91‑1 (MQ=255)

                              gtttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcat                                                              <  1:410132/91‑1 (MQ=255)

                              gtttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcat                                                              <  1:2796674/91‑1 (MQ=255)

                               tttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcatc                                                             <  1:2132156/91‑1 (MQ=255)

                                ttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATAtcatca                                                            >  1:2986040/1‑91 (MQ=255)

                                 tgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATcatcat                                                           >  1:3302393/1‑91 (MQ=255)

                                 tgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATcatcat                                                           >  1:3302264/1‑91 (MQ=255)

                                    ttttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATAtt                                                        >  1:986028/1‑91 (MQ=255)

                                     tttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATAttt                                                       >  1:1218426/1‑91 (MQ=255)

                                      ttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTa                                                      <  1:1042660/91‑1 (MQ=255)

                                      ttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTa                                                      <  1:653985/91‑1 (MQ=255)

                                       tCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTaa                                                     <  1:2821449/91‑1 (MQ=255)

                                       tCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTaa                                                     <  1:3065710/91‑1 (MQ=255)

                                       tCTTTAAATACTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTaa                                                     >  1:2084762/1‑91 (MQ=255)

                                        cTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAAc                                                    >  1:1608673/1‑91 (MQ=255)

                                         tttAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACa                                                   >  1:3315514/1‑91 (MQ=255)

                                          ttAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACAt                                                  <  1:201080/91‑1 (MQ=255)

                                            aaaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTa                                                >  1:1716008/1‑91 (MQ=255)

                                            aaaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTa                                                >  1:1715337/1‑91 (MQ=255)

                                             aaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTaa                                               <  1:3249150/91‑1 (MQ=255)

                                             aaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTaa                                               <  1:2681132/91‑1 (MQ=255)

                                              aTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTaaa                                              <  1:484810/91‑1 (MQ=255)

                                              aTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATGACCTAAAAAAATATCATCATATTTAACATTaaa                                              >  1:3403830/1‑91 (MQ=255)

                                              aTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATCACCTAAAAAAATATCATCATATTTAACATTaaa                                              >  1:3405150/1‑91 (MQ=255)

                                                gCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAg                                            <  1:480183/91‑1 (MQ=255)

                                                 cTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGc                                           >  1:2464133/1‑91 (MQ=255)

                                                 cTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGc                                           >  1:3704860/1‑91 (MQ=255)

                                                   aaaaTAATTGATTGCGTTTTATCATGCGTAAATACGAAATTTTCGTATGCATTACCTAAAAAAATATCAGCACATTTAACATTAAAAGCAc                                         >  1:188030/1‑91 (MQ=255)

                                                    aaaTAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACt                                        <  1:1141703/91‑1 (MQ=255)

                                                    aaaTAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACt                                        >  1:786724/1‑91 (MQ=255)

                                                    aaaTAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACt                                        >  1:1254079/1‑91 (MQ=255)

                                                     aataatTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTc                                       <  1:568937/91‑1 (MQ=255)

                                                     aataatTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTc                                       >  1:1574694/1‑91 (MQ=255)

                                                      ataatTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCa                                      <  1:2787812/91‑1 (MQ=255)

                                                       taatTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCat                                     >  1:1996045/1‑91 (MQ=255)

                                                        aatTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCata                                    <  1:365089/91‑1 (MQ=255)

                                                         atTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCatat                                   <  1:59932/91‑1 (MQ=255)

                                                         atTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCatat                                   >  1:822822/1‑91 (MQ=255)

                                                         atTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCatat                                   >  1:3873862/1‑91 (MQ=255)

                                                          ttGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCAtata                                  >  1:1999636/1‑91 (MQ=255)

                                                           tGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATAc                                 >  1:4085284/1‑91 (MQ=255)

                                                           tGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATAc                                 >  1:4085389/1‑91 (MQ=255)

                                                           tGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATAc                                 <  1:2490634/91‑1 (MQ=255)

                                                            gATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACt                                <  1:369145/91‑1 (MQ=255)

                                                                tCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGa                            <  1:3612841/91‑1 (MQ=255)

                                                                 cTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGaa                           >  1:4127142/1‑91 (MQ=255)

                                                                  tttttATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGaaa                          >  1:1199655/1‑91 (MQ=255)

                                                                    tttATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGt                        >  1:105754/1‑91 (MQ=255)

                                                                     ttATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTa                       >  1:38083/1‑91 (MQ=255)

                                                                     ttATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTa                       >  1:126382/1‑91 (MQ=255)

                                                                       aTCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTaaa                     >  1:4054267/1‑91 (MQ=255)

                                                                       aTCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTaaa                     >  1:887492/1‑91 (MQ=255)

                                                                        tCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAAc                    >  1:1671539/1‑91 (MQ=255)

                                                                        tCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAAc                    >  1:703672/1‑91 (MQ=255)

                                                                         cATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACt                   <  1:2816584/91‑1 (MQ=255)

                                                                          aTTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTa                  >  1:2792662/1‑91 (MQ=255)

                                                                              gTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCtt              <  1:334689/91‑1 (MQ=255)

                                                                              gTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCtt              <  1:941634/91‑1 (MQ=255)

                                                                                aaaTACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTa            <  1:84104/91‑1 (MQ=255)

                                                                                aaaTACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTa            <  1:85275/91‑1 (MQ=255)

                                                                                   tACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATg         >  1:2747682/1‑91 (MQ=255)

                                                                                   tACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATg         >  1:3940805/1‑91 (MQ=255)

                                                                                   tACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATg         >  1:3243201/1‑91 (MQ=255)

                                                                                   tACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATg         >  1:3266805/1‑91 (MQ=255)

                                                                                    aCGAAAATTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGt        >  1:974706/1‑91 (MQ=255)

                                                                                     cGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGtt       >  1:3407513/1‑91 (MQ=255)

                                                                                     cGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGtt       <  1:1883909/91‑1 (MQ=255)

                                                                                     cGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGtt       >  1:2358104/1‑91 (MQ=255)

                                                                                      gAAAATTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGttt      >  1:4008826/1‑91 (MQ=255)

                                                                                      gAAAATTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGttt      >  1:2976242/1‑91 (MQ=255)

                                                                                         aTTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAACGTTTGTa   >  1:294140/1‑91 (MQ=255)

                                                                                          ttttCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTag  <  1:688744/91‑1 (MQ=255)

ATGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTAG  >  USA300TCH1516_ALE/2123454‑2123634

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