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.7%	1311.0 / ‑2.4	350	coding (464/1062 nt)	ypdF	Aminopeptidase YpdF
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (197/152); minor base . (0/1); total (197/153)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 355 total aligned reads displayed.

TACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTG  >  USA300TCH1516_ALE/1646101‑1646275

 aCAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtt                                                                                     >  2:3183359/1‑91 (MQ=255)

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

   aaTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt                                                                                   >  1:4236651/1‑91 (MQ=255)

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

      cGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTaca                                                                                <  2:507040/91‑1 (MQ=255)

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

        tATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                              >  2:609511/1‑91 (MQ=255)

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

          tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAGTATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                            >  2:6235334/1‑91 (MQ=255)

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

            ggatgctgGTCCCTCTCCCGGTAGTTCTAGCGTTTTGCTTTCTAATATTGCCTATAATTCTTTTTCCGTCGTGCCACCTTTTACAACAGtt                                                                          >  2:2695428/3‑91 (MQ=255)

            gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGtt                                                                          >  2:3814709/1‑91 (MQ=255)

            gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGtt                                                                          >  1:2231567/1‑91 (MQ=255)

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

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

               tgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaa                                                                       <  2:4150486/91‑1 (MQ=255)

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

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

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

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

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

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

                      ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCCTTAATTCTTTTTCAGTCATGCATCCTTTTACAACAGTTACAATATAtt                                                                >  2:6837082/1‑91 (MQ=255)

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

                        aTCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCa                                                              <  1:2782999/91‑1 (MQ=255)

                        aTCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCa                                                              >  2:202857/1‑91 (MQ=255)

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

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

                            gCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTATAATATATTCATATg                                                          >  2:5705628/1‑91 (MQ=255)

                             cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTACAATATATTCATATGt                                                         >  2:3883855/1‑91 (MQ=255)

                             cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         <  2:1314454/91‑1 (MQ=255)

                             cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         <  2:1578367/91‑1 (MQ=255)

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

                              tCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGtt                                                        <  1:6796421/91‑1 (MQ=255)

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

                               ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATGCATATGttt                                                       >  2:6658795/1‑91 (MQ=255)

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

                                   aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                         <  2:674809/85‑1 (MQ=255)

                                   aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                   >  2:2198555/1‑91 (MQ=255)

                                    ttCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                  >  2:2082196/1‑91 (MQ=255)

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

                                     tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaa                                                 <  2:5260858/91‑1 (MQ=255)

                                     tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaa                                                 <  1:3592070/91‑1 (MQ=255)

                                     tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaa                                                 <  1:4039933/91‑1 (MQ=255)

                                     tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaa                                                 <  2:3467105/91‑1 (MQ=255)

                                        tgCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaacaa                                              <  2:3996383/90‑1 (MQ=255)

                                        aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaacaa                                              >  1:1520018/1‑91 (MQ=255)

                                         gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTGCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAat                                             >  1:4972712/1‑91 (MQ=255)

                                            ttttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATAtt                                          >  2:4274546/1‑91 (MQ=255)

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

                                             tttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa                                         <  2:2327561/91‑1 (MQ=255)

                                             tttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa                                         >  2:3541439/1‑91 (MQ=255)

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

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

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

                                                 cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctg                                     <  1:6162994/91‑1 (MQ=255)

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

                                                    tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAGAATATATTCATATGTTTCATCAACAATATTTGCTGCtt                                  <  1:6755405/91‑1 (MQ=255)

                                                    tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtt                                  <  2:4295449/91‑1 (MQ=255)

                                                    tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATAGTAGCTGCtt                                  >  1:4582213/1‑91 (MQ=255)

                                                     cTAATATTGCGTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttt                                 <  1:7005541/91‑1 (MQ=255)

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

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

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

                                                        atatTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTg                              >  2:1225216/1‑91 (MQ=255)

                                                        atatTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTg                              <  1:4549512/91‑1 (MQ=255)

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

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

                                                           ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAt                           <  2:1646570/91‑1 (MQ=255)

                                                            tGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAtt                          <  1:2904216/91‑1 (MQ=255)

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

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

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

                                                                tttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                      >  2:3285014/1‑91 (MQ=255)

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

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

                                                                  tAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                    <  1:4328515/91‑1 (MQ=255)

                                                                  tAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                    <  2:3047536/91‑1 (MQ=255)

                                                                   aaTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCaa                   <  2:3672906/91‑1 (MQ=255)

                                                                     ttCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAtt                 >  1:2692211/1‑91 (MQ=255)

                                                                      tCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTAATTAAAGCAAttt                >  2:644604/1‑91 (MQ=255)

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

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

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

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

                                                                         ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTGCATCAACAATATTAGCTGCTTTTGGAATTAAAGCAATTTCGt             >  1:3825032/1‑91 (MQ=255)

                                                                           ttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCa           >  1:3377312/1‑91 (MQ=255)

                                                                            tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTAATTAAAGCAATTTCGTCAg          <  1:6405410/91‑1 (MQ=255)

                                                                             cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTATTGCCATTTCGTCAGc         >  2:7020872/1‑91 (MQ=255)

                                                                             cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc         >  2:4391314/1‑91 (MQ=255)

                                                                             cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc         >  1:4791870/1‑91 (MQ=255)

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

                                                                              aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGAAATTTCGTCAGCa        >  2:2664828/1‑91 (MQ=255)

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

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

                                                                                tCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATc      >  2:1138589/1‑91 (MQ=255)

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

                                                                                  aTGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCtt    <  2:178916/91‑1 (MQ=255)

                                                                                   tGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCttt   >  2:744143/1‑91 (MQ=255)

                                                                                   tGCCTCCTTTTACAACAGTTAAAATATATTCATATATTTCATCAAGAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCttt   >  2:6779947/1‑91 (MQ=255)

TACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTG  >  USA300TCH1516_ALE/1646101‑1646275

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