BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	667,064	T→G	8.1%	D18E (GAT→GAG)	USA300TCH1516_00600 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	667,064	0	T	G	8.1%	408.6 / 22.5	123	D18E (GAT→GAG)	USA300TCH1516_00600	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (41/72); new base G (3/7); total (44/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 127 total aligned reads displayed.

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATT  >  USA300TCH1516_ALE/666974‑667151

gCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGAt                                                                                         >  1:2516069/1‑91 (MQ=255)

gCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGAt                                                                                         >  1:2516302/1‑91 (MQ=255)

gCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGAt                                                                                         >  1:1240009/1‑91 (MQ=255)

gCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGAt                                                                                         <  1:1882738/91‑1 (MQ=255)

 cTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATg                                                                                        <  1:97768/91‑1 (MQ=255)

 cTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATg                                                                                        <  1:145472/91‑1 (MQ=255)

 cTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATg                                                                                        <  1:355935/91‑1 (MQ=255)

 cTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATg                                                                                        <  1:2810174/91‑1 (MQ=255)

  tAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATgg                                                                                       >  1:803808/1‑91 (MQ=255)

  tAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATgg                                                                                       >  1:842934/1‑91 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGt                                                                                      <  1:1782688/91‑1 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGt                                                                                      <  1:1367860/91‑1 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGt                                                                                      <  1:1481613/91‑1 (MQ=255)

    aaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTc                                                                                     <  1:303296/91‑1 (MQ=255)

    aaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTc                                                                                     <  1:303578/91‑1 (MQ=255)

    aaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTc                                                                                     <  1:3473047/91‑1 (MQ=255)

     aaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCa                                                                                    <  1:1671455/91‑1 (MQ=255)

       ttGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATa                                                                                  <  1:2083775/91‑1 (MQ=255)

         gTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATaaa                                                                                >  1:2359098/1‑91 (MQ=255)

         gTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATaaa                                                                                <  1:668340/91‑1 (MQ=255)

           tATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAgg                                                                              >  1:1160620/1‑91 (MQ=255)

           tATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAgg                                                                              <  1:3723909/91‑1 (MQ=255)

             tCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTa                                                                            <  1:311318/91‑1 (MQ=255)

              cGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTAc                                                                           <  1:3822231/91‑1 (MQ=255)

              cGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTAc                                                                           <  1:2912291/91‑1 (MQ=255)

                ttAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACaa                                                                         <  1:326182/91‑1 (MQ=255)

                ttAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACaa                                                                         <  1:783707/91‑1 (MQ=255)

                ttAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACaa                                                                         <  1:783423/91‑1 (MQ=255)

                 tAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACaaa                                                                        <  1:338961/91‑1 (MQ=255)

                  aaaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAat                                                                       <  1:1139527/91‑1 (MQ=255)

                   aaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAata                                                                      <  1:3452037/91‑1 (MQ=255)

                   aaTATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAata                                                                      >  1:2842711/1‑91 (MQ=255)

                    atatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAatat                                                                     >  1:3015677/1‑91 (MQ=255)

                    atatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAatat                                                                     <  1:2312589/91‑1 (MQ=255)

                    atatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAatat                                                                     <  1:247540/91‑1 (MQ=255)

                     tatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATAtt                                                                    <  1:2872101/91‑1 (MQ=255)

                     tatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATAtt                                                                    <  1:87113/91‑1 (MQ=255)

                     tatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATAtt                                                                    <  1:2872218/91‑1 (MQ=255)

                     tatGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATAtt                                                                    <  1:1217436/91‑1 (MQ=255)

                      atGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTg                                                                   >  1:86619/1‑91 (MQ=255)

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCGTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTgg                                                                  <  1:3528360/91‑1 (MQ=255)

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCCTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTgg                                                                  <  1:3528256/91‑1 (MQ=255)

                        ggAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGt                                                                 >  1:538915/1‑91 (MQ=255)

                        ggAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGt                                                                 >  1:4057375/1‑91 (MQ=255)

                         gAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGtt                                                                >  1:534656/1‑91 (MQ=255)

                           ggTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCa                                                              <  1:2833184/91‑1 (MQ=255)

                            gTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCaa                                                             <  1:2963694/91‑1 (MQ=255)

                             ttATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAAt                                                            >  1:831413/1‑91 (MQ=255)

                              tataTCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAAtt                                                           <  1:461751/91‑1 (MQ=255)

                                 atCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGc                                                        <  1:969060/91‑1 (MQ=255)

                                   cATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCtt                                                      <  1:2321287/91‑1 (MQ=255)

                                      tGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTGGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGa                                                   <  1:3591818/91‑1 (MQ=255)

                                      tGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTGGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGa                                                   <  1:3591773/91‑1 (MQ=255)

                                        tCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACa                                                 <  1:1674804/91‑1 (MQ=255)

                                        tCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACa                                                 <  1:2124173/91‑1 (MQ=255)

                                          tGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTGATAAAGGTACAAATATTTGTTCAATTAGCGTAGACAat                                               >  1:520256/1‑91 (MQ=255)

                                            agaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAatat                                             >  1:404889/1‑91 (MQ=255)

                                              aaaaaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTg                                           <  1:971208/91‑1 (MQ=255)

                                               aaaaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGc                                          <  1:1715351/91‑1 (MQ=255)

                                                aaaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCa                                         <  1:1938134/91‑1 (MQ=255)

                                                aaaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCa                                         >  1:2404528/1‑91 (MQ=255)

                                                 aaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAc                                        <  1:3617983/91‑1 (MQ=255)

                                                   aGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATCAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCa                                      >  1:3292533/1‑91 (MQ=255)

                                                   aGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCa                                      <  1:1655659/91‑1 (MQ=255)

                                                   aGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCa                                      >  1:4050315/1‑91 (MQ=255)

                                                   aGGCTTTAATTTCAATATCATAAAAAGTGATCCCATAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCa                                      <  1:1540565/91‑1 (MQ=255)

                                                         taattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGttttt                                >  1:3557190/1‑91 (MQ=255)

                                                         taattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGttttt                                >  1:1331291/1‑91 (MQ=255)

                                                         taattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGttttt                                <  1:1752239/91‑1 (MQ=255)

                                                           attttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTAt                              <  1:3050889/91‑1 (MQ=255)

                                                              ttaatATCAGAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGa                           >  1:2086172/1‑91 (MQ=255)

                                                                aatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATg                         >  1:1855526/1‑91 (MQ=255)

                                                                  tatCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAACATTGCACCAGTTTTTATCGACGtt                       >  1:2431841/1‑91 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                   >  1:2944863/1‑91 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                   >  1:2944827/1‑91 (MQ=255)

                                                                      aTAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTa                   <  1:4110591/91‑1 (MQ=255)

                                                                        aaaaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAAc                 >  1:1361516/1‑91 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAActta              <  1:1251069/91‑4 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAActta              <  1:1250962/91‑4 (MQ=255)

                                                                             aTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaaga            <  1:1084184/91‑1 (MQ=255)

                                                                             aTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaaga            <  1:2637704/91‑1 (MQ=255)

                                                                             aTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaaga            <  1:1086813/91‑1 (MQ=255)

                                                                             aTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaaga            >  1:1494701/1‑91 (MQ=255)

                                                                              tGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaa           <  1:3783378/91‑1 (MQ=255)

                                                                              tGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaa           <  1:1210427/91‑1 (MQ=255)

                                                                               gACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaag          <  1:2012251/91‑1 (MQ=255)

                                                                               gACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaag          >  1:3585493/1‑91 (MQ=255)

                                                                                aCCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGc         >  1:592619/1‑91 (MQ=255)

                                                                                aCCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGc         >  1:1221002/1‑91 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa        >  1:1027574/1‑91 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa        <  1:209125/91‑1 (MQ=255)

                                                                                 cccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTGAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa        <  1:3291032/91‑1 (MQ=255)

                                                                                  catctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt       >  1:2178405/3‑91 (MQ=255)

                                                                                   atctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt      >  1:1078546/2‑91 (MQ=255)

                                                                                    tctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAttt     <  1:2844058/91‑1 (MQ=255)

                                                                                     ctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTa    <  1:666793/91‑1 (MQ=255)

                                                                                     ctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTa    >  1:1473039/1‑91 (MQ=255)

                                                                                      tAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAt   >  1:359799/1‑91 (MQ=255)

                                                                                      tAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAt   <  1:2496346/91‑1 (MQ=255)

                                                                                       aGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTAtt  <  1:822979/91‑1 (MQ=255)

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATT  >  USA300TCH1516_ALE/666974‑667151

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