BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,123,544	T→A	6.9%	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.9%	875.8 / 39.8	262	N66I (AAT→ATT)	USA300HOU_RS10825	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (118/126); new base A (8/10); total (126/136)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 284 total aligned reads displayed.

ATGTTCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTAG  >  USA300TCH1516_ALE/2123454‑2123634

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

   ttCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAattt                                                                                         <  1:3829387/91‑1 (MQ=255)

    tCATTTGTTGTTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTc                                                                                        <  1:3719978/91‑1 (MQ=255)

    tCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTc                                                                                        >  2:169411/1‑91 (MQ=255)

    tCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTc                                                                                        <  2:1148883/91‑1 (MQ=255)

    tCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTc                                                                                        <  1:658195/91‑1 (MQ=255)

    tCATTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTc                                                                                        >  1:2364027/1‑91 (MQ=255)

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

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

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

      aTTTGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGt                                                                                      <  2:2338574/91‑1 (MQ=255)

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

       tttGTTAGTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAATTTCGTa                                                                                     <  1:1947906/91‑1 (MQ=255)

        ttGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTAt                                                                                    >  1:565305/1‑91 (MQ=255)

        ttGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAAATTTCGTAt                                                                                    <  1:2070650/91‑1 (MQ=255)

         tGTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTAtt                                                                                   >  2:934832/1‑91 (MQ=255)

          gTTATTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTATTATCATTCGTAAACACGAAATTTTCGTATTc                                                                                  >  1:3319434/1‑91 (MQ=255)

             aTTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCAtt                                                                               <  2:442712/91‑1 (MQ=255)

             aTTCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCAtt                                                                               <  1:2463517/91‑1 (MQ=255)

              ttCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATTATTCGTAAATACGAAATTTTCGTATATATTa                                                                              >  2:3568386/1‑91 (MQ=255)

              ttCCTCCTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTa                                                                              <  2:3619850/91‑1 (MQ=255)

              ttCCTCCTATTAAGACGTATGTTTTTCCTTAAATGCTAAAACAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTa                                                                              >  2:620161/1‑91 (MQ=255)

               tcctccTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTAc                                                                             >  2:1556786/1‑91 (MQ=255)

               tcctccTATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTAc                                                                             >  2:603307/1‑91 (MQ=255)

                   cctATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaa                                                                         >  2:1662165/1‑91 (MQ=255)

                    ctATTAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaa                                                                        >  1:1138277/1‑91 (MQ=255)

                       ttAAGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTaaaaaa                                                                     >  1:2704267/1‑91 (MQ=255)

                          aGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAata                                                                  <  2:693968/91‑1 (MQ=255)

                          aGATGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAata                                                                  <  1:623454/91‑1 (MQ=255)

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

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

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

                            aTGTTTGTTTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatc                                                                <  2:1919292/91‑1 (MQ=255)

                             tgtttgttTTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatca                                                               >  1:1859743/1‑91 (MQ=255)

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

                               tttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcatc                                                             >  1:1309169/1‑91 (MQ=255)

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

                               tttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcatc                                                             >  2:401712/1‑91 (MQ=255)

                               tttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATatcatc                                                             >  1:3337952/1‑91 (MQ=255)

                                ttgtttTTCTTTAGATGCTAAAATAATTGATTTCTTTTTATCAGTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATAtcatca                                                            <  1:3878289/91‑1 (MQ=255)

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

                                ttgtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATAtcatca                                                            <  1:1231647/91‑1 (MQ=255)

                                  gtttTTCTTTAAATGCTAAGATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCata                                                          <  2:1969199/91‑1 (MQ=255)

                                  gtttTTCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCata                                                          >  2:2613997/1‑91 (MQ=255)

                                  gtatttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCata                                                          <  2:2553353/88‑1 (MQ=255)

                                   tttttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCatat                                                         >  2:2127406/1‑91 (MQ=255)

                                    ttttCTTTAAATGCTAAAATTATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATAtt                                                        >  2:1330728/1‑91 (MQ=255)

                                    ttttCTTTAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATAtt                                                        <  2:187059/91‑1 (MQ=255)

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

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

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

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

                                          ttAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACTAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACAt                                                  <  1:3071902/91‑1 (MQ=255)

                                          ttAAATGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACAt                                                  >  2:1738761/1‑91 (MQ=255)

                                            aaaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTa                                                <  2:3019223/91‑1 (MQ=255)

                                             aaTGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTaa                                               >  1:786708/1‑91 (MQ=255)

                                               tGCTAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTaaaa                                             >  2:1840104/1‑91 (MQ=255)

                                                  tAAAATAATTGATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCa                                          >  1:1541356/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                                                            gATTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACt                                >  1:1481693/1‑91 (MQ=255)

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

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

                                                             aTTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACtt                               >  1:3644876/1‑91 (MQ=255)

                                                             aTTTCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACtt                               <  2:3510132/91‑1 (MQ=255)

                                                              tttCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATAccta                              <  2:846034/91‑4 (MQ=255)

                                                              tttCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTa                              <  2:2445892/91‑1 (MQ=255)

                                                                tCTTTTTATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGa                            >  1:2180976/1‑91 (MQ=255)

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

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

                                                                   ttttATCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAg                         <  2:2148919/91‑1 (MQ=255)

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

                                                                       aTCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTaaa                     <  1:3680407/91‑1 (MQ=255)

                                                                       aTCATTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTaaa                     <  1:249891/91‑1 (MQ=255)

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

                                                                          aTTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTa                  <  1:146296/91‑1 (MQ=255)

                                                                          aTTCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTa                  <  2:2215314/91‑1 (MQ=255)

                                                                           ttCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTAt                 >  1:2812912/1‑91 (MQ=255)

                                                                           ttCGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTAt                 <  1:167562/91‑1 (MQ=255)

                                                                             cGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCt               >  1:3188518/1‑91 (MQ=255)

                                                                             cGTAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCt               >  2:1593929/1‑91 (MQ=255)

                                                                               tAAATACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCttt             <  2:1812180/91‑1 (MQ=255)

                                                                                aaaTACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTa            >  2:1073744/1‑91 (MQ=255)

                                                                                 aaTACGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTaa           <  1:492359/91‑1 (MQ=255)

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

                                                                                    aCGAAATTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGt        >  2:3725027/1‑91 (MQ=255)

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

                                                                                     cGAAAATTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATCTACTTAGAAAGTAAACTATCTTTAATGtt       >  1:138688/1‑91 (MQ=255)

                                                                                     cGAAAATTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGtt       >  1:3387279/1‑91 (MQ=255)

                                                                                       aaaTTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTg     <  1:1729966/91‑1 (MQ=255)

                                                                                        aaTTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCtt              >  2:3408244/1‑81 (MQ=255)

                                                                                         aTTTTCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTa   >  1:2735655/1‑91 (MQ=255)

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

                                                                                          ttttCGTATTCATTACCTAAAAAAATATCATCATATTTAACATTAAAAGCACTCATATACTTAGAAAGTAAACTATCTTTAATGTTTGTag  <  2:3351569/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: `‑`