BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,176	C→G	G182R (GGA→CGA)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,176	0	C	G	99.1%	295.2 / ‑5.6	114	G182R (GGA→CGA)	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); major base G (64/48); minor base A (0/1); total (65/49)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 114 total aligned reads displayed.

TGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑CG‑G‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTT  >  NC_000913/4013110‑4013245

tGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTaa                                                               >  2:4241801/1‑76 (MQ=255)

tGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATACTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTaa                                                               >  1:1251031/1‑76 (MQ=255)

 ggCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTT‑TCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGa                                                             >  1:1132415/1‑76 (MQ=255)

  gcgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGa                                                             <  2:1201165/76‑1 (MQ=255)

  gcgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGa                                                             >  1:3253269/1‑76 (MQ=255)

   cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGAt                                                            >  1:287028/1‑76 (MQ=255)

   cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGAt                                                            >  1:2448743/1‑76 (MQ=255)

    gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATc                                                           >  1:2188217/1‑76 (MQ=255)

    gcATGGTTTCAACGCTCTCCTGCGGAATGCAGTTATCCTGTACAACAGATAAGCCGAGAATC‑GGCG‑CGTGAAGATc                                                           >  2:3180644/1‑76 (MQ=255)

    ccgtGGTTTCAACGCTCTCCAGCGGAATACTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATc                                                           <  1:4596/73‑1 (MQ=255)

      aTGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCgg                                                         >  2:2183341/1‑76 (MQ=255)

       tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGt                                                        >  1:81966/1‑76 (MQ=255)

       tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGt                                                        <  2:1972992/76‑1 (MQ=255)

        ggTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGtt                                                       <  2:1035798/76‑1 (MQ=255)

         gTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTg                                                      >  1:1990901/1‑76 (MQ=255)

          tttCAAGTCTCTCTTGCGCAATTCTTTTCTCCTGACCACCATCTAAGCCGATAGTC‑GGCG‑CGTTAAGATCGGTTGc                                                     <  1:259877/76‑1 (MQ=255)

          tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGc                                                     <  1:4587486/76‑1 (MQ=255)

          tctCAACGCGCTCGTGGGGATTGCTGTTATCCTGCCCACCATATAAGCCGAGAATC‑GGCG‑CTTTAAGCTCGGTTGc                                                     <  1:1634384/74‑1 (MQ=255)

          gttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGc                                                     <  2:874448/75‑1 (MQ=255)

            ttgaCGCTCTCCGGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGa                                                   <  2:385710/73‑1 (MQ=255)

             cAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGat                                                  >  2:2779383/1‑76 (MQ=255)

             cAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGat                                                  <  2:3867619/76‑1 (MQ=255)

              aaCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGata                                                 >  1:3210275/1‑76 (MQ=255)

              aaCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGata                                                 >  2:3334075/1‑76 (MQ=255)

               aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGatat                                                >  1:1522054/1‑76 (MQ=255)

                 gccctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              >  2:1319581/4‑76 (MQ=255)

                 gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              >  2:3052543/1‑76 (MQ=255)

                 gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              >  1:686155/1‑76 (MQ=255)

                 gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              <  2:3081467/76‑1 (MQ=255)

                 gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              >  2:2320873/1‑76 (MQ=255)

                 gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              <  1:1923638/76‑1 (MQ=255)

                 cctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCa                                              <  1:465968/75‑1 (MQ=255)

                  ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCaa                                             <  2:2448742/76‑1 (MQ=255)

                   tctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAAc                                            >  1:1008826/1‑76 (MQ=255)

                    ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACa                                           <  1:79547/76‑1 (MQ=255)

                    ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACa                                           <  1:2054212/76‑1 (MQ=255)

                     tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAg                                          >  1:860167/1‑76 (MQ=255)

                     tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAg                                          >  2:3992034/1‑76 (MQ=255)

                     tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAg                                          >  1:4647997/1‑76 (MQ=255)

                     tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAg                                          >  1:2983544/1‑76 (MQ=255)

                       cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgt                                        >  1:4321909/1‑76 (MQ=255)

                        tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAgggt                                       >  1:3672990/1‑74 (MQ=255)

                        tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgtt                                       >  1:4165679/1‑76 (MQ=255)

                          ctgAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgttgt                                     >  1:4171292/3‑76 (MQ=255)

                          cGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgttgt                                     >  2:3032201/1‑76 (MQ=255)

                           ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgttgtt                                    >  1:2333772/1‑76 (MQ=255)

                           ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGgttgtt                                    >  1:4284157/1‑76 (MQ=255)

                            gAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGttt                                   >  2:982949/1‑76 (MQ=255)

                            gAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGttt                                   >  1:2182416/1‑76 (MQ=255)

                             aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                  <  2:3578940/76‑1 (MQ=255)

                             aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                  >  1:1218777/1‑76 (MQ=255)

                             aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                  >  1:788059/1‑76 (MQ=255)

                             aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                  <  2:4559888/76‑1 (MQ=255)

                             aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                  >  1:1189314/1‑76 (MQ=255)

                              aTGCTGTTATCCTGGCCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGtttag                                 <  2:2742166/76‑3 (MQ=255)

                                gCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                               <  2:2552464/76‑1 (MQ=255)

                                gCTGTTATCCTCACCCCCATATAAGCCGAGACTC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGCTTGTTCCGGt                               <  1:3363791/76‑1 (MQ=255)

                                 cTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTg                              >  1:3040444/1‑76 (MQ=255)

                                 cTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTg                              >  1:2650037/1‑76 (MQ=255)

                                  tggtATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                             <  2:1610803/73‑1 (MQ=255)

                                  tGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                             >  2:3692948/1‑76 (MQ=255)

                                  tGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                             >  2:1729336/1‑76 (MQ=255)

                                  tGTTATCCGGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                             <  2:4490601/76‑1 (MQ=255)

                                   gTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                            >  2:2004975/1‑76 (MQ=255)

                                   gTTATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                            >  1:1510111/1‑76 (MQ=255)

                                    ttATCCTGACCACCATATAAGCTGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                           <  1:1319582/76‑1 (MQ=255)

                                    ttATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                           <  1:2320874/76‑1 (MQ=255)

                                     tATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                          >  2:34736/1‑76 (MQ=255)

                                     tATCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                          >  1:2680732/1‑76 (MQ=255)

                                      aTCCTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                         <  2:408385/76‑1 (MQ=255)

                                       tacTGACCACCATATAAGCCGAGAATCGGG‑G‑CGTGACGATCGGTTGCGATATCAACTGGTTGTTTCGGTGAATTCa                        <  2:227823/74‑1 (MQ=255)

                                       gccTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                        <  2:2212680/75‑1 (MQ=255)

                                        ccTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGACTTccg                       >  2:4866132/1‑74 (MQ=255)

                                         cTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                      <  2:81966/76‑1 (MQ=255)

                                         cTGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                      <  1:781947/76‑1 (MQ=255)

                                          tGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                     <  1:963250/76‑1 (MQ=255)

                                          tGACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                     <  1:3630647/76‑1 (MQ=255)

                                           gACCACCATATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGa                    >  1:2087650/1‑76 (MQ=255)

                                             ccgccaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACt                  <  2:3623231/73‑1 (MQ=255)

                                             ccaccaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACt                  <  2:3966206/76‑1 (MQ=255)

                                              cacgaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                 <  2:2258968/76‑1 (MQ=255)

                                              caccaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                 <  2:1593459/76‑1 (MQ=255)

                                               accaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTg                >  1:1036844/1‑76 (MQ=255)

                                                ccaTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt               >  1:1498417/1‑76 (MQ=255)

                                                 caTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc              <  2:2109816/76‑1 (MQ=255)

                                                 caTATAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc              >  1:3978167/1‑76 (MQ=255)

                                                  atatAAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTcc             >  2:4129753/1‑75 (MQ=255)

                                                    ataAGCCGAGAATC‑GGCG‑CGTTAATATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGcc           <  1:2776856/76‑1 (MQ=255)

                                                     taAGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGccc          >  1:3416409/1‑76 (MQ=255)

                                                      aaGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGGATTCAGCGACTTGTCGCCCg         >  1:2476662/1‑76 (MQ=255)

                                                      aaGCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg         <  1:383956/76‑1 (MQ=255)

                                                        gCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTc       >  1:4593291/1‑76 (MQ=255)

                                                        gCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTc       <  1:1133434/76‑1 (MQ=255)

                                                        gCCGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTc       <  2:4647997/76‑1 (MQ=255)

                                                         ccGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCa      >  1:1831832/1‑76 (MQ=255)

                                                         ccGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGCCa      >  2:4715356/1‑76 (MQ=255)

                                                          caagaATC‑AG‑GACGTTAAGATCGGTTGCTCTATCAACAGGTATTTTCGGGTACTTGAGCGGCTTGTCGCCCGTCAg     <  1:1641796/74‑1 (MQ=255)

                                                          cGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAg     >  1:3293666/1‑76 (MQ=255)

                                                          cGAGAATC‑GGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAg     >  1:4079474/1‑76 (MQ=255)

                                                             gaATC‑CGCG‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCCGttt  >  2:429014/1‑76 (MQ=255)

TGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑CG‑G‑CGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTT  >  NC_000913/4013110‑4013245

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