BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,360,384	T→A	5.5%	H233Q (CAT→CAA)	thrB_2 →	Homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,360,384	0	T	A	5.5%	797.5 / 24.4	235	H233Q (CAT→CAA)	thrB_2	Homoserine kinase
Reads supporting (aligned to +/- strand): ref base T (107/115); new base A (5/8); total (112/123)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 246 total aligned reads displayed.

TACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTATTTTA  >  USA300TCH1516_ALE/1360294‑1360474

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            >  1:2552760/1‑91 (MQ=255)

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            <  1:1715724/91‑1 (MQ=255)

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            <  1:446195/91‑1 (MQ=255)

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            <  1:922700/91‑1 (MQ=255)

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            <  2:1949634/91‑1 (MQ=255)

tACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAt                                                                                            <  2:1789044/91‑1 (MQ=255)

 aCAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAtt                                                                                           <  1:3332033/91‑1 (MQ=255)

   aaTGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTa                                                                                         <  2:967126/91‑1 (MQ=255)

   aaTGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTa                                                                                         <  1:2084465/91‑1 (MQ=255)

   aaTGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTa                                                                                         <  2:3850953/91‑1 (MQ=255)

   aaTGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTa                                                                                         <  2:691566/91‑1 (MQ=255)

     tGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAAt                                                                                       >  2:115971/1‑91 (MQ=255)

      gATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAAtt                                                                                      >  2:920013/1‑91 (MQ=255)

      gATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTaaat                                                                                      <  1:3233044/91‑3 (MQ=255)

       aTTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTg                                                                                     <  2:519123/91‑1 (MQ=255)

        tttGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGc                                                                                    <  1:1454348/91‑1 (MQ=255)

        tttGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGc                                                                                    >  1:3458735/1‑91 (MQ=255)

        tttGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCAATTAAATGc                                                                                    <  2:3160613/91‑1 (MQ=255)

         ttGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCt                                                                                   <  1:2658613/91‑1 (MQ=255)

         ttGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCt                                                                                   >  1:860046/1‑91 (MQ=255)

          tgtgCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTg                                                                                  <  2:3671794/91‑1 (MQ=255)

           gtgCATTAGCACAGCACAATTATGAATTAGCAGGTAACCTCATGCAACAAGCTGGCTTTCATGAACCGTATCGTCAGCATTTAATCGCTGa                                                                                 >  1:3608900/1‑91 (MQ=255)

             gCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAAt                                                                               <  1:1668566/91‑1 (MQ=255)

              cATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAAtt                                                                              <  2:1863391/91‑1 (MQ=255)

              cATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAAtt                                                                              <  2:1863675/91‑1 (MQ=255)

              cATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAAtt                                                                              <  1:2936402/91‑1 (MQ=255)

              cATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAAtt                                                                              <  1:3022770/91‑1 (MQ=255)

                  agcacagcacAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTtgat                                                                          >  1:1097745/1‑91 (MQ=255)

                     acagcacaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGaa                                                                       >  1:2387152/1‑91 (MQ=255)

                       agcacaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGt                                                                     <  1:115971/91‑1 (MQ=255)

                       agcacaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGt                                                                     >  2:2553725/1‑91 (MQ=255)

                        gcacaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTg                                                                    >  2:80110/1‑91 (MQ=255)

                         cacaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGa                                                                   <  2:2650188/91‑1 (MQ=255)

                          acaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGaa                                                                  <  1:691089/91‑1 (MQ=255)

                          acaATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGaa                                                                  <  1:1076139/91‑1 (MQ=255)

                           caATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGaaa                                                                 <  2:1203743/91‑1 (MQ=255)

                             aTTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAAc                                                               <  1:2312361/91‑1 (MQ=255)

                              ttATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACa                                                              <  2:1909441/91‑1 (MQ=255)

                               tATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACaa                                                             <  1:2214194/91‑1 (MQ=255)

                               tATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACaa                                                             <  2:1407338/91‑1 (MQ=255)

                                aTGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAAt                                                            <  2:1228410/91‑1 (MQ=255)

                                 tGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAAtt                                                           <  2:2778981/91‑1 (MQ=255)

                                  gAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTg                                                          <  1:1023747/91‑1 (MQ=255)

                                  gAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTg                                                          >  1:1385083/1‑91 (MQ=255)

                                    aTTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCt                                                        <  2:2673853/91‑1 (MQ=255)

                                     ttAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAtttaaaa                                                       >  2:3285868/1‑87 (MQ=255)

                                     ttAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTa                                                       >  1:1221627/1‑91 (MQ=255)

                                      tAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAg                                                      >  2:1374095/1‑91 (MQ=255)

                                       aGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGt                                                     >  1:538165/1‑91 (MQ=255)

                                         cAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCa                                                   <  2:2545948/91‑1 (MQ=255)

                                         cAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCa                                                   <  1:2181307/91‑1 (MQ=255)

                                         cAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCa                                                   <  1:1590202/91‑1 (MQ=255)

                                           ggTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAAc                                                 >  1:575514/1‑91 (MQ=255)

                                            gTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACa                                                >  1:1911409/1‑91 (MQ=255)

                                              aaaCTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATa                                              >  1:961097/1‑91 (MQ=255)

                                                 cTCATGCAACACGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGAGGAAGTGAAAACAATTGCTAGTCAACATAATg                                           <  2:2961758/91‑1 (MQ=255)

                                                 cTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATg                                           <  2:333121/91‑1 (MQ=255)

                                                 cTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAAAAATTGCTAGTCAACATAATg                                           >  2:2522516/1‑91 (MQ=255)

                                                    aTGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCt                                        <  1:370898/91‑1 (MQ=255)

                                                      gCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTAt                                      >  1:1452141/1‑91 (MQ=255)

                                                       caacaaGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATg                                     >  1:1305291/1‑91 (MQ=255)

                                                       caacaaGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATg                                     <  2:2289307/91‑1 (MQ=255)

                                                          caaGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCaa                                  >  2:3143329/1‑91 (MQ=255)

                                                             gATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTg                               >  2:3523200/1‑91 (MQ=255)

                                                              aTGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGt                              <  2:2608714/91‑1 (MQ=255)

                                                               tGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTa                             >  1:1064443/1‑91 (MQ=255)

                                                               tGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCCATGCAACTGTa                             >  1:3511215/1‑91 (MQ=255)

                                                                ggCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTaa                            >  1:1499215/1‑91 (MQ=255)

                                                                ggCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTaa                            <  2:1876045/91‑1 (MQ=255)

                                                                  cTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAAtt                          >  2:553654/1‑91 (MQ=255)

                                                                   tttCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTa                         >  2:145663/1‑91 (MQ=255)

                                                                    ttCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAg                        <  2:2842405/91‑1 (MQ=255)

                                                                      cATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAgtg                      >  1:842744/1‑91 (MQ=255)

                                                                       aTGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAgtgg                     <  2:549503/91‑1 (MQ=255)

                                                                        tGAACCGCATCGTCAGCATTTAATTGCTGAATTTGATGAACTGAAAAGAATTGCTAGTCAACATAACGCCTATGCAACTGTAATTAgtggc                    >  2:3612822/1‑90 (MQ=255)

                                                                          aaCCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGc                  >  2:2764478/1‑91 (MQ=255)

                                                                           aCCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCt                 >  1:2875786/1‑91 (MQ=255)

                                                                           aCCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCt                 >  1:448002/1‑91 (MQ=255)

                                                                            ccGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTg                >  1:2013341/1‑91 (MQ=255)

                                                                              gTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGa              <  1:1377748/91‑1 (MQ=255)

                                                                               tATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGAc             <  1:2478372/91‑1 (MQ=255)

                                                                               tATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGAc             >  2:2100819/1‑91 (MQ=255)

                                                                                aTCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGAcc            >  2:3558021/1‑91 (MQ=255)

                                                                                 tCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCa           <  1:3800215/91‑1 (MQ=255)

                                                                                 tCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCa           >  2:1651781/1‑91 (MQ=255)

                                                                                 tCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCa           >  2:2178059/1‑91 (MQ=255)

                                                                                 aattCAGCAATTAAATGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCa           >  2:2463207/4‑91 (MQ=255)

                                                                                  cGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCaa          >  1:644913/1‑91 (MQ=255)

                                                                                   ttCAGCAATTAAATGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         >  1:368777/2‑91 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTg                <  1:3141216/84‑1 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         >  1:2156988/1‑91 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         <  2:860573/91‑1 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         <  1:3026224/91‑1 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         >  1:1723635/1‑91 (MQ=255)

                                                                                   gTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAAc         <  1:374211/91‑1 (MQ=255)

                                                                                     cAGCAATTAAATGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTa       >  1:533208/1‑91 (MQ=255)

                                                                                      aGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTAt      <  2:2318105/91‑1 (MQ=255)

                                                                                      aGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTAt      <  2:842744/91‑1 (MQ=255)

                                                                                        cATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTAttt    <  2:3739630/91‑1 (MQ=255)

                                                                                         aTTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTAtttt   <  1:3481928/91‑1 (MQ=255)

TACAATGATTTGTGCATTAGCACAGCACAATTATGAATTAGCAGGTAAACTCATGCAACAAGATGGCTTTCATGAACCGTATCGTCAGCATTTAATTGCTGAATTTGATGAAGTGAAAACAATTGCTAGTCAACATAATGCCTATGCAACTGTAATTAGTGGTGCTGGACCAACTATTTTA  >  USA300TCH1516_ALE/1360294‑1360474

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