BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	693,253	G→C	5.2%	Y66* (TAC→TAG)	USA300HOU_RS03345 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	693,253	0	G	C	5.2%	904.6 / 28.7	252	Y66* (TAC→TAG)	USA300HOU_RS03345	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (113/126); new base C (4/9); total (117/135)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 259 total aligned reads displayed.

GCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGA  >  USA300TCH1516_ALE/693163‑693342

tccAATATAAATACAATCTTCCAAATTCACTTGTAGGATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCg                                                                                           <  2:3837346/90‑1 (MQ=255)

gCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCg                                                                                           <  2:1147257/91‑1 (MQ=255)

gCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCg                                                                                           <  1:84856/91‑1 (MQ=255)

gCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCg                                                                                           >  2:1869565/1‑91 (MQ=255)

 ccAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGt                                                                                          <  2:1303034/91‑1 (MQ=255)

  cAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTa                                                                                         <  1:2566690/91‑1 (MQ=255)

  cAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTa                                                                                         <  2:2800370/91‑1 (MQ=255)

   aaTATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTaa                                                                                        <  2:2168565/91‑1 (MQ=255)

   aaTATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTaa                                                                                        <  1:1581020/91‑1 (MQ=255)

   aaTATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTaa                                                                                        >  2:1081859/1‑91 (MQ=255)

      ataAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTaaa                                                                                       <  2:3210424/89‑1 (MQ=255)

      ataAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTaaa                                                                                       >  1:3210424/1‑89 (MQ=255)

       taAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAttt                                                                                    >  1:271859/1‑91 (MQ=255)

       taAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAttt                                                                                    <  2:2504923/91‑1 (MQ=255)

        aaaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAAtttt                                                                                   >  1:2407425/1‑91 (MQ=255)

         aaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTa                                                                                  >  1:1324984/1‑91 (MQ=255)

         aaTACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTa                                                                                  >  1:2329298/1‑91 (MQ=255)

           tACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAgg                                                                                <  1:3673624/91‑1 (MQ=255)

             cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTacgct                                                                              <  1:2106360/91‑5 (MQ=38)

             cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTacgct                                                                              <  2:1193985/91‑5 (MQ=38)

             cAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCGGCCTAAAATTTACGct                                                                              <  1:3572758/91‑2 (MQ=38)

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

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

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

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

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

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

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

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

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

                       aaTTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTAc                                                                    >  2:1392963/1‑91 (MQ=255)

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

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

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

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

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

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

                              ttGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAttt                                                             <  2:1599821/91‑1 (MQ=255)

                              ttGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTAttt                                                             >  2:2214281/1‑91 (MQ=255)

                                 tAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCTTACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCt                                                          >  2:2286330/1‑91 (MQ=255)

                                  aGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCtt                                                         >  1:1365497/1‑91 (MQ=255)

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

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

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

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

                                      ttCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAAc                                                     >  1:2183236/1‑91 (MQ=255)

                                       tCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACg                                                    >  1:1740912/1‑91 (MQ=255)

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

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

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

                                              tGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTAt                                             <  2:841112/91‑1 (MQ=255)

                                               gTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACttcc                                            >  2:3586224/1‑88 (MQ=255)

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

                                                        cACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAAttt                                   >  1:3289887/1‑91 (MQ=255)

                                                         aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTa                                  >  1:794077/1‑91 (MQ=255)

                                                         aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTa                                  >  1:715487/1‑91 (MQ=255)

                                                         aCTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTa                                  >  2:1694510/1‑91 (MQ=255)

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

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

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

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

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

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

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

                                                              aGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTg                             >  1:3497651/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)

                                                                  cATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTAtt                         >  1:730970/1‑91 (MQ=255)

                                                                      cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                     >  2:3577094/1‑91 (MQ=255)

                                                                      cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                     >  1:1422644/1‑91 (MQ=255)

                                                                      cacaATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCtt                     <  1:939703/91‑1 (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               >  1:2773037/1‑91 (MQ=255)

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

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

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

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

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

                                                                              aCAAAGGTCAGCGTAAA‑TTTAGGCCGACCTTTGACCGGATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAg            >  2:3085907/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                    gTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCAtt       >  1:2154386/1‑91 (MQ=255)

                                                                                      cAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa     <  1:1648135/91‑1 (MQ=255)

                                                                                      cAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTaa     >  2:3428305/1‑91 (MQ=255)

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

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

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

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

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

GCCAATATAAATACAATCTTCCAAATTCACTTGTAGCATTCTGCCATGTCCTGCTTCACTCAAGCTCATACACAATAGACAAAGGTCAGCGTAAATTTTAGGCCGACCTTTTACCGTATTTTCTTTAACGTAACTATCATTAAATTTAGACTGTATTTCTTTAGCTAAAGTCATTAATGA  >  USA300TCH1516_ALE/693163‑693342

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