BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,002,366	G→T	A935S (GCT→TCT)	flu →	antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,002,366	0	G	T	100.0%	304.2 / NA	81	A935S (GCT→TCT)	flu	antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (25/56); total (25/56)

Only 100 of 121 total aligned reads displayed.

GGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTCCCGG  >  REL606/3002293‑3002438

ggAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGctct                                                                          >  2:3253749/1‑76 (MQ=255)

 gAAGCCCGTGTCCGGGAAAATATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTg                                                                         <  2:2391328/76‑1 (MQ=255)

 gAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTg                                                                         <  2:4497587/76‑1 (MQ=255)

 gAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTg                                                                         <  1:2503695/76‑1 (MQ=255)

    gCCCGTGTCCGGGAAAATATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAg                                                                      <  2:884691/76‑1 (MQ=255)

    gCCCGTGTCCGGGAAAATATCCCCCTGGGCGTGCAGGCCGGTTATGCCCGCAGCGTCAGCGGCAGCAGCTCTGAAg                                                                      <  2:3796690/76‑1 (MQ=255)

    gCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAg                                                                      <  1:4064735/76‑1 (MQ=255)

      ccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCTCTGGAGGc                                                                    >  1:3400540/1‑76 (MQ=255)

       cGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGCGCTCTGAAAGCt                                                                   >  2:3968687/1‑76 (MQ=255)

        gtgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCta                                                                  <  1:2431711/76‑1 (MQ=255)

        gtgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTCAAGGCta                                                                  >  1:4140206/1‑76 (MQ=255)

         tgtCGGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCtat                                                                 <  2:2902992/76‑1 (MQ=255)

         tgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGCAGGCtat                                                                 >  1:4023027/1‑76 (MQ=255)

         tgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAAGCtat                                                                 >  1:3564051/1‑76 (MQ=255)

          gtCCGGGAAAATATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCtat‑a                                                               <  1:902830/76‑1 (MQ=255)

          gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTGTGGACGCtat‑a                                                               >  1:4756792/1‑76 (MQ=255)

          gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCtat‑a                                                               <  2:3241002/76‑1 (MQ=255)

           tCCGGGAAAATATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑aa                                                              <  1:1351906/76‑1 (MQ=255)

           tCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCTCTGAAGGCTAT‑aa                                                              >  1:4606204/1‑76 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑aaa                                                             >  1:2745117/1‑75 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             <  2:1301392/76‑1 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             <  1:3962892/76‑1 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             >  1:3101750/1‑76 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             <  2:740631/76‑1 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             <  2:2516236/76‑1 (MQ=255)

            ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAc                                                             <  1:5011065/76‑1 (MQ=255)

               ggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AAcagg                                                          >  2:1667489/1‑73 (MQ=255)

                gAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTc                                                         <  1:3174989/76‑1 (MQ=255)

                 aaaaTATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCa                                                        <  1:2572529/76‑1 (MQ=255)

                 aaaaTATCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCa                                                        <  1:791573/76‑1 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCCGCTCTGAAGGGTAT‑AACGGGCA‑g                                                      >  1:1638632/1‑76 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑g                                                      <  2:4197361/76‑1 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑g                                                      >  2:3678754/1‑76 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑g                                                      >  1:2937028/1‑76 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑g                                                      >  2:4925750/1‑76 (MQ=255)

                  aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑g                                                      >  1:1477376/1‑76 (MQ=255)

                   aaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCCCTGAAGGCTAT‑AACGGGCA‑gg                                                     >  1:2666652/1‑76 (MQ=255)

                   aaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑gg                                                     >  1:882055/1‑76 (MQ=255)

                   aaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCCGGAGGCTAT‑AACGGTCA‑gg                                                     >  1:4574082/1‑76 (MQ=255)

                     tatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCTCTGAAGGCTTTAAACGTTCA‑GGc                                                    >  2:4274197/1‑76 (MQ=255)

                     tatCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGcc                                                   >  2:4852037/1‑76 (MQ=255)

                        cACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGccccc                                                >  1:2085507/1‑73 (MQ=255)

                        cACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGGCA‑GGccccc                                                >  2:193290/1‑73 (MQ=255)

                        cACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGCTCA‑AACGGTCC‑GGGCaac                                                >  1:4057428/1‑74 (MQ=255)

                         gcccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGc                                               <  1:1145581/75‑1 (MQ=255)

                          cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCCCTGTAGGGTAT‑AAAGGTCA‑GGCCACGCt                                              >  2:933342/1‑76 (MQ=255)

                          cccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCCGAAGGCTAT‑AACCGTCA‑GGCCACGcc                                              >  1:3242781/1‑75 (MQ=255)

                           ccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCCCTGGAGGCTAT‑TACCGGCC‑AGCCACGccg                                             >  2:4164784/1‑74 (MQ=255)

                           ccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTg                                             <  1:1147799/76‑1 (MQ=255)

                           ccTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTg                                             >  1:2475049/1‑76 (MQ=255)

                              gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAAt                                          <  1:772389/76‑1 (MQ=255)

                              gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGGCA‑GGCCACGCCGaaa                                          >  1:4221909/1‑75 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGCGCTCTGAAGGCCAT‑AACGGTCC‑AGGCACGCGGAAagtc                                       >  2:1267041/1‑75 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCCCGCTGAAtgtg                                       >  1:2390995/1‑76 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAAtgtg                                       <  2:1184149/76‑1 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAAtgtg                                       <  2:4310648/76‑1 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAAtgtg                                       <  2:4441317/76‑1 (MQ=255)

                                 cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAAtgtg                                       <  2:2805764/76‑1 (MQ=255)

                                   ttCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGAc                                     >  1:4683629/1‑76 (MQ=255)

                                   ttCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCCGAAtgtggc                                     >  1:443348/1‑74 (MQ=255)

                                    tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGGAGCTCTTAAGGGTAT‑AACCGGCA‑GGGCCCGCTGAAtgtggac                                    >  2:66251/1‑73 (MQ=255)

                                    tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGAcc                                    >  1:710800/1‑76 (MQ=255)

                                     cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCCGCTCCGAAGGCTAT‑AAAGGTCC‑GGCCCCACCCAATGTGAGCt                                   >  2:2931028/1‑76 (MQ=255)

                                     cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCt                                   <  1:3366712/76‑1 (MQ=255)

                                     cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCCGAAGGGTTT‑AACGGTCA‑AGCCACCCTGAATATGAACt                                   >  2:3386831/1‑76 (MQ=255)

                                      aGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCtt                                  <  2:55482/76‑1 (MQ=255)

                                          cGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGa                              <  2:1146977/76‑1 (MQ=255)

                                          cGGTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGa                              >  1:4183408/1‑76 (MQ=255)

                                           ggTTATGCCCACAGCGTCAGCGGCAGGAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGAc                             >  1:4074/1‑76 (MQ=255)

                                           ggTTATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGGACTTCTGAc                             >  1:4127830/1‑76 (MQ=255)

                                             ttATGCCCACAGCGTCAGCGGCAGGAGCTCTGAAGGCTAT‑AAAGGGCA‑GGCCCCGCTGAATGTGGACTTCTTACaa                           >  2:2867836/1‑75 (MQ=255)

                                             ttATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGGGAACTTCTGACaa                           >  1:2704576/1‑75 (MQ=255)

                                             ttATGCCCACAGCGTCAGCGGCAGCAGCTCCGGAGGCTAT‑AAAGGTCAGGGCCAC‑CTGAATTTTACCTTCTGACaa                           >  2:420022/1‑75 (MQ=255)

                                              tATGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGa                          <  2:545100/76‑1 (MQ=255)

                                               aTGCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCCGAATGTGACCTTCTGACAGaa                         >  1:2825841/1‑76 (MQ=255)

                                               aTGCCCACAGCGTCAGCGGCAGCAGCCCTGAAGGCTAT‑AACGGGCA‑GGCCACGCTGAATGGGACCTTCTGACAGaa                         >  2:629037/1‑76 (MQ=255)

                                                 gCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAAcc                       <  2:3501909/76‑1 (MQ=255)

                                                 gCCCACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCCGAATGTGACCTTCTGGCAAAAcc                       >  2:4534333/1‑76 (MQ=255)

                                                  cccACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCa                      >  1:367876/1‑76 (MQ=255)

                                                   ccACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCAt                     <  2:958371/76‑1 (MQ=255)

                                                   ccACAGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCAt                     <  2:1175977/76‑1 (MQ=255)

                                                     acaGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCg                   <  2:2659827/76‑1 (MQ=255)

                                                     acaGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCg                   <  1:121432/76‑1 (MQ=255)

                                                     acaGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCg                   <  1:2700107/76‑1 (MQ=255)

                                                     acaGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCg                   <  2:3674294/76‑1 (MQ=255)

                                                      caGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGc                  <  1:388856/76‑1 (MQ=255)

                                                      caGCGTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGc                  <  2:4186821/76‑1 (MQ=255)

                                                          gTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCctct              <  2:625416/76‑1 (MQ=255)

                                                          gTCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCctct              >  1:1914157/1‑76 (MQ=255)

                                                           tCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGGCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCtt             >  2:328692/1‑76 (MQ=255)

                                                           tCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCtt             >  2:3851890/1‑76 (MQ=255)

                                                           tCAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑AGCCACCCCGAAAGTGACCTTCTGACAGAACCCACGCCCCtt             >  2:2198059/1‑76 (MQ=255)

                                                            cAGCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCttt            >  1:1210263/1‑76 (MQ=255)

                                                              gCGGCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTtgt          >  2:532436/1‑76 (MQ=255)

                                                                ggCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTgg        >  1:3604071/1‑76 (MQ=255)

                                                                ggCAGCAGCTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTgg        <  1:1858143/76‑1 (MQ=255)

                                                                   agcagcTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTcc     <  2:1692191/76‑1 (MQ=255)

                                                                   agcagcTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTcc     <  2:737584/76‑1 (MQ=255)

                                                                   agcagcTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTcc     >  2:914234/1‑76 (MQ=255)

                                                                      agcTCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTCCCgg  >  2:2552410/1‑76 (MQ=255)

GGAAGCCCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGCTAT‑AACGGTCA‑GGCCACGCTGAATGTGACCTTCTGACAGAACCATCGCCTCTTTGTGGTCCCGG  >  REL606/3002293‑3002438

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