BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	100.0%	227.9 / NA	92	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (50/42); total (50/42)

Only 100 of 102 total aligned reads displayed.

TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG  >  NC_000913/3554013‑3554157

tGGCGTAACAGCGTGGTTTGCTGATGAGTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                          >  1:3295970/1‑76 (MQ=255)

tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                          <  2:3898284/76‑1 (MQ=255)

 ggCTTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgc                                                                         <  1:958288/76‑2 (MQ=255)

  gCGTAGCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                        <  1:5201573/76‑3 (MQ=255)

   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCGCCCGGCGGATAAAATCCgcag                                                                       <  2:402364/76‑4 (MQ=255)

   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                       <  1:1175433/76‑4 (MQ=255)

   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                       <  1:5096802/76‑4 (MQ=255)

   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAAACCgcag                                                                       >  1:353988/1‑73 (MQ=255)

      aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                    <  1:1656419/76‑1 (MQ=255)

       aCAGCGTGGTTTGCTGATGAATCCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                   <  1:4812792/76‑1 (MQ=255)

        cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                  <  1:1976977/76‑1 (MQ=255)

        cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                  <  1:4362633/76‑1 (MQ=255)

        cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                  <  1:3787068/76‑1 (MQ=255)

         ggCGTGGTTTGGTGATGAATTCCTATTTTCGTCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                 <  2:5676351/75‑1 (MQ=255)

         aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAAACCGCAggcggcc                                                                 >  2:2298912/1‑75 (MQ=255)

         aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgccgcg                                                                 >  2:2274251/1‑76 (MQ=255)

          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCCGCACCCGGCGGATAAAATCCGCACGCGgcgc                                                                >  2:5721926/1‑76 (MQ=255)

          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATGAAATCCGCAGGCGgcgc                                                                >  2:4104492/1‑76 (MQ=255)

          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                >  2:5057189/1‑76 (MQ=255)

          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                >  1:2639754/1‑76 (MQ=255)

          gCGTGGTTCGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAAACCGCACGCGgcgc                                                                >  2:4321760/1‑76 (MQ=255)

           cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGcc                                                               <  1:955174/76‑1 (MQ=255)

            gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                              <  2:2768272/76‑1 (MQ=255)

             tGGTTTGCTGATGAATTCCTATTTTCGCCATCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGa                                                             <  1:632939/76‑1 (MQ=255)

             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGa                                                             >  1:2892822/1‑76 (MQ=255)

              ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                            <  1:2651084/76‑1 (MQ=255)

              gatttGCTCATTCATTCTTACTAGCGCCAGCAGCGACTCCAGCACCCGGCGTATAAAATCCGCAGGCGGCGCCGAc                                                            <  1:1979856/74‑1 (MQ=255)

               gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                           >  1:2553306/1‑76 (MQ=255)

               gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                           <  2:5036077/76‑1 (MQ=255)

               gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCCACg                                                           >  2:4109703/1‑76 (MQ=255)

                tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                          >  1:1157933/1‑76 (MQ=255)

                tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                          >  2:1184339/1‑76 (MQ=255)

                 ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCCGGCGGCGCCGACggg                                                         >  1:154798/1‑76 (MQ=255)

                 ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                         <  2:3660057/76‑1 (MQ=255)

                 ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                         <  1:5612681/76‑1 (MQ=255)

                  tGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCCGCGGCGACcggg                                                        >  2:4300821/1‑75 (MQ=255)

                   gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggggc                                                       >  2:4921413/1‑74 (MQ=255)

                    gtgatgcATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                      <  2:235551/75‑1 (MQ=255)

                    cTGATGAATTCCTATTTTCGCCAGCAGCTGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                      <  1:3235741/76‑1 (MQ=255)

                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCCGCGGACAAAATCCGCCGGCGGCGCCCACTGGGTATCg                                                   >  2:5074667/1‑76 (MQ=255)

                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGGGCCGGCGGGGGATCg                                                   >  1:2750722/1‑76 (MQ=255)

                       atgaATTCCTATTTCCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                   <  1:469243/76‑1 (MQ=255)

                        tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCCGGGGTATCgg                                                  >  2:4698431/1‑76 (MQ=255)

                        tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                  >  1:3754572/1‑76 (MQ=255)

                         gatttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                 <  1:1644977/73‑1 (MQ=255)

                         gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                 <  2:4085044/76‑1 (MQ=255)

                           aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                               <  2:3963088/76‑3 (MQ=255)

                              ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                          >  1:3245206/1‑76 (MQ=255)

                              ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                          >  2:3180926/1‑76 (MQ=255)

                               cgaTTTGCCCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                         <  2:4167793/74‑2 (MQ=255)

                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                         <  1:4030910/76‑2 (MQ=255)

                               cTATTTGCGCCAGCAGCGGCGCCAGCACCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                         <  1:4139832/76‑2 (MQ=255)

                                tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGt‑gccg‑ccc                                        >  2:1476032/1‑74 (MQ=255)

                                     ccGCCAGCAGCGGCTCCAGCACCCGGCTGATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCAc                                   <  1:1905043/75‑1 (MQ=255)

                                     acGCCAGCAGCGGCTCGAGCGCCCGGTGGATAACATCCGCAGGCGGCTCCGACGTGTTATCGGT‑GCCGCC‑GCTCAc                                   <  1:4424497/75‑1 (MQ=255)

                                        ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCCCACTTc                                >  2:1197616/1‑76 (MQ=255)

                                         cagcggGGGCCCCAGCCCCCGGGGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCa                               <  2:1474502/76‑1 (MQ=255)

                                         cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCa                               <  2:1157933/76‑1 (MQ=255)

                                         cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCa                               >  2:4675287/1‑76 (MQ=255)

                                         cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCa                               >  1:2015550/1‑76 (MQ=255)

                                           gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAac                             <  1:1184339/76‑1 (MQ=255)

                                            cagcCGCGCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCAGT‑GCCGCC‑GCTCACTTCAaca                            <  1:1521357/76‑1 (MQ=255)

                                             agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCCacac                           >  1:5158571/1‑76 (MQ=255)

                                               gggCTCCGCCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgc                         <  1:2450077/75‑1 (MQ=255)

                                               cGGCCCCAGGCCCCGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCAACAcgc                         <  2:4756446/76‑1 (MQ=255)

                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCCACAcgcg                        >  2:4031848/1‑76 (MQ=255)

                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                        >  1:973893/1‑76 (MQ=255)

                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                        >  1:4032361/1‑76 (MQ=255)

                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                        >  1:2620736/1‑76 (MQ=255)

                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAAGCGGCGCCGACGGGCTATCGGT‑GCCGCC‑GCTCTCTCCGACAcgcg                        >  1:430747/1‑76 (MQ=255)

                                                  cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGaa                      <  2:2604784/76‑1 (MQ=255)

                                                   tCCAGCACCGGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAg                     <  1:3246191/76‑1 (MQ=255)

                                                   tCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAg                     >  2:2708615/1‑76 (MQ=255)

                                                    ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAgg                    >  1:3415682/1‑76 (MQ=255)

                                                    ccAGCACCCGGCGGAGAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑GGCCG‑CAGCTCACCTCAACACGCGAAgg                    >  1:765461/1‑76 (MQ=255)

                                                     cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGt                   >  1:2253356/1‑76 (MQ=255)

                                                      aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTc                  >  1:2635715/1‑76 (MQ=255)

                                                      aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑CGCCGCC‑GCTCCCTTCAACACGCGAAggcc                  >  2:5465629/1‑74 (MQ=255)

                                                       gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTcc                 <  1:3948417/76‑1 (MQ=255)

                                                         aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCac               >  2:1476794/1‑75 (MQ=255)

                                                          cgcgggggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATc              <  2:2668259/74‑1 (MQ=255)

                                                          cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                        >  2:345115/1‑48 (MQ=255)

                                                          cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                        <  1:345115/50‑3 (MQ=255)

                                                          cccGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATc              <  1:837798/76‑1 (MQ=255)

                                                           ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCCACACGCGAAGGTCCATCg             >  1:3360784/1‑76 (MQ=255)

                                                            cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCgg            <  2:3754572/76‑1 (MQ=255)

                                                            cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCATccg            >  2:2795281/1‑74 (MQ=255)

                                                             ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc           <  1:5001649/76‑1 (MQ=255)

                                                             ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCATCCGc           >  1:2817253/1‑76 (MQ=255)

                                                              gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑GGCCGCC‑GCTCACTTCAACACGCGAAGGGCCATCggga          >  1:148181/1‑74 (MQ=255)

                                                              gcggATAAAAGCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCa          <  2:2708224/76‑1 (MQ=255)

                                                                ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCaaa        >  1:1384032/1‑76 (MQ=255)

                                                                  aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc      <  1:2343869/76‑1 (MQ=255)

                                                                   tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGAC‑GCTCAACTCAACACGCGAAGGTCCATCAGCAAACca     >  1:2532808/1‑76 (MQ=255)

                                                                    aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac    <  1:3791634/76‑1 (MQ=255)

                                                                    aaaaTCCGCAGGCGGAGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACATCAACACGCGAAGGACCAGCGGCAAACcac    >  2:3689105/1‑76 (MQ=255)

                                                                     aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca   <  1:2063379/76‑1 (MQ=255)

                                                                     aaaTCCGCAGGCGGCGCCGACAGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACccca   >  2:3558371/1‑76 (MQ=255)

                                                                      aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCCTCGGCAAACCAAAg  >  1:5702219/1‑76 (MQ=255)

                                                                      aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAACCCACAg  >  2:3393108/1‑76 (MQ=255)

TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAG  >  NC_000913/3554013‑3554157

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