BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,138,306	A→T	5.2%	F104Y (TTT→TAT) ‡	USA300HOU_RS10905 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,138,306	0	A	T	5.2%	864.9 / 26.5	250	F104Y (TTT→TAT) ‡	USA300HOU_RS10905	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (123/114); new base T (10/3); total (133/117)
Fisher's exact test for biased strand distribution p-value = 9.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 261 total aligned reads displayed.

CTATTCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGT  >  USA300TCH1516_ALE/2138216‑2138395

  aTTCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATaataaat                                                                                         >  2:1666683/1‑91 (MQ=255)

   ttCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATaataaata                                                                                        <  1:776846/91‑1 (MQ=255)

      ttttGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAAtatta                                                                                     <  1:1502045/91‑1 (MQ=255)

       tttGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAAtattat                                                                                    <  2:1634317/91‑1 (MQ=255)

         tGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataa                                                                                  >  1:2589600/1‑91 (MQ=255)

          gAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataat                                                                                 <  2:3880293/91‑1 (MQ=255)

          gAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTataat                                                                                 <  2:1873337/91‑1 (MQ=255)

            aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATATTTATTataatac                                                                               <  1:3167819/91‑2 (MQ=255)

            aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                               <  1:395198/91‑1 (MQ=255)

            aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                               <  2:192010/91‑1 (MQ=255)

            aaaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAAATAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAg                                                                               <  2:21895/91‑1 (MQ=255)

             aaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATATTTATTataatact                                                                              <  1:1138650/91‑3 (MQ=39)

             aaaTATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGc                                                                              <  1:3472170/91‑1 (MQ=255)

                tatCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTAt                                                                           <  1:1245142/91‑1 (MQ=255)

                tatCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTAt                                                                           >  2:1549380/1‑91 (MQ=255)

                  tCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCg                                                                         <  1:2472208/91‑1 (MQ=255)

                   cTGCATATCTTCACAGTAATTTTGTAACAGATTAATAAAATTATAAAAACTTAACAACAAGTATTATAATAAATATTATAATAGCTATCGt                                                                        <  2:3607926/91‑1 (MQ=255)

                   cTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGt                                                                        >  2:3562681/1‑91 (MQ=255)

                     gCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTAt                                                                      >  1:2869583/1‑91 (MQ=255)

                      cATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTAtt                                                                     <  2:2063782/91‑1 (MQ=255)

                          tCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAAt                                                                 <  1:2908184/91‑1 (MQ=255)

                           cTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATa                                                                <  1:1130506/91‑1 (MQ=255)

                              cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                             <  2:1787352/91‑1 (MQ=255)

                              cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                             >  1:712511/1‑91 (MQ=255)

                              cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                             >  2:2459664/1‑91 (MQ=255)

                              cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                             >  2:2924704/1‑91 (MQ=255)

                              cacaGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGtt                                                             >  1:3416422/1‑91 (MQ=255)

                                caGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTa                                                           >  2:2885013/1‑91 (MQ=255)

                                caGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTa                                                           >  2:954824/1‑91 (MQ=255)

                                 aGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAg                                                          <  1:3393645/91‑1 (MQ=255)

                                    aaTTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAat                                                       <  1:2718824/91‑1 (MQ=255)

                                      ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                     >  2:3502203/1‑91 (MQ=255)

                                      ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                     <  2:125565/91‑1 (MQ=255)

                                      ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                     <  2:125215/91‑1 (MQ=255)

                                      ttttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAatat                                                     <  2:3726880/91‑1 (MQ=255)

                                        ttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGc                                                   <  1:2298644/91‑1 (MQ=255)

                                        ttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGc                                                   >  2:1468460/1‑91 (MQ=255)

                                        ttATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGc                                                   <  1:3797593/91‑1 (MQ=255)

                                          ataACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCaa                                                 >  2:59604/1‑91 (MQ=255)

                                          ataACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCaa                                                 >  2:2316343/1‑91 (MQ=255)

                                           taACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCaaa                                                <  2:454277/91‑1 (MQ=255)

                                            aaCAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAAc                                               >  1:1724426/1‑91 (MQ=255)

                                             aCAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACt                                              <  1:1291432/91‑1 (MQ=255)

                                              cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACt                                              >  2:599988/1‑90 (MQ=255)

                                              cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACt                                              <  1:599332/90‑1 (MQ=255)

                                              cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTg                                             <  1:289369/91‑1 (MQ=255)

                                              cAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTg                                             >  1:2049843/1‑91 (MQ=255)

                                               aGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGa                                            <  2:2511378/91‑1 (MQ=255)

                                               aGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGa                                            <  1:1623597/91‑1 (MQ=255)

                                                  ttAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaa                                         >  1:2138266/1‑91 (MQ=255)

                                                    aataaaaataaaAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaat                                       <  2:1531543/91‑1 (MQ=255)

                                                    aataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaat                                       >  1:547297/1‑91 (MQ=255)

                                                     ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                      >  2:3409299/1‑91 (MQ=255)

                                                     ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCAATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                      >  2:2872844/1‑91 (MQ=255)

                                                     ataaAATTAAAAAAACATAACAACAAGTATTATAATAAATACTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaata                                      >  2:2952364/1‑91 (MQ=255)

                                                       aaaaTTAGAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaa                                    >  2:2730025/1‑91 (MQ=255)

                                                        aaaTTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGaaaaataaaa                                   >  1:99559/1‑91 (MQ=255)

                                                            tAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATca                               <  2:1795413/91‑1 (MQ=255)

                                                             aaaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcac                              >  2:3785868/1‑91 (MQ=255)

                                                             aaaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcac                              >  2:1119445/1‑91 (MQ=255)

                                                              aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcaca                             >  1:115771/1‑91 (MQ=255)

                                                              aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcaca                             >  2:873398/1‑91 (MQ=255)

                                                              aaaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATcaca                             >  2:113776/1‑91 (MQ=255)

                                                               aaaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACaa                            <  2:3696649/91‑1 (MQ=255)

                                                                aaaaCATAACAACAAGTATTATAATAAATCTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                           >  2:2760278/1‑91 (MQ=255)

                                                                aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                           >  2:175681/1‑91 (MQ=255)

                                                                aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                           <  2:3252659/91‑1 (MQ=255)

                                                                aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                           >  2:722912/1‑91 (MQ=255)

                                                                aaaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAt                           >  1:1895507/1‑91 (MQ=255)

                                                                 aaaCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAAtt                          <  2:3042359/91‑1 (MQ=255)

                                                                   aCATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTaa                        >  1:625945/1‑91 (MQ=255)

                                                                    cATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTaaa                       >  1:492813/1‑91 (MQ=255)

                                                                     aTAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTaaaa                      <  2:2425122/91‑1 (MQ=255)

                                                                     aTAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTaaaa                      <  1:1618588/91‑1 (MQ=255)

                                                                      tAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                     >  2:3459276/1‑91 (MQ=255)

                                                                      tAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                     <  1:2144979/91‑1 (MQ=255)

                                                                       aacaacAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATa                    >  1:1668619/1‑91 (MQ=255)

                                                                       aacaacAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATa                    <  1:1902243/91‑1 (MQ=255)

                                                                        acaacaAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATAc                   <  2:2971165/91‑1 (MQ=255)

                                                                           acaaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTaa                <  1:910975/91‑1 (MQ=255)

                                                                             aaGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGt              >  1:100661/1‑91 (MQ=255)

                                                                               gTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTaa            <  1:415814/91‑1 (MQ=255)

                                                                               gTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTaa            <  2:522670/91‑1 (MQ=255)

                                                                                tattatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAAt           >  1:3288203/1‑91 (MQ=255)

                                                                                 attatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg          >  2:3809258/1‑91 (MQ=255)

                                                                                 attatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg          >  2:2647639/1‑91 (MQ=255)

                                                                                 attatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg          <  2:629789/91‑1 (MQ=255)

                                                                                 attatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATg          <  1:2796948/91‑1 (MQ=255)

                                                                                  ttatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt         >  2:2792576/1‑91 (MQ=255)

                                                                                  ttatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt         >  1:3309935/1‑91 (MQ=255)

                                                                                  ttatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt         >  1:824187/1‑91 (MQ=255)

                                                                                  ttatAATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt         >  1:1758091/1‑91 (MQ=255)

                                                                                  ttatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAAt                     >  1:445627/1‑79 (MQ=255)

                                                                                  ttatAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGt         <  1:3705154/91‑1 (MQ=255)

                                                                                   tataATATTTATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTa        >  2:1203695/1‑91 (MQ=255)

                                                                                   tataATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATACAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTa        >  2:3477777/1‑91 (MQ=255)

                                                                                    ataataAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTaa       >  2:459714/1‑91 (MQ=255)

                                                                                     taataaATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTaaa      <  2:684482/91‑1 (MQ=255)

                                                                                      aataaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAAt     <  2:791215/91‑1 (MQ=255)

                                                                                       ataaataTTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATg    <  1:3502203/91‑1 (MQ=255)

CTATTCTTTTGAAAAATATCTGCATATCTTCACAGTAATTTTATAACAGATTAATAAAATTAAAAAAACATAACAACAAGTATTATAATAAATATTATAATAGCTATCGTATTAAATAGTTTAGAATATGCAAACTGAAAAATAAAAATCACAATTAAAATACTAAGTAATGTAAATGGT  >  USA300TCH1516_ALE/2138216‑2138395

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