BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,430,695	C→G	6.7%	intergenic (+804/‑56)	ybiT → / → lysC	putative ABC transporter ATP‑binding protein YbiT/Aspartokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,430,695	0	C	G	6.7%	733.5 / 32.7	211	intergenic (+804/‑56)	ybiT/lysC	putative ABC transporter ATP‑binding protein YbiT/Aspartokinase
Reads supporting (aligned to +/- strand): ref base C (106/90); new base G (10/4); total (116/95)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 215 total aligned reads displayed.

TGGACTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTGGCGGATC  >  USA300TCH1516_ALE/1430605‑1430785

tgGACTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTc                                                                                             <  1:2549045/91‑1 (MQ=255)

  gACTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCg                                                                                           >  1:1663886/1‑91 (MQ=255)

   aCTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGa                                                                                          >  2:3221709/1‑91 (MQ=255)

   aCTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGa                                                                                          >  1:612432/1‑91 (MQ=255)

   aCTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGa                                                                                          <  2:3356655/91‑1 (MQ=255)

    ctgtttGGTTGAGGGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGGGCAAGTTCGCATTCCGaa                                                                                         <  2:1890211/88‑1 (MQ=255)

    cTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGaa                                                                                         >  2:3239907/1‑91 (MQ=255)

     tttttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAAc                                                                                        >  2:825117/1‑91 (MQ=255)

      ttttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACt                                                                                       <  1:971042/91‑1 (MQ=255)

      ttttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACt                                                                                       <  1:1701595/91‑1 (MQ=255)

      ttttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACt                                                                                       >  1:3478823/1‑91 (MQ=255)

       tttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACtt                                                                                      <  1:3862561/91‑1 (MQ=255)

       tttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACtt                                                                                      >  1:3570725/1‑91 (MQ=255)

       tttGGCTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAACAAGTTCGGAATGCGAACtt                                                                                      <  2:3577434/91‑1 (MQ=38)

        ttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTg                                                                                     >  2:3736620/1‑91 (MQ=255)

        ttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTg                                                                                     >  1:3770335/1‑91 (MQ=255)

        ttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTg                                                                                     <  2:1599425/91‑1 (MQ=255)

        ttGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGCGCATCACTTGTATATAGATTAAGAACAAGTTCGGAATGCGAACTTg                                                                                     <  1:3074903/91‑1 (MQ=37)

         tGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGt                                                                                    >  1:1796823/1‑91 (MQ=255)

         tGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGt                                                                                    <  2:3827231/91‑1 (MQ=255)

          ggTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGtt                                                                                   >  1:1204462/1‑91 (MQ=255)

           gTTAATAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑c                                                                                 <  2:3146288/91‑1 (MQ=255)

           gTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑c                                                                                 >  1:249799/1‑91 (MQ=255)

           gTTAAGAGCTGACAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑c                                                                                 >  1:2970734/1‑91 (MQ=255)

            ttAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑Ct                                                                                <  1:284893/91‑1 (MQ=255)

             tAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑Ctt                                                                               <  1:2051053/91‑1 (MQ=255)

              aaGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAACAAGTTCGGAATGCGAACTTGCT‑Ctta                                                                              <  1:3558704/91‑2 (MQ=25)

               agagCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑Ctttt                                                                             <  2:207544/91‑1 (MQ=255)

                gagCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑Cttttt                                                                            >  1:2354543/1‑91 (MQ=255)

                    tGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTttat                                                                        >  2:3544434/1‑91 (MQ=255)

                    tGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACCTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTttat                                                                        >  2:3186335/1‑91 (MQ=255)

                     gagaGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTttatt                                                                       >  2:876067/1‑91 (MQ=255)

                     gagaGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTttatt                                                                       <  1:1430517/91‑1 (MQ=255)

                         gTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATtg                                                                   >  2:2283594/1‑91 (MQ=255)

                          tttGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATtgt                                                                  <  1:3531799/91‑1 (MQ=255)

                           ttGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATtgtg                                                                 <  2:781332/91‑1 (MQ=255)

                            tGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATtgtgt                                                                >  1:2870066/1‑91 (MQ=255)

                              tCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGc                                                              <  2:2242452/91‑1 (MQ=255)

                              tCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGc                                                              <  1:1262675/91‑1 (MQ=255)

                              tCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGc                                                              <  1:3818012/91‑1 (MQ=255)

                              tCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGc                                                              <  1:1678158/91‑1 (MQ=255)

                              tCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGc                                                              <  2:1777127/91‑1 (MQ=255)

                               cATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGcc                                                             <  1:3247558/91‑1 (MQ=255)

                               cATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGcc                                                             >  1:3405589/1‑91 (MQ=255)

                                   attTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTc                                                         >  2:2882872/1‑91 (MQ=255)

                                    tttAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTCGTT‑CTTTTTTTATTATTGTGTGCCCTTcc                                                        <  2:3385043/91‑1 (MQ=255)

                                     ttAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCt                                                       <  1:792086/91‑1 (MQ=255)

                                      tttaaaTAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTa                                                      >  1:549819/4‑91 (MQ=255)

                                      tttaaaTAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTa                                                      >  1:1982113/4‑91 (MQ=255)

                                      tttaaaTAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTa                                                      >  1:2007735/4‑91 (MQ=255)

                                      tAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTa                                                      >  1:1831344/1‑91 (MQ=255)

                                        aaaaTAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAAc                                                    >  1:769290/1‑91 (MQ=255)

                                         aaaTAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACa                                                   >  1:2760242/1‑91 (MQ=255)

                                                ggAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTagga                                            >  2:2066568/1‑91 (MQ=255)

                                                ggAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTagga                                            >  1:1984670/1‑91 (MQ=255)

                                                  aGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTaggagg                                          <  1:1731650/91‑1 (MQ=255)

                                                   gTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAggagga                                         <  1:2212454/91‑1 (MQ=255)

                                                    tGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGAt                                        >  1:379851/1‑91 (MQ=255)

                                                          aCTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATg                                  >  1:3865714/1‑91 (MQ=255)

                                                           cTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATgg                                 >  1:1557649/1‑91 (MQ=255)

                                                            ttGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGt                                <  1:2283594/91‑1 (MQ=255)

                                                              gTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTTCCTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTa                               <  2:3650728/91‑1 (MQ=255)

                                                              gTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTaa                              >  1:810155/1‑91 (MQ=255)

                                                               tatataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAAc                             <  1:2009948/91‑1 (MQ=255)

                                                               tatataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAAc                             <  2:53910/91‑1 (MQ=255)

                                                                atataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACa                            <  2:2469542/91‑1 (MQ=255)

                                                                 tataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaa                           <  2:1932044/91‑1 (MQ=255)

                                                                 tataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaa                           <  1:3177101/91‑1 (MQ=255)

                                                                  ataGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaag                          <  1:401399/91‑1 (MQ=255)

                                                                   taGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaaga                         >  2:1403706/1‑91 (MQ=255)

                                                                   taGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaaga                         <  2:1391209/91‑1 (MQ=255)

                                                                   taGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACaaga                         <  1:1477391/91‑1 (MQ=255)

                                                                      aTTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgt                      <  1:2608927/91‑1 (MQ=255)

                                                                      aTTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgt                      <  2:2671560/91‑1 (MQ=255)

                                                                       ttAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtg                     <  1:1106087/91‑1 (MQ=255)

                                                                       ttAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtg                     >  2:1722312/1‑91 (MQ=255)

                                                                         aaGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtg                   >  1:1917194/1‑91 (MQ=255)

                                                                         aaGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtg                   <  2:2505213/91‑1 (MQ=255)

                                                                         aaGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtg                   <  2:1878342/91‑1 (MQ=255)

                                                                         aaGAACAAGTTCGGAATGCGAACTTGCT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGCGGCTTATATGGTAACAAGAAgtgtg                   >  1:1397617/1‑91 (MQ=25)

                                                                         aaGAACAAGTTCGGAATGCGAACTTGCT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtg                   >  1:580607/1‑91 (MQ=25)

                                                                         aaGAACAAGTTCGGAATGCGAACTTGCT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtg                   >  2:1434518/1‑91 (MQ=25)

                                                                          agagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtgt                  <  1:1436848/91‑1 (MQ=255)

                                                                          agagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtgt                  >  2:1145092/1‑91 (MQ=255)

                                                                          agagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAgtgtgt                  <  2:3755278/91‑1 (MQ=255)

                                                                           gagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGtt                 >  2:2691286/1‑91 (MQ=255)

                                                                           gagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGtt                 <  1:2185575/91‑1 (MQ=255)

                                                                           gagCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGtt                 <  2:1518745/91‑1 (MQ=255)

                                                                            agCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTg                >  2:1920228/1‑91 (MQ=255)

                                                                             gCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGa               <  2:1464580/91‑1 (MQ=255)

                                                                              cAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGaa              >  2:1173069/1‑91 (MQ=255)

                                                                               aaGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGaaa             <  2:2352437/91‑1 (MQ=255)

                                                                                   tCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTGGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTg         <  2:1033359/91‑1 (MQ=255)

                                                                                   tCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTg         >  2:718743/1‑91 (MQ=255)

                                                                                     gCATTCCGAACTTGTT‑CTTTTTTTCTTATTGTGTGCTCTCCCTAACAATCAGGTGGATTATATGGTAATAAGAAGTGTGTTGAACTTTGGc       >  2:1217820/1‑91 (MQ=255)

                                                                                     gCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTGGc       >  2:890731/1‑91 (MQ=255)

                                                                                      cATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTGGCg      >  1:2451283/1‑91 (MQ=255)

                                                                                         tCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATCAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTGGCGGat   >  2:3321750/1‑91 (MQ=255)

                                                                                          ccGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTg                <  1:2377243/77‑1 (MQ=255)

                                                                                          ccGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTg                >  2:2377243/1‑77 (MQ=255)

TGGACTTTTTGGTTAAGAGCTGAGAGTTTGTCATTATTTAAAAATAATGGAGTGCATCACTTGTATATAGATTAAGAGCAAGTTCGCATTCCGAACTTGTT‑CTTTTTTTATTATTGTGTGCCCTTCCTAACAATTAGGAGGATTATATGGTAACAAGAAGTGTGTTGAAATTTGGCGGATC  >  USA300TCH1516_ALE/1430605‑1430785

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