BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,646,189	G→C	A155G (GCA→GGA)	ypdF ←	Aminopeptidase YpdF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,646,189	0	G	C	99.5%	696.0 / NA	191	A155G (GCA→GGA)	ypdF	Aminopeptidase YpdF
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (89/101); minor base T (1/0); total (90/101)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

CAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTG  >  USA300TCH1516_ALE/1646103‑1646275

cAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                    >  2:2570297/1‑91 (MQ=255)

 aaTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt                                                                                   <  1:384115/91‑1 (MQ=255)

   tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                 <  2:100874/91‑1 (MQ=255)

   tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                 <  2:800843/91‑1 (MQ=255)

     gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                               >  1:2330025/1‑91 (MQ=255)

     gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                               >  1:178125/1‑91 (MQ=255)

     gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                               <  2:2428055/91‑1 (MQ=255)

       aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaaca                                                                             >  1:229581/1‑91 (MQ=255)

        tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                            >  1:1363973/1‑91 (MQ=255)

        tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                            <  2:3282550/91‑1 (MQ=255)

         cGAATGATGGTCCATCTGCTCCTAATTCTAGCATTGTGCTTTCTATTATTCCGTTTAATTGTTTTTCAGTCATGCCTCCTTTTATAACAGt                                                                           >  2:2689612/1‑91 (MQ=255)

           aaTGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTa                                                                         <  2:2224455/91‑1 (MQ=255)

           aaTGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTa                                                                         <  2:2224448/91‑1 (MQ=255)

            atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaa                                                                        <  1:2588673/91‑1 (MQ=255)

            atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaa                                                                        <  1:1617722/91‑1 (MQ=255)

              gatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                      <  2:3207131/91‑1 (MQ=255)

              gatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                      >  2:1167224/1‑91 (MQ=255)

              gatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCACGCCTCCTTTTACAACAGTTaaaa                                                                      <  2:2185183/91‑1 (MQ=255)

               atgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAat                                                                     >  2:2127613/1‑91 (MQ=255)

                tgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                    <  1:1098432/91‑1 (MQ=255)

                   tCCATCTGCTCCTAATTCTAGCATTTTGCTTTTTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                 >  2:376358/1‑91 (MQ=255)

                   tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                 >  2:3355923/1‑91 (MQ=255)

                   tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                 >  1:340781/1‑91 (MQ=255)

                   tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                 <  1:1010379/91‑1 (MQ=255)

                    ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATAtt                                                                >  2:1122542/1‑91 (MQ=255)

                    ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATAtt                                                                >  2:1148885/1‑91 (MQ=255)

                    ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATAtt                                                                >  2:501006/1‑91 (MQ=255)

                     cATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTc                                                               <  2:1390342/91‑1 (MQ=255)

                       tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                             <  2:3390046/91‑1 (MQ=255)

                       tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                             <  1:1350930/91‑1 (MQ=255)

                       tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                             <  2:2505343/91‑1 (MQ=255)

                       tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                             <  2:1853290/91‑1 (MQ=255)

                         tGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCatat                                                           >  1:129064/1‑91 (MQ=255)

                          gCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATg                                                          <  1:769151/91‑1 (MQ=255)

                           cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         >  2:2178289/1‑91 (MQ=255)

                           cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         >  1:1328831/1‑91 (MQ=255)

                           cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         >  1:210408/1‑91 (MQ=255)

                           cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         >  1:211806/1‑91 (MQ=255)

                            tCCTAATTCTAGCATTTTGCTTTCTAATATTCCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAGCAGTTAAAATACATTCATATGtt                                                        >  2:2774771/1‑91 (MQ=255)

                             ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGttt                                                       >  2:1627391/1‑91 (MQ=255)

                             ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGttt                                                       >  1:3236542/1‑91 (MQ=255)

                             ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGttt                                                       >  1:995615/1‑91 (MQ=255)

                                aaTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcat                                                    <  1:830790/91‑1 (MQ=255)

                                 aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                   <  1:2542566/91‑1 (MQ=255)

                                  ttCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                  <  1:891343/91‑1 (MQ=255)

                                  ttCTAGCATTTTGCTTTCTAATATTGCCTTTAAGTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                  <  1:573622/91‑1 (MQ=255)

                                    cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAGAATATATTCATATGTTTCATcaac                                                <  2:543605/91‑1 (MQ=255)

                                    cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaac                                                >  1:3118071/1‑91 (MQ=255)

                                    cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaac                                                <  1:1646909/91‑1 (MQ=255)

                                    cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaac                                                >  1:912416/1‑91 (MQ=255)

                                     tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaaca                                               <  2:2424478/91‑1 (MQ=255)

                                      aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaacaa                                              <  1:3204609/91‑1 (MQ=255)

                                       gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACCACAGTTAAAATATATTCATATGTTTCATCAACAat                                             >  2:2640148/1‑91 (MQ=255)

                                       gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAat                                             >  1:330433/1‑91 (MQ=255)

                                       gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAat                                             >  1:153466/1‑91 (MQ=255)

                                        cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAata                                            >  2:83132/1‑91 (MQ=255)

                                          ttttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATAtt                                          <  1:2031749/91‑1 (MQ=255)

                                             tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                       >  1:1583203/1‑91 (MQ=255)

                                             tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                       >  1:1879096/1‑91 (MQ=255)

                                             tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                       <  2:2933158/91‑1 (MQ=255)

                                              gCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgct                                      >  1:1100093/1‑91 (MQ=255)

                                               cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctg                                     >  2:2284813/1‑91 (MQ=255)

                                               cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctg                                     >  2:2836556/1‑91 (MQ=255)

                                                tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                    <  2:968915/91‑1 (MQ=255)

                                                 ttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                   >  2:2913094/1‑91 (MQ=255)

                                                 ttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                   <  2:788190/91‑1 (MQ=255)

                                                 ttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                   <  2:456182/91‑1 (MQ=255)

                                                   cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttt                                 <  2:232027/91‑1 (MQ=255)

                                                    tAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                <  1:3342560/91‑1 (MQ=255)

                                                       tatTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGa                             <  2:1691088/91‑1 (MQ=255)

                                                       tatTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGa                             >  1:1014519/1‑91 (MQ=255)

                                                         ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCCTCAACAATATTAGCTGCTTTTTGAAt                           >  2:1331007/1‑91 (MQ=255)

                                                         ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAt                           >  2:3320051/1‑91 (MQ=255)

                                                           gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                         >  2:2559038/1‑91 (MQ=255)

                                                           gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                         <  1:1867881/91‑1 (MQ=255)

                                                            ccTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                        <  1:501006/91‑1 (MQ=255)

                                                             cTTTAATTCGTTCTCAGTCAGGCCTCCTATTACAACCGTTAAGATATATTCATGTGTGTCATCTGCAATATTAGCTGCTTTTTGAGTTcaa                       <  2:1498902/91‑1 (MQ=255)

                                                              tttAATTCTTTTTTTTTTTTTCTTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                      >  2:2963770/1‑91 (MQ=255)

                                                              tttAATTCTTTTTCAGTCCTGCCTCCTTTTACAACAGTGAAAATATATTCATATGTTTCATCAACAAGATTAGCTGCTTTTTGAATTAAAg                      <  1:2108277/91‑1 (MQ=255)

                                                              tttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                      <  1:3538574/91‑1 (MQ=255)

                                                              tttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                      <  1:1180244/91‑1 (MQ=255)

                                                               ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                     <  2:354680/91‑1 (MQ=255)

                                                               ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                     <  2:2198106/91‑1 (MQ=255)

                                                                tAATTCTTTTTCAGTCATGCCTCCTATTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                    <  1:3376605/91‑1 (MQ=255)

                                                                     cTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTc               >  1:2178655/1‑91 (MQ=255)

                                                                      tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg              <  2:882239/91‑1 (MQ=255)

                                                                      tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg              >  2:1396931/1‑91 (MQ=255)

                                                                      tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg              <  1:2965844/91‑1 (MQ=255)

                                                                      tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTCCg              >  2:1549527/1‑91 (MQ=255)

                                                                       ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGt             >  2:3093362/1‑91 (MQ=255)

                                                                          tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAg          >  1:736477/1‑91 (MQ=255)

                                                                          tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGACAg          >  1:2920729/1‑91 (MQ=255)

                                                                           cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc         <  2:1284874/91‑1 (MQ=255)

                                                                             gTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCAt       <  2:3590839/91‑1 (MQ=255)

                                                                             gTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCAt       >  2:422022/1‑91 (MQ=255)

                                                                             gTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACAt       >  1:2812376/1‑91 (MQ=255)

                                                                              tCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCAt       >  1:1248407/1‑90 (MQ=255)

                                                                              tCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATc      <  1:1911399/91‑1 (MQ=255)

                                                                               cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt     >  1:1846519/1‑91 (MQ=255)

                                                                               cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt     <  1:1299070/91‑1 (MQ=255)

CAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTG  >  USA300TCH1516_ALE/1646103‑1646275

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