BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	682,085	Δ1 bp	5.9%	coding (74/504 nt)	USA300TCH1516_00616 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	682,085	0	C	.	5.9%	949.1 / 49.9	273	coding (74/504 nt)	USA300TCH1516_00616	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (138/117); new base . (13/3); total (152/121)
Fisher's exact test for biased strand distribution p-value = 3.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 287 total aligned reads displayed.

TAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACT  >  USA300TCH1516_ALE/681995‑682175

taAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAAc                                                                                            <  1:1181494/91‑1 (MQ=255)

taAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAAc                                                                                            >  1:3526841/1‑91 (MQ=255)

 aaaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAAcc                                                                                           >  1:1052714/1‑91 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCCCGACCTAAACTCTTAGAAATCTGCCTGAT‑GTTTAAACca                                                                                          <  1:3144227/91‑1 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACca                                                                                          >  1:3364979/1‑91 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACca                                                                                          >  1:2454436/1‑91 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACca                                                                                          >  1:3822936/1‑91 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACca                                                                                          <  2:2130672/91‑1 (MQ=255)

  aaCACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACca                                                                                          >  2:400580/1‑91 (MQ=255)

    cacCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACcatc                                                                                        <  1:851486/91‑1 (MQ=255)

    cacCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACcatc                                                                                        >  2:3564241/1‑91 (MQ=255)

     acCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACcatca                                                                                       >  2:3355303/1‑91 (MQ=255)

       cATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCAtt                                                                                     <  1:2039534/91‑1 (MQ=255)

         tCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTa                                                                                   <  1:3596078/91‑1 (MQ=255)

         ccAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGGTTTAAA‑CATCAttca                                                                                   <  2:302633/90‑3 (MQ=255)

            aCAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAAt                                                                                >  2:1497147/1‑91 (MQ=255)

            aCAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAAt                                                                                >  2:3273209/1‑91 (MQ=255)

            aCAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAAt                                                                                >  2:1939208/1‑91 (MQ=255)

              aGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATAc                                                                              <  2:658482/91‑1 (MQ=255)

               gAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACt                                                                             >  2:2284256/1‑91 (MQ=255)

                aaCGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACtt                                                                            >  2:100849/1‑91 (MQ=255)

                 aCGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACttt                                                                           >  1:2068065/1‑91 (MQ=255)

                 aCGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACttt                                                                           >  2:147442/1‑91 (MQ=255)

                  cGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTg                                                                          >  2:3895378/1‑91 (MQ=255)

                    cTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCg                                                                        <  2:650466/91‑1 (MQ=255)

                     ttATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGa                                                                       <  1:2238746/91‑1 (MQ=255)

                      tATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGAt                                                                      >  1:130441/1‑91 (MQ=255)

                       aTTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATa                                                                     >  1:1932353/1‑91 (MQ=255)

                       aTTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATa                                                                     >  2:1621075/1‑91 (MQ=255)

                        tttACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAg                                                                    >  1:1014256/1‑91 (MQ=255)

                         ttACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGc                                                                   >  1:3451080/1‑91 (MQ=255)

                          tACTGCTCTATAAAGAACTGCATCATCACGACCTAAGATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCt                                                                  <  2:1789545/91‑1 (MQ=255)

                             ttCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGa                                                               <  2:1354113/91‑1 (MQ=255)

                             ttCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGa                                                               <  2:3090455/91‑1 (MQ=255)

                              tctcTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGAc                                                              >  2:2297023/1‑91 (MQ=255)

                               ctctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACg                                                             <  1:647958/91‑1 (MQ=255)

                                tctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgt                                                            >  2:1676857/1‑91 (MQ=255)

                                tctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgt                                                            <  1:1576052/91‑1 (MQ=255)

                                tctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGATAGACgt                                                            >  1:1963374/1‑91 (MQ=255)

                                 ctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgtg                                                           >  2:2680956/1‑91 (MQ=255)

                                 ctATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgtg                                                           <  1:1475286/91‑1 (MQ=255)

                                  tataAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgtgt                                                          <  2:55420/91‑1 (MQ=255)

                                  tataAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACgtgt                                                          >  2:909916/1‑91 (MQ=255)

                                    taAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGttt                                                        >  1:2645677/1‑91 (MQ=255)

                                    taAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGttt                                                        >  2:1992188/1‑91 (MQ=255)

                                      aaGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTta                                                      <  2:2227050/91‑1 (MQ=255)

                                           cTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTc                                                 <  2:2179851/91‑1 (MQ=255)

                                                   tcaCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCtt                                         <  2:351190/91‑1 (MQ=255)

                                                   tcaCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCtt                                         <  2:931759/91‑1 (MQ=255)

                                                     aCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTa                                       >  2:2425262/1‑91 (MQ=255)

                                                     aCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTa                                       <  2:1746383/91‑1 (MQ=255)

                                                     aCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTa                                       <  2:1569191/91‑1 (MQ=255)

                                                     aCGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTa                                       <  1:703839/91‑1 (MQ=255)

                                                       gACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACg                                     >  1:2614727/1‑91 (MQ=255)

                                                       gACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACg                                     >  1:1868978/1‑91 (MQ=255)

                                                        aCCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGt                                    <  2:2893062/91‑1 (MQ=255)

                                                         ccTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtt                                   <  2:3833970/91‑1 (MQ=255)

                                                          cTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGttt                                  >  1:2340651/1‑91 (MQ=255)

                                                           tAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttt                                 <  1:446473/91‑1 (MQ=255)

                                                           tAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttt                                 <  1:3064111/91‑1 (MQ=255)

                                                            aaaaTCTTTGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGttttt                                >  2:1841773/1‑91 (MQ=255)

                                                            aaaaTCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGttttt                                <  1:2756581/91‑1 (MQ=255)

                                                              aaTCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGttttttt                              >  2:847295/1‑91 (MQ=255)

                                                               aTCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttc                             >  1:3804381/1‑91 (MQ=255)

                                                                tCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGttttttctt                            <  2:28882/91‑4 (MQ=255)

                                                                  ttAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttt                          <  1:941153/91‑1 (MQ=255)

                                                                  ttAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttt                          <  2:1107329/91‑1 (MQ=255)

                                                                  ttAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttt                          <  1:1389321/91‑1 (MQ=255)

                                                                   tAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGTTAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttctttt                         >  1:2935728/1‑91 (MQ=255)

                                                                   tAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttctttt                         <  2:1224897/91‑1 (MQ=255)

                                                                    aGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttt                        >  1:2017784/1‑91 (MQ=255)

                                                                    aGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttt                        >  2:2020384/1‑91 (MQ=255)

                                                                    aGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttt                        >  1:2509162/1‑91 (MQ=255)

                                                                     gAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttctttttt                       <  1:486438/91‑1 (MQ=255)

                                                                      aaaTCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttttt                      >  1:2974387/1‑91 (MQ=255)

                                                                      aaaTCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttttt                      >  2:2650498/1‑91 (MQ=255)

                                                                      aaaTCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGtttttttcttttttt                      >  2:858481/1‑91 (MQ=255)

                                                                        aTCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGTGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACG‑TTTTTTCTTTTTTTCAt                   >  1:329201/1‑91 (MQ=255)

                                                                          cTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATg                  <  2:2777182/91‑1 (MQ=255)

                                                                           taaaTGATGGTTTAAA‑CATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                 >  1:3961303/3‑91 (MQ=255)

                                                                           taaaTGAT‑GGTTTAAACATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                 >  1:1533927/3‑91 (MQ=255)

                                                                           tGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                 >  2:316356/1‑91 (MQ=255)

                                                                           tGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                 >  1:2181433/1‑91 (MQ=255)

                                                                            gAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGta                <  2:1057191/91‑1 (MQ=255)

                                                                              atgatgGTTTAAA‑CATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtata              >  1:2428132/1‑91 (MQ=255)

                                                                              atgat‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGGCGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtata              <  2:3469335/91‑1 (MQ=255)

                                                                              atgat‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtata              <  2:284816/91‑1 (MQ=255)

                                                                               tgat‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTatat             >  1:869589/1‑91 (MQ=255)

                                                                                 at‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcc           >  2:854808/1‑91 (MQ=255)

                                                                                 at‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcc           >  2:2163143/1‑91 (MQ=255)

                                                                                  t‑gTTTAAACCATCATTTATTTACTTTGCGATAGCTAACCGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTATTTCATGTATAtcct          <  2:670176/91‑1 (MQ=255)

                                                                                  t‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcct          >  1:2140388/1‑91 (MQ=255)

                                                                                  t‑gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcct          >  1:1836819/1‑91 (MQ=255)

                                                                                    gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcctc         <  2:2515323/91‑1 (MQ=255)

                                                                                    gTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcctc         <  2:3711671/91‑1 (MQ=255)

                                                                                     tttAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcctct        >  2:3942456/1‑90 (MQ=255)

                                                                                     tttAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcctcc        >  1:3648602/1‑91 (MQ=255)

                                                                                     tttAAA‑CATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCTTGTATATcctcct       >  1:1669083/1‑91 (MQ=255)

                                                                                      ttAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATcctcct       >  2:1292958/1‑91 (MQ=255)

                                                                                      ttAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATcctcct       >  2:1165199/1‑91 (MQ=255)

TAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGAT‑GTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACT  >  USA300TCH1516_ALE/681995‑682175

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