BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	576,908	Δ1 bp	7.7%	intergenic (+331/‑59)	gltX → / → cysE	Glutamate‑‑tRNA ligase/Serine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	576,908	0	A	.	7.7%	652.8 / 44.5	195	intergenic (+331/‑59)	gltX/cysE	Glutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand): ref base A (78/102); new base . (4/11); total (82/113)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 202 total aligned reads displayed.

CTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGA  >  USA300TCH1516_ALE/576818‑576998

 ttAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACaa                                                                                            <  2:311617/91‑1 (MQ=255)

  tAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAAt                                                                                           <  1:2795316/91‑1 (MQ=255)

  tAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAAt                                                                                           <  2:2373727/91‑1 (MQ=255)

   aGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATa                                                                                          >  1:3592943/1‑91 (MQ=255)

   aGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATa                                                                                          >  1:575712/1‑91 (MQ=255)

   aGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCGAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATa                                                                                          <  2:3190151/91‑1 (MQ=255)

   aGTTTCCTTGATAATTAGGAAGCGCCCGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATa                                                                                          <  2:2976699/91‑1 (MQ=255)

    gTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATAc                                                                                         <  1:3880788/91‑1 (MQ=255)

      ttCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACaa                                                                                       <  1:3722690/91‑1 (MQ=255)

          ttGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGt                                                                                  >  1:1179764/1‑91 (MQ=38)

           tGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGt                                                                                  <  1:893901/91‑1 (MQ=255)

            gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAAC‑ATACAAGTTGTA‑t                                                                               <  1:2001498/91‑1 (MQ=255)

             aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGTA‑tt                                                                              <  2:1232844/91‑2 (MQ=37)

             aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGTA‑tt                                                                              <  2:1233071/91‑2 (MQ=38)

             aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGTA‑tt                                                                              <  1:2124677/91‑2 (MQ=37)

              tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑Tg                                                                              <  2:738069/91‑1 (MQ=255)

                 tttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTa                                                                           <  2:2357549/91‑1 (MQ=255)

                  ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAg                                                                          >  1:308843/1‑91 (MQ=255)

                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTT‑TGTTAgt                                                                         >  1:1943570/1‑91 (MQ=255)

                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAgt                                                                         >  1:2423613/1‑91 (MQ=255)

                    ggAAGCGCCCGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAATTTGTA‑TGTTAgtg                                                                        <  2:713057/91‑1 (MQ=255)

                    ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAgtg                                                                        <  1:255445/91‑1 (MQ=255)

                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAgtggt                                                                      >  1:2895973/1‑91 (MQ=255)

                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGTATTGTTAgaggc                                                                      <  1:363549/91‑5 (MQ=21)

                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAAC‑ATACAAGTTGTATTGTTAgaggc                                                                      <  1:921732/91‑5 (MQ=21)

                       aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAgtggtg                                                                     >  2:1367473/1‑91 (MQ=255)

                       aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAgtggtg                                                                     >  1:3212963/1‑91 (MQ=255)

                        gcgcCAGCAATATTATTAATGTTATGTGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGc                                                                    <  2:3359810/91‑1 (MQ=255)

                          gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCtt                                                                  >  2:1997862/1‑91 (MQ=255)

                           ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCttt                                                                 <  1:2781113/91‑1 (MQ=255)

                           ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCttt                                                                 >  1:1965916/1‑91 (MQ=255)

                                aaTATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGa                                                            >  1:1658478/1‑91 (MQ=255)

                                 atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGat                                                           <  2:2160939/91‑1 (MQ=255)

                                 atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGat                                                           <  1:2335420/91‑1 (MQ=255)

                                 atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGat                                                           >  1:395772/1‑91 (MQ=255)

                                 atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGat                                                           <  1:2335261/91‑1 (MQ=255)

                                   attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGatat                                                         <  1:2817584/91‑1 (MQ=255)

                                     tattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATAttt                                                       >  2:1181416/1‑91 (MQ=255)

                                     tattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATAttt                                                       >  2:3039908/1‑91 (MQ=255)

                                      attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTa                                                      <  1:1438081/91‑1 (MQ=255)

                                      attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTa                                                      >  1:2072284/1‑91 (MQ=255)

                                       ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTaa                                                     <  2:1045998/91‑1 (MQ=255)

                                        taATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAAt                                                    >  2:116907/1‑91 (MQ=255)

                                          aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAAttt                                                  <  2:1204319/91‑1 (MQ=255)

                                           tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                 >  1:70002/1‑91 (MQ=255)

                                           tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                 <  1:1997418/91‑1 (MQ=255)

                                           tGTCATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                 <  1:3494685/91‑1 (MQ=255)

                                            gTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCg                                                <  2:3181315/91‑1 (MQ=255)

                                                ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                            <  1:2264044/91‑1 (MQ=255)

                                                ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                            <  1:152220/91‑1 (MQ=255)

                                                ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                            <  1:3777481/91‑1 (MQ=255)

                                                ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                            >  2:392277/1‑91 (MQ=255)

                                                 tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                           >  1:3815753/1‑91 (MQ=255)

                                                  cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTaa                                                     >  2:1762714/1‑80 (MQ=255)

                                                   gatAAATTCAGAGTGGAACCGTGCGTAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTTGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                         <  2:264786/89‑1 (MQ=255)

                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                         >  2:3816800/1‑91 (MQ=255)

                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                         <  1:1543177/91‑1 (MQ=255)

                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATACTTAATTTCGCAGGTAc                                         >  2:3714402/1‑91 (MQ=255)

                                                     tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACtt                                       >  2:2071032/1‑91 (MQ=255)

                                                      aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                      >  1:1298872/1‑91 (MQ=255)

                                                      aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                      <  2:1725775/91‑1 (MQ=255)

                                                       aaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCa                                     <  1:3562684/91‑1 (MQ=255)

                                                         ttCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAt                                   <  1:2880079/91‑1 (MQ=255)

                                                          tCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAAtt                                  <  1:2179008/91‑1 (MQ=255)

                                                           cAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTa                                 >  2:3282456/1‑91 (MQ=255)

                                                            agagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTAt                                >  2:1834737/1‑91 (MQ=255)

                                                             gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                               >  1:3689742/1‑91 (MQ=255)

                                                              agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                              >  1:1052749/1‑91 (MQ=255)

                                                                tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                            <  2:3338353/91‑1 (MQ=255)

                                                                tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                            >  1:1139661/1‑91 (MQ=255)

                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  1:3180700/91‑1 (MQ=255)

                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  1:2692263/91‑1 (MQ=255)

                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  2:3507074/91‑1 (MQ=255)

                                                                  gAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAttt                          <  1:1310189/91‑1 (MQ=255)

                                                                  gAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAttt                          >  2:2531945/1‑91 (MQ=255)

                                                                   aaCCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTa                         >  1:2031218/1‑91 (MQ=255)

                                                                    aCCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTAc                        <  2:440613/91‑1 (MQ=255)

                                                                     ccGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACg                       >  2:639147/1‑91 (MQ=255)

                                                                      cGTGCGGAAGCGCCTCTAACAATACAGCTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                      <  2:1960403/91‑1 (MQ=255)

                                                                      cGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                      <  1:999354/91‑1 (MQ=255)

                                                                        tGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                    <  2:1522816/91‑1 (MQ=255)

                                                                         gCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                   <  2:2794/91‑1 (MQ=255)

                                                                          cGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGtt                  >  2:1082101/1‑91 (MQ=255)

                                                                           ggAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                 >  2:2265608/1‑91 (MQ=255)

                                                                           ggAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                 <  1:1347896/91‑1 (MQ=255)

                                                                            gAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                >  2:28137/1‑91 (MQ=255)

                                                                            gAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                >  2:1535597/1‑91 (MQ=255)

                                                                             aaGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCaa               <  1:3426629/91‑1 (MQ=255)

                                                                              aGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt              <  1:1729521/91‑1 (MQ=255)

                                                                              aGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt              <  2:3584791/91‑1 (MQ=255)

                                                                                    actAAC‑ATACAAGTTGTATTGTTAGAGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTga        >  1:224525/2‑91 (MQ=37)

                                                                                     ctAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag       <  2:1711465/91‑1 (MQ=255)

                                                                                      tAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg      >  1:1847370/1‑91 (MQ=255)

                                                                                      tAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg      >  1:946374/1‑91 (MQ=255)

                                                                                        aCAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag    <  2:1688427/91‑1 (MQ=255)

                                                                                        aCAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag    <  1:850182/91‑1 (MQ=255)

                                                                                        aCAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag    <  1:2373459/91‑1 (MQ=255)

                                                                                         cAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGaggagg   <  1:946987/91‑1 (MQ=255)

                                                                                          aaTACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga  <  1:2451395/91‑1 (MQ=255)

                                                                                          aaTACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga  <  1:1022618/91‑1 (MQ=255)

CTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTA‑TGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGA  >  USA300TCH1516_ALE/576818‑576998

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