BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,010,363	C→T	100%	C146Y (TGT→TAT)	USA300HOU_RS10125 ←	Putative multidrug export ATP‑binding/permease protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,010,363	0	C	T	98.2%	391.5 / ‑2.2	110	C146Y (TGT→TAT)	USA300HOU_RS10125	Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (43/64); total (46/64)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 111 total aligned reads displayed.

ATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCT  >  USA300TCH1516_ALE/2010276‑2010449

atgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAt                                                                                     >  1:372523/1‑91 (MQ=255)

 tgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATc                                                                                    >  1:3104467/1‑91 (MQ=255)

  gtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATCTCACATCTAAAAATAACATTATGGATAGTGCAATAATAATTGTTATAGAATCt                                                                                   >  1:3860209/1‑91 (MQ=255)

     aaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAc                                                                                <  1:100309/91‑1 (MQ=255)

     aaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAc                                                                                <  1:264717/91‑1 (MQ=255)

     aaaggGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATTTAAc                                                                                <  1:3267913/91‑1 (MQ=255)

      aatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAcc                                                                               >  1:2194431/1‑91 (MQ=255)

       atgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCa                                                                              >  1:1220305/1‑91 (MQ=255)

       atgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCa                                                                              <  1:3001785/91‑1 (MQ=255)

        tgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCaa                                                                             >  1:2672988/1‑91 (MQ=255)

         gggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCaaa                                                                            <  1:885741/91‑1 (MQ=255)

          ggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAat                                                                           >  1:3860003/1‑91 (MQ=255)

          ggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAat                                                                           <  1:3605693/91‑1 (MQ=255)

            aaGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAatat                                                                         >  1:3253202/1‑91 (MQ=255)

            aaGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAatat                                                                         >  1:1520848/1‑91 (MQ=255)

              gATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTc                                                                       <  1:1425568/91‑1 (MQ=255)

               aTAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                      <  1:140378/91‑1 (MQ=255)

                tAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCAt                                                                     >  1:2936494/1‑91 (MQ=255)

                    cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                 <  1:1551791/91‑1 (MQ=255)

                    cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                 >  1:2302433/1‑91 (MQ=255)

                    cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                 >  1:2301910/1‑91 (MQ=255)

                    cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                 <  1:190957/91‑1 (MQ=255)

                     aGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAcc                                                                <  1:776205/91‑1 (MQ=255)

                     aGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAcc                                                                >  1:1896767/1‑91 (MQ=255)

                       tgctgcTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCg                                                              <  1:2343186/91‑1 (MQ=255)

                         ctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGt                                                            <  1:674286/91‑1 (MQ=255)

                          tgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGtt                                                           <  1:3353643/91‑1 (MQ=255)

                          tgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGtt                                                           <  1:418735/91‑1 (MQ=255)

                           gctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTa                                                          >  1:1037348/1‑91 (MQ=255)

                                  tCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaa                                                   >  1:2345288/1‑91 (MQ=255)

                                  tCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaa                                                   >  1:2345267/1‑91 (MQ=255)

                                    aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAt                                                 >  1:1202423/1‑91 (MQ=255)

                                    aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAt                                                 <  1:250228/91‑1 (MQ=255)

                                     aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATc                                                >  1:2376198/1‑91 (MQ=255)

                                     aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATc                                                >  1:2375912/1‑91 (MQ=255)

                                     aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATc                                                >  1:1149601/1‑91 (MQ=255)

                                         cacaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                            >  1:1329665/1‑91 (MQ=255)

                                         cacaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                            >  1:1136777/1‑91 (MQ=255)

                                         cacaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                            <  1:2987067/91‑1 (MQ=255)

                                          acaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                           <  1:3374810/91‑1 (MQ=255)

                                          acaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                           <  1:2129821/91‑1 (MQ=255)

                                          acaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                           >  1:3809282/1‑91 (MQ=255)

                                             tCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                        >  1:1903861/1‑91 (MQ=255)

                                             tCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                        >  1:3409661/1‑91 (MQ=255)

                                             tCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                        >  1:3409713/1‑91 (MQ=255)

                                                aaaaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTtgtttgtt                                     >  1:2174441/1‑91 (MQ=255)

                                                aaaaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTtgtttgtt                                     <  1:151377/91‑1 (MQ=255)

                                                   aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCaa                                  <  1:3216860/91‑1 (MQ=255)

                                                   aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCaa                                  <  1:685041/91‑1 (MQ=255)

                                                   aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCaa                                  >  1:1284865/1‑91 (MQ=255)

                                                   aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCaa                                  <  1:681546/91‑1 (MQ=255)

                                                    aGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAc                                 <  1:35916/91‑1 (MQ=255)

                                                    aGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAc                                 <  1:287010/91‑1 (MQ=255)

                                                    aGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAc                                 <  1:1486976/91‑1 (MQ=255)

                                                    aGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAc                                 <  1:3213855/91‑1 (MQ=255)

                                                       aCATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatc                              >  1:1892401/1‑91 (MQ=255)

                                                          ttATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCAtt                           >  1:2428519/1‑91 (MQ=255)

                                                          ttATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCAtt                           >  1:2428403/1‑91 (MQ=255)

                                                            aTGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTaa                         <  1:3997089/91‑1 (MQ=255)

                                                            aTGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTaa                         <  1:749771/91‑1 (MQ=255)

                                                            aTGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTaa                         >  1:2657861/1‑91 (MQ=255)

                                                             tGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAAt                        <  1:2353959/91‑1 (MQ=255)

                                                               gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCa                      <  1:1508676/91‑1 (MQ=255)

                                                                  aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                   <  1:2646101/91‑1 (MQ=255)

                                                                  aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                   <  1:328868/91‑1 (MQ=255)

                                                                  aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                   <  1:329069/91‑1 (MQ=255)

                                                                  aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                   <  1:329082/91‑1 (MQ=255)

                                                                  aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                   <  1:3896877/91‑1 (MQ=255)

                                                                   gTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                  <  1:3582662/91‑1 (MQ=255)

                                                                    tGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTa                 >  1:2959400/1‑91 (MQ=255)

                                                                     gCAATAATAATTGTTATATAATCTAACCAAATATTCATTATCCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAg                <  1:2046230/91‑1 (MQ=255)

                                                                     gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAg                <  1:3042106/91‑1 (MQ=255)

                                                                     gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAg                <  1:178143/91‑1 (MQ=255)

                                                                     gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAg                <  1:179802/91‑1 (MQ=255)

                                                                     gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAg                <  1:620785/91‑1 (MQ=255)

                                                                      cAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGTAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa               <  1:3494073/91‑1 (MQ=255)

                                                                      cAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa               <  1:2693342/91‑1 (MQ=255)

                                                                      cAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa               <  1:2693243/91‑1 (MQ=255)

                                                                      cAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa               <  1:1339662/91‑1 (MQ=255)

                                                                        ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata             <  1:4026948/91‑1 (MQ=255)

                                                                        ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata             >  1:1156457/1‑91 (MQ=255)

                                                                        ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata             <  1:1789372/91‑1 (MQ=255)

                                                                         taataatTGTTGTATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat            <  1:3763165/91‑1 (MQ=255)

                                                                         taataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat            >  1:2756898/1‑91 (MQ=255)

                                                                          aataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATAtt           >  1:3638219/1‑91 (MQ=255)

                                                                            taatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTAc         <  1:3098827/91‑1 (MQ=255)

                                                                            taatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTAc         <  1:2814267/91‑1 (MQ=255)

                                                                             aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACt        >  1:2114702/1‑91 (MQ=255)

                                                                             aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAGCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACt        <  1:3600930/91‑1 (MQ=255)

                                                                              atTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACtt       <  1:3724650/91‑1 (MQ=255)

                                                                               ttGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAGTGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTg      >  1:3594347/1‑91 (MQ=255)

                                                                                tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa     <  1:3820722/91‑1 (MQ=255)

                                                                                tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa     >  1:459097/1‑91 (MQ=255)

                                                                                tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa     <  1:1616874/91‑1 (MQ=255)

                                                                                tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa     <  1:1616848/91‑1 (MQ=255)

                                                                                 gTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAc    <  1:2209716/91‑1 (MQ=255)

                                                                                  ttATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc   <  1:1103552/91‑1 (MQ=255)

                                                                                  ttATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc   <  1:1102874/91‑1 (MQ=255)

                                                                                  ttATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc   <  1:2123435/91‑1 (MQ=255)

                                                                                   tatataATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCt  >  1:4073710/1‑91 (MQ=255)

ATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCT  >  USA300TCH1516_ALE/2010276‑2010449

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