BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,415,013	C→G	5.1%	intergenic (+24/+24)	yitU_2 → / ← mqo	5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU/malate:quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,415,013	0	C	G	5.1%	423.2 / 9.6	118	intergenic (+24/+24)	yitU_2/mqo	5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (60/52); new base G (3/3); total (63/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 120 total aligned reads displayed.

TAAAATCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCATTTCTCCA  >  USA300TCH1516_ALE/1414926‑1415103

tAAAATCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCAtt                                                                                         >  1:729483/1‑91 (MQ=255)

gaaaaTCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAACTAAAAACACCCCTATCAAATGATAATCAtt                                                                                         <  1:2382154/90‑1 (MQ=255)

   aaTCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATc                                                                                      <  1:901716/91‑1 (MQ=255)

     tCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCaa                                                                                    >  1:3410207/1‑91 (MQ=255)

       aaaCACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCGATTGATAATGATTATCattt                                                                                  <  1:2093084/91‑4 (MQ=38)

       aaaCACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCGATTGATAATGATTATCattt                                                                                  <  1:1586906/91‑4 (MQ=38)

        aaCACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCg                                                                                 <  1:1655097/91‑1 (MQ=255)

         acacATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGa                                                                                >  1:200785/1‑91 (MQ=255)

         acacATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGa                                                                                >  1:1345469/1‑91 (MQ=255)

         acacATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGa                                                                                >  1:3764619/1‑91 (MQ=255)

          cacaTCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGAt                                                                               >  1:2674695/1‑91 (MQ=255)

            caTCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCGATTGATAATGATTATCATTTGATAg                                                                             <  1:77233/91‑1 (MQ=25)

            caTCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAg                                                                             <  1:3943008/91‑1 (MQ=255)

               cAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAgggg                                                                          <  1:132359/91‑1 (MQ=255)

               aaataataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGGTAgggg                                                                          <  1:3219996/90‑1 (MQ=255)

                aataataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGc                                                                         >  1:264850/1‑91 (MQ=255)

                aataataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGc                                                                         >  1:4009102/1‑91 (MQ=255)

                   aataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAt                                                                      <  1:414298/91‑1 (MQ=255)

                   aataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAt                                                                      <  1:414429/91‑1 (MQ=255)

                   aataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAt                                                                      <  1:1139097/91‑1 (MQ=255)

                    ataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAtt                                                                     <  1:3756391/91‑1 (MQ=255)

                    ataatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAtt                                                                     >  1:2932613/1‑91 (MQ=255)

                      aatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAtttt                                                                   >  1:3928233/1‑91 (MQ=255)

                      aatGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTAtttt                                                                   >  1:4035753/1‑91 (MQ=255)

                       atGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTta                                                                  <  1:3028007/91‑1 (MQ=255)

                       atGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTta                                                                  >  1:882282/1‑91 (MQ=255)

                       atGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTta                                                                  <  1:2171326/91‑1 (MQ=255)

                          cacaAGCTATTTGCAAAGTCTTAAAATCATATATATACAAACACCCATATCAAATGATAATCATTATCAATCGATAGGGGCTATTTtaata                                                               >  1:3305776/1‑91 (MQ=255)

                           acaAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTtaataa                                                              <  1:2752538/91‑1 (MQ=255)

                           acaAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTtaataa                                                              >  1:1551202/1‑91 (MQ=255)

                             aaGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATaaaa                                                            <  1:3060299/91‑1 (MQ=255)

                               gCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAAtt                                                          >  1:3892205/1‑91 (MQ=255)

                                  aTTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGt                                                       <  1:1338519/91‑1 (MQ=255)

                                  aTTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGt                                                       <  1:2406344/91‑1 (MQ=255)

                                  aTTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGt                                                       >  1:591669/1‑91 (MQ=255)

                                  aTTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGt                                                       <  1:1339100/91‑1 (MQ=255)

                                   tttGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTc                                                      >  1:1391358/1‑91 (MQ=255)

                                    ttGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTcc                                                     >  1:2568831/1‑91 (MQ=255)

                                     tGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCt                                                    >  1:235384/1‑91 (MQ=255)

                                      gCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTc                                                   <  1:2961623/91‑1 (MQ=255)

                                       cAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCg                                                  <  1:474442/91‑1 (MQ=255)

                                       cAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCg                                                  <  1:474394/91‑1 (MQ=255)

                                        aaaGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGa                                                 >  1:126778/1‑91 (MQ=255)

                                        aaaGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGa                                                 >  1:2103054/1‑91 (MQ=255)

                                         aaGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGaa                                                <  1:1405143/91‑1 (MQ=255)

                                         aaGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGaa                                                <  1:1306735/91‑1 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          <  1:3992324/91‑1 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          <  1:1162443/91‑1 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          >  1:224571/1‑91 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          >  1:2502349/1‑91 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          <  1:2731212/91‑1 (MQ=255)

                                               tAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttc                                          >  1:2502386/1‑91 (MQ=255)

                                                aaaaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttct                                         <  1:3057623/91‑1 (MQ=255)

                                                 aaaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCACTCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttctt                                        >  1:3291247/1‑91 (MQ=255)

                                                 aaaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttctt                                        <  1:2841516/91‑1 (MQ=255)

                                                 aaaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttctt                                        >  1:3136493/1‑91 (MQ=255)

                                                 aaaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttctt                                        <  1:546643/91‑1 (MQ=255)

                                                  aaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttcttc                                       <  1:1497831/91‑1 (MQ=255)

                                                  aaTCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATttcttc                                       <  1:1690007/91‑1 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    >  1:1190941/1‑91 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    >  1:1191416/1‑91 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    >  1:842914/1‑91 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    <  1:3887274/91‑1 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    <  1:2633278/91‑1 (MQ=255)

                                                     cATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctc                                    >  1:2737482/1‑91 (MQ=255)

                                                      atatAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCctct                                   <  1:2187078/91‑1 (MQ=255)

                                                       tataAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCtt                                  >  1:3103295/1‑91 (MQ=255)

                                                       tataAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCtt                                  <  1:2636611/91‑1 (MQ=255)

                                                       tataAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCtt                                  >  1:3123509/1‑91 (MQ=255)

                                                         taaataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCa                                >  1:802122/1‑91 (MQ=255)

                                                         taaataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCa                                >  1:2315568/1‑91 (MQ=255)

                                                           aataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATc                              >  1:2031333/1‑91 (MQ=255)

                                                           aataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATc                              <  1:1382307/91‑1 (MQ=255)

                                                           aataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATc                              <  1:1382617/91‑1 (MQ=255)

                                                            ataaaAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCt                             <  1:4107904/91‑1 (MQ=255)

                                                              aaaaaCACCCCTAGCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTaa                           <  1:4151744/91‑1 (MQ=255)

                                                               aaaaCACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAAt                          >  1:477322/1‑91 (MQ=255)

                                                                aaaCACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATc                         <  1:4100045/91‑1 (MQ=255)

                                                                aaaCACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATc                         >  1:3374223/1‑91 (MQ=255)

                                                                aaaCACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATc                         <  1:3441078/91‑1 (MQ=255)

                                                                  acacCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAGTTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCa                       <  1:3511518/91‑1 (MQ=255)

                                                                  acacCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCa                       <  1:3511469/91‑1 (MQ=255)

                                                                   cacCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCaa                      >  1:2988546/1‑91 (MQ=255)

                                                                   cacCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCaa                      >  1:2203667/1‑91 (MQ=255)

                                                                    acCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCaaa                     <  1:380834/91‑1 (MQ=255)

                                                                     ccccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAaat                    >  1:3494670/1‑91 (MQ=255)

                                                                     ccccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAaat                    <  1:1793884/91‑1 (MQ=255)

                                                                     ccccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAaat                    >  1:3495659/1‑91 (MQ=255)

                                                                     ccccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCAAATCCAaat                    >  1:3976042/1‑91 (MQ=255)

                                                                      cccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAaata                   >  1:891845/1‑91 (MQ=255)

                                                                      cccTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAaata                   >  1:891217/1‑91 (MQ=255)

                                                                             gattGATAATGATTATCATTTGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGc            >  1:283733/4‑91 (MQ=38)

                                                                             aaaTGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCGAAATAATTCTGc            <  1:3188498/91‑1 (MQ=255)

                                                                              attGATAATGATTATCATTTGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGcc           >  1:3036833/3‑91 (MQ=38)

                                                                                 gATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCAtt        >  1:3422063/1‑91 (MQ=255)

                                                                                 gATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCAtt        >  1:3417421/1‑91 (MQ=255)

                                                                                      tCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCACATAATTCTGCCGTTTCTcc   >  1:2082250/1‑91 (MQ=255)

                                                                                      tCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCATTTCTcc   >  1:3426264/1‑91 (MQ=255)

                                                                                       cATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCATTTCTCCa  >  1:4047051/1‑91 (MQ=255)

                                                                                       cATTATCAATCGATAGGGGCTATTTTAATAAAATTCGACCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCATTTCTCCa  >  1:3494998/1‑91 (MQ=255)

TAAAATCAAACACATCAATAATAATGCACAAGCTATTTGCAAAGTCTTAAAATCATATAAATAAAAACACCCCTATCAAATGATAATCATTATCAATCGATAGGGGCTATTTTAATAAAATTCGTCCTCGAACATTTCTTCCTCTTCATCTAATCCAAATAATTCTGCCATTTCTCCA  >  USA300TCH1516_ALE/1414926‑1415103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`