BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,905,775	G→C	12.3%	E128Q (GAA→CAA)	sigS →	RNA polymerase sigma factor SigS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,905,775	0	G	C	12.3%	442.0 / 46.3	138	E128Q (GAA→CAA)	sigS	RNA polymerase sigma factor SigS
Reads supporting (aligned to +/- strand): ref base G (69/52); new base C (6/11); total (75/63)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 144 total aligned reads displayed.

CATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGT  >  USA300TCH1516_ALE/1905685‑1905865

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGACGGATACAAGCAATTTg                                                                                            >  1:379453/1‑91 (MQ=255)

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTg                                                                                            >  1:3389602/1‑91 (MQ=255)

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTg                                                                                            >  1:3389255/1‑91 (MQ=255)

cATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTg                                                                                            >  1:2673234/1‑91 (MQ=255)

   gACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCaaa                                                                                         <  1:2764059/91‑1 (MQ=255)

   gACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCaaa                                                                                         <  1:2764002/91‑1 (MQ=255)

   gACCTTCGTTTACAAGATATCTTCAAGCTGTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCaaa                                                                                         <  1:2764060/91‑1 (MQ=255)

    aCCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAt                                                                                        >  1:4064945/1‑91 (MQ=255)

    aCCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCAAAt                                                                                        <  1:253036/91‑1 (MQ=255)

    aCCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTCAAAt                                                                                        <  1:523825/91‑1 (MQ=255)

     ccTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGGGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAtt                                                                                       <  1:368587/91‑1 (MQ=255)

     ccTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAtt                                                                                       <  1:108222/91‑1 (MQ=255)

     ccTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAAtt                                                                                       <  1:1294544/91‑1 (MQ=255)

        tCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCt                                                                                    >  1:113600/1‑91 (MQ=255)

         cGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTg                                                                                   <  1:367583/91‑1 (MQ=255)

          gTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGa                                                                                  <  1:1267306/91‑1 (MQ=255)

          gTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGa                                                                                  >  1:1746097/1‑91 (MQ=255)

            ttACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGaaa                                                                                >  1:1631285/1‑91 (MQ=255)

              aCAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATc                                                                              <  1:2251179/91‑1 (MQ=255)

               cAAGATATCTTCTAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCa                                                                             >  1:3186770/1‑91 (MQ=255)

                aaGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCAt                                                                            <  1:3073492/91‑1 (MQ=255)

                 aGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATg                                                                           <  1:2604795/91‑1 (MQ=255)

                    tatCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGACATCATGTCa                                                                        >  1:3974842/1‑91 (MQ=255)

                    tatCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCa                                                                        >  1:2606950/1‑91 (MQ=255)

                     atCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATTTCAt                                                                       <  1:3272517/91‑1 (MQ=255)

                      tCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCAtt                                                                      <  1:1641719/91‑1 (MQ=255)

                       cTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTa                                                                     >  1:2643981/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:2106498/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    <  1:2439666/91‑1 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:1314806/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:4137905/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:3404760/1‑91 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    <  1:3752551/91‑1 (MQ=255)

                        ttCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTAt                                                                    >  1:3404852/1‑91 (MQ=255)

                         tCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATc                                                                   <  1:1174475/91‑1 (MQ=255)

                            aGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCt                                                                >  1:3541150/1‑91 (MQ=255)

                            aGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCt                                                                >  1:3735890/1‑91 (MQ=255)

                                tttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCa                                                            >  1:3762845/1‑91 (MQ=255)

                                tttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCa                                                            <  1:2712606/91‑1 (MQ=255)

                                tttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCa                                                            >  1:440968/1‑91 (MQ=255)

                                 ttAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCaa                                                           >  1:177610/1‑91 (MQ=255)

                                  tAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAAc                                                          <  1:452941/91‑1 (MQ=255)

                                  tAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAACAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAAc                                                          >  1:3525228/1‑91 (MQ=255)

                                    aaTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGa                                                        <  1:321981/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       <  1:2728612/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       <  1:40316/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       <  1:1228216/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       <  1:1228020/91‑1 (MQ=255)

                                     aTCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGat                                                       >  1:3405484/1‑91 (MQ=255)

                                       cAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGatta                                                     >  1:3880227/1‑91 (MQ=255)

                                        aaaGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaa                                                    >  1:1760279/1‑91 (MQ=255)

                                         aaGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaa                                                   <  1:4106664/91‑1 (MQ=255)

                                          agagAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaat                                                  >  1:2260656/1‑91 (MQ=255)

                                           gagaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatt                                                 <  1:1541919/91‑1 (MQ=255)

                                           gagaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatt                                                 <  1:1658089/91‑1 (MQ=255)

                                            agaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaatta                                                >  1:1330271/1‑91 (MQ=255)

                                             gaAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGattaaattaa                                               <  1:1472376/91‑1 (MQ=255)

                                                aGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAGATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAattaatta                                            <  1:3294265/91‑1 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           >  1:2818562/1‑91 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           >  1:56763/1‑91 (MQ=255)

                                                 gACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAAttaattaa                                           >  1:61033/1‑91 (MQ=255)

                                                  aCTATGGCTCAAACTATACCTTGAAGGATGCAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAg                                          <  1:3330254/91‑1 (MQ=255)

                                                   cTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGa                                         <  1:1889198/91‑1 (MQ=255)

                                                       ggCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTc                                     >  1:876322/1‑91 (MQ=255)

                                                       ggCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTc                                     >  1:876593/1‑91 (MQ=255)

                                                        gCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTcc                                    <  1:484554/91‑1 (MQ=255)

                                                        gCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTcc                                    <  1:1442620/91‑1 (MQ=255)

                                                        gCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTcc                                    >  1:3131390/1‑91 (MQ=255)

                                                         cTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCg                                   >  1:955314/1‑91 (MQ=255)

                                                          tCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGt                                  >  1:3654171/1‑91 (MQ=255)

                                                            aaaCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTa                                <  1:1919241/91‑1 (MQ=255)

                                                              aCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAg                              <  1:1090944/91‑1 (MQ=255)

                                                              aCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAg                              <  1:216699/91‑1 (MQ=255)

                                                              aCTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAg                              <  1:3647456/91‑1 (MQ=255)

                                                               cTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGc                             >  1:2940961/1‑91 (MQ=255)

                                                                      ttGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATg                      <  1:2966658/91‑1 (MQ=255)

                                                                       tGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGc                     <  1:2183972/91‑1 (MQ=255)

                                                                        gAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGcc                    >  1:1743433/1‑91 (MQ=255)

                                                                        gAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGcc                    >  1:3337802/1‑91 (MQ=255)

                                                                           ggATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAAc                 >  1:1795503/1‑91 (MQ=255)

                                                                            gATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTACATGCCAACa                >  1:3967679/1‑91 (MQ=255)

                                                                            gATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACa                >  1:894724/1‑91 (MQ=255)

                                                                             aTACAGGCAATTTGAAATTGCTGAACTCATGTCATTATCGCTTTCAACGATTAGATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACAt               <  1:3859347/91‑1 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:4010225/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:1950592/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:1951279/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:4010140/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:2925399/1‑91 (MQ=255)

                                                                              tACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATa              >  1:2925439/1‑91 (MQ=255)

                                                                               aCAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATaa             >  1:710283/1‑91 (MQ=255)

                                                                               aCAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATaa             >  1:2702379/1‑91 (MQ=255)

                                                                               aCAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATaa             >  1:3154334/1‑91 (MQ=255)

                                                                               aCAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATaa             >  1:1407196/1‑91 (MQ=255)

                                                                                 caGCAATTTCAAATTGCTTAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAAtt           >  1:3715220/2‑91 (MQ=255)

                                                                                  aGCAATTTCAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAAttt          >  1:1540556/1‑91 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  1:1712638/91‑1 (MQ=255)

                                                                                      aTTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      <  1:2419085/91‑1 (MQ=255)

                                                                                      aTTTCAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAAt      >  1:1093091/1‑91 (MQ=255)

                                                                                        ttGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTa    >  1:1912343/1‑91 (MQ=255)

                                                                                          gAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGt  >  1:4078809/1‑91 (MQ=255)

CATGACCTTCGTTTACAAGATATCTTCAAGCTTTTAAATCAAAGAGAAAGACTATGGCTCAAACTATACCTTGAAGGATACAAGCAATTTGAAATTGCTGAAATCATGTCATTATCGCTTTCAACGATTAAATTAATTAAGATGTCCGTTAAGCGTAAATGCCAACATAATTTTAATTAGT  >  USA300TCH1516_ALE/1905685‑1905865

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