BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	109,970	T→C	7.6%	intergenic (+49/‑172)	plc → / → USA300HOU_RS00515	1‑phosphatidylinositol phosphodiesterase/putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	109,970	0	T	C	7.6%	463.3 / 19.2	144	intergenic (+49/‑172)	plc/USA300HOU_RS00515	1‑phosphatidylinositol phosphodiesterase/putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (88/45); new base C (6/5); total (94/50)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 157 total aligned reads displayed.

TATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGATAATA  >  USA300TCH1516_ALE/109880‑110060

tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAat                                                                                            >  1:2383496/1‑91 (MQ=255)

tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAat                                                                                            >  1:2382577/1‑91 (MQ=255)

  tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAatta                                                                                          >  1:330357/1‑91 (MQ=255)

  tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAatta                                                                                          >  1:329016/1‑91 (MQ=255)

   gatgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAattat                                                                                         <  1:1824146/91‑1 (MQ=255)

    atgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAattatt                                                                                        >  1:337989/1‑91 (MQ=255)

    atgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAattatt                                                                                        <  1:1658210/91‑1 (MQ=255)

     tgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTAttt                                                                                       >  1:1140706/1‑91 (MQ=255)

       atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTc                                                                                     >  1:1002141/1‑91 (MQ=255)

        tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTc                                                                                     <  1:1010776/91‑1 (MQ=25)

        tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTc                                                                                     <  1:3071598/91‑1 (MQ=25)

          tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCact                                                                                  >  1:2732255/1‑91 (MQ=255)

            gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCactac                                                                                >  1:367500/1‑91 (MQ=255)

            gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCactac                                                                                >  1:1533650/1‑91 (MQ=255)

             tAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCatact                                                                                <  1:3559161/91‑5 (MQ=21)

              aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTCAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTc                                                                              >  1:3291041/1‑91 (MQ=255)

              aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTCAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTc                                                                              >  1:3291055/1‑91 (MQ=255)

               aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                                             >  1:3835912/1‑91 (MQ=255)

                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                            <  1:3955203/91‑1 (MQ=255)

                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                            <  1:3955394/91‑1 (MQ=255)

                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                            <  1:1940488/91‑1 (MQ=255)

                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                            <  1:3958032/91‑1 (MQ=255)

                 tGAAATTATAGATAGTAATAAATAAGGATTCATTAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTGTTTTCACTACTCGtt                                                                           <  1:3538135/91‑1 (MQ=255)

                 tGAAATTATAGATAGTAATAAATAAGGATTCAGTAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTGTTTTCACTACTCGtt                                                                           <  1:3538768/91‑1 (MQ=255)

                 tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                                           >  1:3074368/1‑91 (MQ=255)

                  gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                                          >  1:1230103/1‑91 (MQ=255)

                  gAAATTATAGATAGGAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                                          <  1:840379/91‑1 (MQ=255)

                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                         >  1:2769747/1‑91 (MQ=255)

                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                         <  1:729417/91‑1 (MQ=255)

                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                         <  1:2990780/91‑1 (MQ=255)

                    aaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttttt                                                                        >  1:870984/1‑91 (MQ=255)

                      ttATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAt                                                                      >  1:2118573/1‑91 (MQ=255)

                        atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                                    >  1:1415452/1‑91 (MQ=255)

                        atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                                    >  1:2832641/1‑91 (MQ=255)

                        atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                                    >  1:1415356/1‑91 (MQ=255)

                          agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTg                                                                  >  1:4118879/1‑91 (MQ=255)

                           gatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGa                                                                 >  1:63267/1‑91 (MQ=255)

                             tagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaaa                                                               >  1:287191/1‑91 (MQ=255)

                             tagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaaa                                                               >  1:286666/1‑91 (MQ=255)

                                 aataaataAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                                           <  1:36723/91‑1 (MQ=255)

                                 aataaataAGGAGTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                                           >  1:3193895/1‑91 (MQ=255)

                                   taaataaGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGt                                                         >  1:605821/1‑91 (MQ=255)

                                       taaGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaata                                                     <  1:502889/91‑1 (MQ=255)

                                          ggATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTACTTATTTTCACTACTCGTTTTTATTTTGAACATAAGTAATAAtt                                                  >  1:4147620/1‑91 (MQ=255)

                                           gATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                                 >  1:1646860/1‑91 (MQ=255)

                                           gATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                                 >  1:1647580/1‑91 (MQ=255)

                                            aTTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTca                                                >  1:930397/1‑91 (MQ=255)

                                               cAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTcaaca                                             >  1:1862521/1‑91 (MQ=255)

                                                 ataatGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                           >  1:3660392/1‑91 (MQ=255)

                                                 ataatGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                           >  1:3456600/1‑91 (MQ=255)

                                                 ataatGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                           <  1:2470475/91‑1 (MQ=255)

                                                 ataatGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                           >  1:1540504/1‑91 (MQ=255)

                                                    atGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTACTTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                        >  1:3509405/1‑91 (MQ=255)

                                                    atGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                        >  1:1249323/1‑91 (MQ=255)

                                                     tGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                       <  1:2979958/91‑1 (MQ=255)

                                                      gATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtattat                                      <  1:60882/91‑1 (MQ=255)

                                                      gATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtattat                                      <  1:60688/91‑1 (MQ=255)

                                                         atTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                   >  1:2428699/1‑91 (MQ=255)

                                                         atTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                   >  1:2136340/1‑91 (MQ=255)

                                                          ttAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAt                                  >  1:747022/1‑91 (MQ=255)

                                                          ttAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAt                                  >  1:1598575/1‑91 (MQ=255)

                                                             aGACGAGTA‑TGAAAATAGTTAGGTTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                               <  1:80979/91‑1 (MQ=255)

                                                             aGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                               <  1:80964/91‑1 (MQ=255)

                                                              gACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                              >  1:1847722/1‑91 (MQ=255)

                                                              gACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                              >  1:2371560/1‑91 (MQ=255)

                                                              gACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                              <  1:3895057/91‑1 (MQ=255)

                                                              gACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                              >  1:2806287/1‑91 (MQ=255)

                                                              aaCGAGTAGTGAAAATAATTAGAATCTAACTATTTTCA‑TACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                              >  1:3818921/2‑91 (MQ=18)

                                                               aCGAGTAGTGAAAATAATTAGAATCTAACTATTTTCA‑TACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                             >  1:3741761/1‑91 (MQ=18)

                                                                cGAGTAGTGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                             >  1:1168100/1‑91 (MQ=25)

                                                                cGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACa                            <  1:2332434/91‑1 (MQ=255)

                                                                 gAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg                           <  1:3621462/91‑1 (MQ=255)

                                                                 gAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg                           <  1:3621439/91‑1 (MQ=255)

                                                                  aGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa                          >  1:2885651/1‑91 (MQ=255)

                                                                   gTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAt                         >  1:1653732/1‑91 (MQ=255)

                                                                   gTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAt                         >  1:675070/1‑91 (MQ=255)

                                                                     a‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                       <  1:3681186/91‑1 (MQ=255)

                                                                      gtGAAAATAATTAGAATCTAACTATTTTCA‑TACTCGTCTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgt                      >  1:1447344/2‑91 (MQ=21)

                                                                         aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt                    >  1:2738669/1‑91 (MQ=255)

                                                                           aataatTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttgt                  >  1:1990459/1‑91 (MQ=38)

                                                                           aaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttgt                  <  1:4097309/91‑1 (MQ=255)

                                                                            ataattAGAATCTAACTATTTTCACTACTCGCTTGTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAAAGATTGCTtgtg                 >  1:3879712/5‑91 (MQ=37)

                                                                            aTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTtgtg                 <  1:774969/91‑1 (MQ=255)

                                                                             tAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGa                >  1:426919/1‑91 (MQ=255)

                                                                              aGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGaa               >  1:2247784/1‑91 (MQ=255)

                                                                              aGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGaa               >  1:2512298/1‑91 (MQ=255)

                                                                              aGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGaa               >  1:3521721/1‑91 (MQ=255)

                                                                                 tAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtt            >  1:3082281/1‑91 (MQ=255)

                                                                                 tAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtt            >  1:2337738/1‑91 (MQ=255)

                                                                                   gATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtttt          >  1:1854456/1‑91 (MQ=255)

                                                                                    aTTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAttttt         <  1:584714/91‑1 (MQ=255)

                                                                                    aTTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAttttt         >  1:536/1‑91 (MQ=255)

                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat      <  1:2543698/91‑1 (MQ=255)

                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat      <  1:261454/91‑1 (MQ=255)

                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat      <  1:2630339/91‑1 (MQ=255)

                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat      <  1:2890749/91‑1 (MQ=255)

                                                                                         aaTTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGataa    >  1:192264/1‑91 (MQ=255)

                                                                                         aaTTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGataa    >  1:4094848/1‑91 (MQ=255)

                                                                                          attattTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGataat   <  1:2655237/91‑1 (MQ=255)

                                                                                           ttattTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGataata  <  1:2228717/91‑1 (MQ=255)

TATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTA‑TGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGATAATA  >  USA300TCH1516_ALE/109880‑110060

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