BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	667,064	T→G	11.9%	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	11.9%	649.3 / 71.9	211	D18E (GAT→GAG)	USA300TCH1516_00600	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (89/96); new base G (13/12); total (103/108)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 214 total aligned reads displayed.

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTGAT  >  USA300TCH1516_ALE/666974‑667154

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

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

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

  tAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATgg                                                                                          <  2:3108578/91‑1 (MQ=255)

   aaaaTTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGt                                                                                         >  1:2244980/1‑91 (MQ=255)

   aaaaTTGTTATAGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGt                                                                                         <  2:1807669/91‑1 (MQ=255)

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

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

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

        tGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTCAttt                                                                                    <  1:1905682/91‑3 (MQ=255)

          ttATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAg                                                                                  >  2:2763694/1‑91 (MQ=255)

          ttATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTCAttttt                                                                                  <  2:3576256/91‑5 (MQ=39)

          ttATCGTTAAATATGGAGGTTATATCATTGACTGAAAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCATCTAGAGGGTCAttttt                                                                                  <  2:1621915/91‑5 (MQ=38)

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

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

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

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

             tCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTa                                                                               >  2:3229771/1‑91 (MQ=255)

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

               gTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACa                                                                             <  1:2687418/91‑1 (MQ=255)

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

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

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

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTgg                                                                     >  1:1242625/1‑91 (MQ=255)

                       tGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTgg                                                                     >  2:946455/1‑91 (MQ=255)

                        ggAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGt                                                                    <  1:1440416/91‑1 (MQ=255)

                        ggAGGTTATATCAATGTATGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACTCTCTAGATGGTCATAAAGGTACAAATATTGGt                                                                    <  2:2242867/91‑1 (MQ=255)

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

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

                          aGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTc                                                                  <  1:1188878/91‑1 (MQ=255)

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

                             ttATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAAt                                                               <  1:2032267/91‑1 (MQ=255)

                             ttATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAAt                                                               <  1:1827440/91‑1 (MQ=255)

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

                               atatCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTa                                                             >  1:781517/1‑91 (MQ=255)

                                tatCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAg                                                            <  2:3068106/91‑1 (MQ=255)

                                  tCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCt                                                          <  2:1394078/91‑1 (MQ=255)

                                   cATTGTCTGAGAAAACAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCtt                                                         >  2:3168723/1‑91 (MQ=255)

                                   cATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCtt                                                         >  2:3169020/1‑91 (MQ=255)

                                    aTTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTa                                                        <  2:2514144/91‑1 (MQ=255)

                                     ttGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       >  2:3752075/1‑91 (MQ=255)

                                     ttGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATCGTCATAAAGGTACAAATATTGGTTCAATTAGCTTaa                                                       >  2:3839481/1‑90 (MQ=255)

                                     atGTCAGAAAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAg                                                       <  2:1456031/90‑1 (MQ=255)

                                      tGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGa                                                      <  2:2895759/91‑1 (MQ=255)

                                       gTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGAc                                                     >  1:2965329/1‑91 (MQ=255)

                                          tGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAat                                                  >  1:1079227/1‑91 (MQ=255)

                                            agaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAatat                                                <  1:1340685/91‑1 (MQ=255)

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

                                             gaAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATAtt                                               <  1:535723/91‑1 (MQ=255)

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

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

                                               aaaaaGGCTTTAATTTTAATATAATAAAAAATGACCCTCTTGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGc                                             <  2:127104/91‑1 (MQ=255)

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

                                                 aaaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAc                                           >  2:1174917/1‑91 (MQ=255)

                                                 aaaGGCTTTAATTTTAATATCATAAAAAATGAACCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAc                                           <  1:3293786/91‑1 (MQ=255)

                                                  aaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAcc                                          <  2:661983/91‑1 (MQ=255)

                                                  aaGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCAcc                                          <  1:2216063/91‑1 (MQ=255)

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

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

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

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        >  1:2641280/1‑91 (MQ=255)

                                                    ggCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAg                                        >  2:538091/1‑91 (MQ=255)

                                                     gCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGt                                       <  2:1036015/91‑1 (MQ=255)

                                                     gCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGt                                       >  2:213760/1‑91 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    <  2:91556/91‑1 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    <  1:3056542/91‑1 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    <  2:3680171/91‑1 (MQ=255)

                                                        ttaattttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGtttt                                    >  2:1362178/1‑91 (MQ=255)

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

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

                                                              ttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGa                              <  1:1931694/91‑1 (MQ=255)

                                                              ttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGa                              >  1:2533795/1‑91 (MQ=255)

                                                              ttaatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGa                              <  2:3290240/91‑1 (MQ=255)

                                                               taatATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGAt                             <  2:3069899/91‑1 (MQ=255)

                                                                 atatCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGt                           >  2:2784623/1‑91 (MQ=255)

                                                                   atCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTg                         >  1:1600218/1‑91 (MQ=255)

                                                                     cATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCt                       <  1:2478986/91‑1 (MQ=255)

                                                                       tAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTaa                     <  1:3391722/91‑1 (MQ=255)

                                                                         aaaaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACa                   >  1:1480074/1‑91 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACaaa                 >  2:2553688/1‑91 (MQ=255)

                                                                           aaaTGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACaaa                 >  1:2897196/1‑91 (MQ=255)

                                                                              tGACCATCTAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAaagaa              >  1:2719395/1‑91 (MQ=255)

                                                                               gACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGc                        <  1:2374700/80‑1 (MQ=255)

                                                                               gACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGc                        >  2:2374700/1‑80 (MQ=255)

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

                                                                                aCCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGc            <  2:1920225/91‑1 (MQ=255)

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

                                                                                 ccatctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCa           >  2:171910/4‑91 (MQ=255)

                                                                                  ccTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt          >  1:2598579/1‑91 (MQ=255)

                                                                                  catctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt          >  2:2939481/3‑91 (MQ=255)

                                                                                  catctAGAGGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAt          >  2:3525594/3‑91 (MQ=255)

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

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

                                                                                   ctctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:502018/1‑91 (MQ=255)

                                                                                   ctctAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCATAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCAtt         >  1:3357669/1‑91 (MQ=255)

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

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

                                                                                      tAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTTTTTTTGTTTCTAACAAAGAAGCATTTAt      >  2:2583239/1‑91 (MQ=255)

                                                                                        gATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTg    <  1:443964/91‑1 (MQ=255)

                                                                                          tGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTGat  <  2:1663311/91‑1 (MQ=255)

GCTAAAATTGTTATCGTTAAATATGGAGGTTATATCATTGTCTGAGAAAAAAGGCTTTAATTTTAATATCATAAAAAATGACCCTCTAGATGGTCATAAAGGTACAAATATTGGTTCAATTAGCTTAGACAATATTGCACCAGTTTTTATCGATGTTGCTAACAAAGAAGCATTTATTGAT  >  USA300TCH1516_ALE/666974‑667154

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