BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	693,263	G→C	5.4%	P63R (CCT→CGT)	USA300HOU_RS03345 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	693,263	0	G	C	5.4%	843.2 / 25.9	240	P63R (CCT→CGT)	USA300HOU_RS03345	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (108/119); new base C (11/2); total (119/121)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 248 total aligned reads displayed.

ATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGC  >  USA300TCH1516_ALE/693173‑693353

aTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAg                                                                                            >  2:1788415/1‑91 (MQ=255)

  aCAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGc                                                                                          <  1:3554926/91‑1 (MQ=255)

   cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGcc                                                                                         <  1:3728729/91‑1 (MQ=255)

     aTCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGa                                                                                       <  2:3408173/91‑1 (MQ=255)

       cTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGAcc                                                                                     >  1:3773225/1‑91 (MQ=255)

        ttCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCt                                                                                    >  2:712111/1‑91 (MQ=255)

          ccAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCttt                                                                                  >  2:3772841/1‑91 (MQ=255)

           cAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCtttt                                                                                 <  1:328892/91‑1 (MQ=255)

           cAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCtttt                                                                                 <  2:271859/91‑1 (MQ=255)

              aTTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTAcc                                                                              <  2:1345409/91‑1 (MQ=255)

              aTTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTAcc                                                                              <  2:1721380/91‑1 (MQ=255)

               ttCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                             <  1:2368040/91‑1 (MQ=255)

               ttCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCg                                                                             <  1:768561/91‑1 (MQ=255)

                 cACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTa                                                                           <  2:2253114/91‑1 (MQ=255)

                 cACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTa                                                                           >  1:3553590/1‑91 (MQ=255)

                  aCTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAt                                                                          <  2:3565580/91‑1 (MQ=255)

                     tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAtttt                                                                       >  1:722283/1‑91 (MQ=255)

                     tGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAtttt                                                                       >  1:3674013/1‑91 (MQ=255)

                      gTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTc                                                                      >  2:2201061/1‑91 (MQ=255)

                      gTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTc                                                                      >  1:926467/1‑91 (MQ=255)

                         gCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCttt                                                                   <  1:2172518/91‑1 (MQ=255)

                          cATTCTGCCATGTCCTGCTTCACTCACGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTa                                                                  >  2:3177268/1‑91 (MQ=255)

                           aTTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTaa                                                                 <  2:1284608/91‑1 (MQ=255)

                            ttCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAAc                                                                <  1:1666703/91‑1 (MQ=255)

                             tCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACg                                                               <  1:1461483/91‑1 (MQ=255)

                               tGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTa                                                             >  1:3375922/1‑91 (MQ=255)

                                gCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTaa                                                            >  2:3566675/1‑91 (MQ=255)

                                gCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTaa                                                            >  1:1735440/1‑91 (MQ=255)

                                  cATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACt                                                          >  2:2771980/1‑91 (MQ=255)

                                   aTGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTa                                                         <  1:3580369/91‑1 (MQ=255)

                                     gTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATc                                                       >  2:3268483/1‑91 (MQ=255)

                                     gTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATc                                                       >  2:361975/1‑91 (MQ=255)

                                           cTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTaaa                                                 <  1:1788415/91‑1 (MQ=255)

                                           cTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTaaa                                                 >  2:1647715/1‑91 (MQ=255)

                                              cACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAttt                                              <  1:1182666/91‑1 (MQ=255)

                                                 tCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGa                                           >  2:1837168/1‑91 (MQ=255)

                                                  cAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGAc                                          >  2:2577405/1‑91 (MQ=255)

                                                  cAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGAc                                          <  2:1832234/91‑1 (MQ=255)

                                                   aaGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACt                                         >  2:2619768/1‑91 (MQ=255)

                                                   aaGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACt                                         >  2:2075006/1‑91 (MQ=255)

                                                   aaGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACt                                         >  2:489483/1‑91 (MQ=255)

                                                    aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                        <  1:347078/91‑1 (MQ=255)

                                                    aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                        >  1:2057019/1‑91 (MQ=255)

                                                    aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                        >  1:1015550/1‑91 (MQ=255)

                                                    aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                        <  1:3834077/91‑1 (MQ=255)

                                                    aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                                        >  1:451934/1‑91 (MQ=255)

                                                      cTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTa                                      <  2:2407425/91‑1 (MQ=255)

                                                       tCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTAt                                     <  1:2863026/91‑1 (MQ=255)

                                                         aTACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTAttt                                   >  2:1099553/1‑91 (MQ=255)

                                                            cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                                <  1:2618239/91‑1 (MQ=255)

                                                            cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                                <  2:663813/91‑1 (MQ=255)

                                                             acaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCttt                               <  1:759307/91‑1 (MQ=255)

                                                                aTAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGc                            <  2:1964264/91‑1 (MQ=255)

                                                                 tAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCt                           <  1:2214281/91‑1 (MQ=255)

                                                                 tAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTCCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCt                           >  2:3475405/1‑91 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          >  1:2773037/1‑91 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          >  2:3655916/1‑91 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          >  2:3198422/1‑91 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          >  2:2928840/1‑91 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          <  1:2562573/91‑1 (MQ=255)

                                                                  aGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTa                          >  2:2927394/1‑91 (MQ=255)

                                                                   gACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaa                         >  1:2452237/1‑91 (MQ=255)

                                                                   gACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaa                         >  2:727565/1‑91 (MQ=255)

                                                                    taaaaGGTCGGCCTAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaaa                        >  1:1296719/3‑91 (MQ=25)

                                                                    aCAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaaa                        >  2:1303557/1‑91 (MQ=255)

                                                                    aCAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTaaa                        <  1:2471589/91‑1 (MQ=255)

                                                                     cAAAGGTCAGCGTAGATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGGCTGTGTTTCTTTAGCTagag                       <  2:58935/91‑1 (MQ=255)

                                                                     cAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAg                       >  2:3831300/1‑91 (MQ=255)

                                                                     cAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCACTAAATTTAGACTGTATTTCTTTAGCTAAAg                       >  2:399727/1‑91 (MQ=255)

                                                                      aaaGGTCGGCCTAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                      >  1:2059987/1‑91 (MQ=25)

                                                                      aaaGGTCGGCCTAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                      >  2:3027103/1‑91 (MQ=25)

                                                                      aaaGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAga                      >  2:2014077/1‑90 (MQ=255)

                                                                      aaaGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGt                      >  2:797027/1‑91 (MQ=255)

                                                                       aaGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTc                     <  1:1375124/91‑1 (MQ=255)

                                                                        aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa                    <  2:2053932/91‑1 (MQ=255)

                                                                        aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa                    <  2:2745552/91‑1 (MQ=255)

                                                                        aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa                    >  1:372355/1‑91 (MQ=255)

                                                                        aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa                    <  2:1611199/91‑1 (MQ=255)

                                                                        aGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCa                    >  2:1554284/1‑91 (MQ=255)

                                                                         ggTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAt                   <  1:1906598/91‑1 (MQ=255)

                                                                         ggTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAt                   <  2:876636/91‑1 (MQ=255)

                                                                             ggcctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAAt               >  1:787021/5‑91 (MQ=38)

                                                                             ggcctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAAt               >  1:3162765/5‑91 (MQ=38)

                                                                             aGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAAt               >  1:653918/1‑91 (MQ=255)

                                                                              gCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATg              <  1:2928631/91‑1 (MQ=255)

                                                                               cctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGa             >  2:2365033/3‑91 (MQ=255)

                                                                               cGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGa             <  1:1643361/91‑1 (MQ=255)

                                                                               cGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGa             >  1:1342724/1‑91 (MQ=255)

                                                                                ctAAAATTTACGCTGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGaa            >  1:964384/2‑91 (MQ=255)

                                                                                  aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc          <  2:1182341/91‑1 (MQ=255)

                                                                                  aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc          <  1:241656/91‑1 (MQ=255)

                                                                                  aaaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAAc          <  1:1264151/91‑1 (MQ=255)

                                                                                   aaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACt         <  2:1421937/91‑1 (MQ=255)

                                                                                   aaTTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACt         >  1:1576311/1‑91 (MQ=255)

                                                                                     ttttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACttt       >  1:198745/1‑91 (MQ=255)

                                                                                     ttttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACttt       >  1:198624/1‑91 (MQ=255)

                                                                                      tttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTa      <  1:1087531/91‑1 (MQ=255)

                                                                                       ttAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTAt     >  1:838386/1‑91 (MQ=255)

                                                                                          ggCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGc  >  2:115990/1‑91 (MQ=255)

                                                                                          ggCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGc  <  1:1442751/91‑1 (MQ=255)

ATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGAAACTTTATCGC  >  USA300TCH1516_ALE/693173‑693353

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