BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,440,007	Δ1 bp	coding (34/201 nt)	cspA_2 ←	Cold shock protein CspA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,440,007	0	C	.	96.0%	939.2 / 34.0	226	coding (34/201 nt)	cspA_2	Cold shock protein CspA
Reads supporting (aligned to +/- strand): ref base C (3/1); major base . (103/110); minor base T (3/6); total (109/117)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACA  >  USA300TCH1516_ALE/1439921‑1440094

tGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGt                                                                                    >  1:1308917/1‑91 (MQ=255)

 gATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGtt                                                                                   >  1:2490329/1‑91 (MQ=255)

 gATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGtt                                                                                   <  1:1882525/91‑1 (MQ=255)

   tttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaa                                                                                 <  2:862962/91‑1 (MQ=255)

   tttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaa                                                                                 <  1:1272837/91‑1 (MQ=255)

     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                               >  1:371661/1‑91 (MQ=255)

     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                               <  1:2398380/91‑1 (MQ=255)

     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                               <  2:1881498/91‑1 (MQ=255)

     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                               >  2:2626214/1‑91 (MQ=255)

      gTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAcc                                                                              <  1:632366/91‑1 (MQ=255)

       tAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCa                                                                             >  2:1309760/1‑91 (MQ=255)

        aaCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAt                                                                            >  1:1632766/1‑91 (MQ=255)

         aCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAtt                                                                           <  2:425354/91‑1 (MQ=255)

           cATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTa                                                                         >  1:1478473/1‑91 (MQ=255)

            aagatGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTaa                                                                        >  1:1391278/5‑91 (MQ=255)

            aTCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTaa                                                                        <  2:3205237/91‑1 (MQ=255)

             tCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAAc                                                                       >  1:2345859/1‑91 (MQ=255)

               ttGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTg                                                                     <  1:2607758/91‑1 (MQ=255)

               ttGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTg                                                                     <  1:383809/91‑1 (MQ=255)

               ttGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTg                                                                     <  2:1569103/91‑1 (MQ=255)

                tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGt                                                                    >  1:3505076/1‑91 (MQ=255)

                tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGt                                                                    >  2:1941567/1‑91 (MQ=255)

                  gTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTATT‑AGCGTTAAACCATTTAACTGTAc                                                                  >  2:2953575/1‑91 (MQ=255)

                      aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTg                                                              <  1:3390572/91‑1 (MQ=255)

                      aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTg                                                              <  1:1870236/91‑1 (MQ=255)

                      aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTg                                                              <  1:2794113/91‑1 (MQ=255)

                        tGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGtt                                                            <  2:1002795/91‑1 (MQ=255)

                         gCTGAAAAATGTACGAATACGTCATTTTCTTCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                           >  2:1043371/1‑91 (MQ=39)

                         gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                           >  2:420085/1‑91 (MQ=255)

                         gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                           <  2:2908216/91‑1 (MQ=255)

                          cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                          <  1:110022/91‑1 (MQ=255)

                          cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                          <  1:577963/91‑1 (MQ=255)

                          cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                          <  1:111355/91‑1 (MQ=255)

                          cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                          >  2:1245624/1‑91 (MQ=255)

                          cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCCTTTTCAGCGTTAAACCATTTAACTGTACCTTGttt                                                           >  1:685764/1‑91 (MQ=255)

                            gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCAt                                                        <  1:3015904/91‑1 (MQ=255)

                             aaaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATa                                                       >  1:608270/1‑91 (MQ=255)

                               aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAAt                                                     <  2:1046897/91‑1 (MQ=255)

                                aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                    >  1:1453349/1‑91 (MQ=255)

                                aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                    <  1:969399/91‑1 (MQ=255)

                                aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                    <  2:3271419/91‑1 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑ATCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 >  2:2170394/1‑91 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 <  1:3351752/91‑1 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 >  1:3317422/1‑91 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 <  1:2437860/91‑1 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 <  2:2573860/91‑1 (MQ=255)

                                   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                 <  1:2111956/91‑1 (MQ=255)

                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                               >  2:2083263/1‑91 (MQ=255)

                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                               <  2:2252202/91‑1 (MQ=255)

                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                               <  2:3537070/91‑1 (MQ=255)

                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                               >  1:1648247/1‑91 (MQ=255)

                                      cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAAc                                              >  2:776103/1‑91 (MQ=255)

                                      cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAAc                                              >  1:568087/1‑91 (MQ=255)

                                              tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGa                                      >  2:1476738/1‑91 (MQ=255)

                                              tCACTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGa                                      <  1:2855235/91‑1 (MQ=255)

                                               cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                     <  1:3213273/91‑1 (MQ=255)

                                               cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                     >  1:2671998/1‑91 (MQ=255)

                                                aTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACt                                    <  1:3348233/91‑1 (MQ=255)

                                                aTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACt                                    >  1:2745285/1‑91 (MQ=255)

                                                     ctcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAat                               <  1:361396/91‑1 (MQ=255)

                                                     ctcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAat                               >  2:3439611/1‑91 (MQ=255)

                                                      tcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAatt                              >  2:2106157/1‑91 (MQ=255)

                                                        cTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattca                            <  1:3132715/91‑1 (MQ=255)

                                                             aCTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCaa                       <  1:2957517/91‑1 (MQ=255)

                                                              cTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAat                      <  1:2590193/91‑1 (MQ=255)

                                                               ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                     <  1:3122482/91‑1 (MQ=255)

                                                               ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                     >  1:884481/1‑91 (MQ=255)

                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                   >  2:2434916/1‑91 (MQ=255)

                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                   <  1:2311368/91‑1 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTACAATTCATTCAATATGc                  >  2:2234836/1‑91 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                  <  1:1056009/91‑1 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                  <  1:1495693/91‑1 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                  <  2:2671998/91‑1 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                  >  1:3430830/1‑91 (MQ=255)

                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                  >  2:1579462/1‑91 (MQ=255)

                                                                   aTAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCt                 >  1:1995799/1‑91 (MQ=255)

                                                                    tAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCtt                >  2:3285715/1‑91 (MQ=255)

                                                                    tAAAGCCGAAGCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCtt                <  1:3357161/91‑1 (MQ=255)

                                                                     aaaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa               <  1:2234859/91‑1 (MQ=255)

                                                                      aaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAt              >  1:152613/1‑91 (MQ=255)

                                                                       aGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt             >  2:466379/1‑91 (MQ=255)

                                                                       aGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt             >  2:2050564/1‑91 (MQ=255)

                                                                       aGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt             <  2:2877752/91‑1 (MQ=255)

                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc            >  1:3218579/1‑91 (MQ=255)

                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc            >  2:3630285/1‑91 (MQ=255)

                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAAGCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc            >  2:3190891/1‑91 (MQ=255)

                                                                          cGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGc          >  2:206260/1‑91 (MQ=255)

                                                                          cGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGc          >  2:1331253/1‑91 (MQ=255)

                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa         >  1:291090/1‑91 (MQ=255)

                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa         >  2:2284716/1‑91 (MQ=255)

                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa         >  1:3274186/1‑91 (MQ=255)

                                                                            aaTCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCat        >  2:98779/1‑91 (MQ=255)

                                                                            aaTCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCat        <  1:1065464/91‑1 (MQ=255)

                                                                               gccttttTAGCGTTAAACCATTTAGCTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCatat      <  2:2825247/88‑1 (MQ=255)

                                                                               ccTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt     >  1:1772659/1‑91 (MQ=255)

                                                                               ccTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt     <  1:1174945/91‑1 (MQ=255)

                                                                                 cttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTa    >  2:589959/2‑91 (MQ=255)

                                                                                 cttttTAGCGTTAAACCATTTAACTGCACCTTGTTTCTTAATCTGAAACCCCCAAGTCTCAAATTCATTCAATATGCTTATTCGCATATTa    >  2:1689656/2‑91 (MQ=255)

                                                                                  ttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc   <  2:2589561/91‑1 (MQ=255)

                                                                                  ttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc   <  1:1919221/91‑1 (MQ=255)

TGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACA  >  USA300TCH1516_ALE/1439921‑1440094

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